Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3580TG44.6% 83.0 / 23.2 67T66P (ACC→CCC) yafYlipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (4/32);  new base G (29/0);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.83e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266214‑266506
                                                                                                                                                   |                                                                                                                                                      
cctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGAc                                                                                                                                                  <  1:799178/150‑1 (MQ=255)
   ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGa                                                                                                                                              <  1:566832/151‑1 (MQ=255)
   ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACAGCAAGGGGGGGcc                                                                                                                                                 >  2:821058/1‑148 (MQ=255)
    gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGAt                                                                                                                                             <  2:747824/151‑1 (MQ=255)
     tATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAggggggggcgg                                                                                                                                                >  2:305706/1‑143 (MQ=255)
      aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATa                                                                                                                                            <  1:1312052/150‑1 (MQ=255)
      aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGGAAAAACCGCAAAGGGGGGCCCGaga                                                                                                                                            >  1:561292/1‑148 (MQ=255)
       ttCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAGGGAAAAACCTCAAAGGGGGGCCCgggga                                                                                                                                            >  1:378043/1‑147 (MQ=255)
         ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGATGTAAAAACTCCAAAGGGGGGGCGGAGaa                                                                                                                                           >  1:1297816/1‑148 (MQ=255)
          cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACTGCAAGGGGGGACCGGATa                                                                                                                                            >  1:763100/1‑146 (MQ=255)
               tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                                   <  2:76439/150‑1 (MQ=255)
                gATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAAcc                                                                                                                                 <  1:453990/151‑1 (MQ=255)
                     tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtgtg                                                                                                                             <  1:1169239/150‑1 (MQ=255)
                      aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGATAAGAAGGCAACCTGTGAc                                                                                                                           >  2:220240/1‑151 (MQ=255)
                       tCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACAGCCAAGGGGGGCCGGGTAAAAAGTCAAACCGTGGcc                                                                                                                          >  1:1057395/1‑151 (MQ=255)
                        cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAAGGGGGGCCGGGGAAGTAGGCAACCTGTGAccc                                                                                                                         >  2:817747/1‑150 (MQ=255)
                        cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGAAAAAACCACAAGGGGGGGCCGAGTAAGTGGGGAAACTGTGagc                                                                                                                          >  2:480287/1‑148 (MQ=255)
                             aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttc                                                                                                                      <  2:1224372/149‑1 (MQ=255)
                                        aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGGCCGGATAAGTAGTCAACCTGTGACCATTTTTCGCTGAcatc                                                                                                         >  2:596051/1‑150 (MQ=255)
                                        aTTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGGAAGAAGTGAAACTGTGAACAttcttttctggca                                                                                                            >  2:45942/1‑140 (MQ=255)
                                          tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAACGGAAAACTTCCTCGGGGGGGCCGGGAAAGAGTCAAACCCTGTGCACACTTCTCTGAcatcctt                                                                                                      >  1:1072962/1‑148 (MQ=255)
                                           catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTTAAAACTGCAAAGGGGGGCCGGGGAAATAGTCAAACCGTGACCaattttttct                                                                                                                >  1:267918/1‑132 (MQ=255)
                                            atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACAAAAAAAACTCAAGGGGGGGCCGGGTAAATAGTCAAACTGTGACCCTTCTTCGCTGAcataaaccc                                                                                                    >  2:743448/1‑145 (MQ=255)
                                                 tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCAGCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCtggt                                                                                                 <  1:637968/149‑1 (MQ=255)
                                                         cGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACcac                                                                                       <  2:71097/151‑1 (MQ=255)
                                                           gTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGCCGGGGAAAAAGTCAACCCGTGACCCTTCTTCTCTGACATCATCCTGGGGGGTGACcac                                                                                       >  2:931008/1‑149 (MQ=255)
                                                             tttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACCGAAAAACCGCCAGGGGGGGCCGGGTAAAAAGGTAAACCGTGACCATTCTTTCCCGACATCATCCTTGGGGCTGGcccc                                                                                       >  2:204996/1‑145 (MQ=255)
                                                               tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAAATGCAAGGGGGGACCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGAcc                                                                                         >  2:410753/1‑143 (MQ=21)
                                                                ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGTACTAAAACCCGCAAGGGGGGAACCGGAAAAAAATCAAACTGTTAACCTTTTTTTCTTTCATCATCCtgctggtctgccccaaccctt                                                                                >  2:74457/1‑136 (MQ=255)
                                                                ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGTCTGACCACATCACt                                                                                 <  1:1289039/150‑1 (MQ=255)
                                                                ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACt                                                                                 <  2:258286/150‑1 (MQ=255)
                                                                 tCCTACATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTg                                                                                <  1:76264/150‑1 (MQ=255)
                                                                  ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                              <  2:1258681/151‑1 (MQ=255)
                                                                     gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAGTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaatc                                                                            <  1:177423/150‑3 (MQ=255)
                                                                      cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgc                                                                          <  1:495806/151‑1 (MQ=255)
                                                                       atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGGAAAAACTGCAAAGAGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGACCCCATCCCTCAAAcac                                                                          >  1:1201521/1‑148 (MQ=18)
                                                                        tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGGAAAGTAGTCAACCTGTGGCCCTTCTTTGCCTACCTCATTCTTGTGGGTGACCACATCCCTGaaaaaggg                                                                        >  2:131218/1‑146 (MQ=255)
                                                                           ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTgg                                                                     <  2:1096822/151‑1 (MQ=255)
                                                                           ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGGCGGGGAAAAAGTCAACCCGTGGACATTCTTTGCTGACATCATCCCTGGGGGTGACCAAAACACTGAAAAGCGgggg                                                                     >  2:1297641/1‑151 (MQ=255)
                                                                             aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACAGCAAAGGGGGAACGGGTAAGAAAGCAACCTTTGAACATTTTTCTGCTGCAATCATCTTGTGGCTGACCACCTTACTAAAAAACGgggggg                                                                   >  1:898842/1‑151 (MQ=255)
                                                                               gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCCCACGGGGGGGGCGGATAAGAAGTCAACCCGTGAACCTTCTTCGCTGTCATCATCCTTGTGGCTGGCCACACCACTGAAAAGCGgggggggg                                                                 >  1:971876/1‑149 (MQ=255)
                                                                               gAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGCGGTAAAAACTGAAAGTAGGGACCCGGAAAGCGATAAACCCCTCACCATTTTTCTCTTCGACAATCCCGGGGGCTGGCCCCATCCCTGAAAACCGGTggggg                                                                  >  2:1110177/1‑149 (MQ=255)
                                                                                aGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                                <  2:125851/151‑1 (MQ=255)
                                                                                 ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAg                                                               <  2:852972/151‑1 (MQ=255)
                                                                                 ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGACCGGAAAAAAAGTCAAACTGTGACCATTCTTCCCTGACAACATCCCGGTGGCTGACCACATCACTGAAACGCGgggggggg                                                                 >  2:907149/1‑147 (MQ=255)
                                                                                    agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTc                                                            <  2:139417/151‑1 (MQ=255)
                                                                                      tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACTGCAAGGGGGGGCCGGGGAAGTAGGCAAACCGTGACCCTTCTTCGCCGACATCATCCCGGGGGGTGGGCCCATCACTGAAACGCGggggggggggggcc                                                           >  1:1121858/1‑142 (MQ=255)
                                                                                        gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAAGGGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGACCACATCACTTAAAAGCGGGGGGGTGGGGGCCtt                                                         >  2:514947/1‑150 (MQ=18)
                                                                                          aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                       <  1:200497/150‑1 (MQ=255)
                                                                                          aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTaa                                                       <  2:635894/150‑1 (MQ=255)
                                                                                               ccATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAGGGGGGGCCGGGTAAGAAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGCGGGGGGAGGGCCTTTAAttt                                                   >  2:26040/1‑149 (MQ=18)
                                                                                                 aTCCTCGAAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTg                                               <  2:1087161/151‑1 (MQ=255)
                                                                                                    cTCGCAGCCATCGGGCGACTGGCGGGTACGTCAAACCTGCAAGGTGGGGCCGGATAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCTGGTGGGTGACCACAACACTGAAACGCGGGGGGGGGGGGGCCTTTAATTTACt                                                >  2:504208/1‑147 (MQ=255)
                                                                                                        cAGCCATCGGGCGACTGGCGGGGACGTAAAACCCGCAAGGGGGGAGCGGATAAGAAGTCAACCCGTGACCCTTCTTCGCTTACCTCCTCCCGGTGGGTGGCCACATCACTGAAACGCGGGGGGGGGGGGGGCTTTAATTTACTGATgggc                                         >  1:849572/1‑148 (MQ=255)
                                                                                                          gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                                      <  1:1197946/151‑1 (MQ=255)
                                                                                                                gggCGACTGGCGGGTACGTAAAAACCGCAAAGGGGGGCCGGGTAAAGAGGCAAACCGTGGCCATTCTTTGCTGGCATCATCCCGGGGGGTGGCCCCAACACTTAAAAGGGGGGGGGGGGGGGGCTTTAATTTTCTGATgggccttagg                                   >  1:1265569/1‑140 (MQ=255)
                                                                                                                 ggCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat                               <  1:408560/151‑1 (MQ=255)
                                                                                                                     aCTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACt                           <  2:1130843/151‑1 (MQ=255)
                                                                                                                       tGGCGGGTACGTAAAAACCGCAAGGGGGGACCCGGTAAGTAGGTAACCTGCGGCCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGGGGTGGGGTGTGGTCCTTAAATTTACTGATCGCCTTCCGGtcatc                              >  1:797539/1‑146 (MQ=18)
                                                                                                                        ggCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt                        <  2:428753/151‑1 (MQ=255)
                                                                                                                         gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                       <  1:1254423/151‑1 (MQ=255)
                                                                                                                         gCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTg                       <  1:784059/151‑1 (MQ=255)
                                                                                                                             ggTACGTAAAACCTGCAAGGTGTGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa                    <  1:707042/151‑1 (MQ=255)
                                                                                                                                      aaaaCCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGAGCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCAt              <  1:1063727/150‑1 (MQ=255)
                                                                                                                                             gCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATg        >  2:1216246/1‑149 (MQ=255)
                                                                                                                                                aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaaa   <  1:907992/151‑1 (MQ=255)
                                                                                                                                                aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa    >  1:1332301/1‑150 (MQ=255)
                                                                                                                                                  gTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  2:560895/1‑150 (MQ=255)
                                                                                                                                                   tGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  2:254591/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266214‑266506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: