Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA23.2% 158.1 / 21.1 69intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (46/7);  new base A (0/16);  total (46/23)
Fisher's exact test for biased strand distribution p-value = 1.28e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752991‑753264
                                                                                                                                      |                                                                                                                                           
ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                        >  1:2881273/1‑140 (MQ=255)
ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                        >  1:2309272/1‑140 (MQ=255)
ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                                        >  1:1173020/1‑140 (MQ=255)
    aTGCGACTACCATGAATTTTTTAATTAAAATTTTTGGTTTGTTAATAATTTGAGCTTTTTTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAACCac                                                                                                                                    <  2:39902/140‑2 (MQ=255)
     tGCGATTACAATGAATTTTTTAATTAAAAGTATTGGTTTGTTGATAATTTGAGCTTTTTTTTTTTTTTTAAATATCTATATAGTTTTTTTAATGGTTTTTATTTTTAAATTTTTTTTTGTTTAGGCATATAAAAATCAAc                                                                                                                                   <  2:2741215/140‑1 (MQ=255)
      gCGACTACCATGAATTTTTTAATTGAAAGTTTTGGTTTGCTAATAATTTGAGTTTTTTTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAAcc                                                                                                                                  <  2:1029197/140‑1 (MQ=255)
       cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                            >  1:368155/1‑129 (MQ=255)
       cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCAtata                                                                                                                                            <  2:368155/129‑1 (MQ=255)
       cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                                 <  2:2881273/140‑1 (MQ=255)
       cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                                 >  2:3060592/1‑140 (MQ=255)
         aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                               >  2:463279/1‑140 (MQ=255)
           tACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCa                                                                                                                             >  1:2204599/1‑140 (MQ=255)
                   atttctaaaTTAAAAGTATTGGGTTGCTGATATTTTGAGTCTTTCTTTTTTTTTTAAATATCTATATAGGTTTGTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCTTATGAAcg                                                                                                                     <  2:2447217/133‑1 (MQ=255)
                     tttttAATTAAAATTTTTGGTTTGCTAATAATTATGGCTTTTTTTTTCTTTTAAAATATCTATATAGGTCTTTTATGGTTTTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAAACACCCCGCCATATGAAcggc                                                                                                                   <  2:606772/140‑1 (MQ=255)
                       tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTGGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                                                                 >  1:1528979/1‑140 (MQ=255)
                        ttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                                >  2:1617970/1‑140 (MQ=255)
                          aaTTAAAAGTATTGGGTTGCTGATATTTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGtt                                                                                                              <  2:2204599/140‑1 (MQ=255)
                             taaaaGTTTTGGTTTGCTGATATTTTGAGTTTTTTTTTTTTTTTAAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTAAATTTTTTTGTTTTTGGGCAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTaaa                                                                                                           <  2:1528979/138‑1 (MQ=255)
                               aatttATTGGTTTGTTGATAATTTGGCGTTTTCTTTTCTTTTTAATTCTCTAAAGAGGTGTTTAAAGGGTTTTTATTTTTTAAATTTTTTTCTTTTAGGTATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAat                                                                                                         <  2:1999085/136‑1 (MQ=255)
                                          ttGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGTTTTTTTAATGGATTTTTTTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                              <  2:2077461/140‑1 (MQ=255)
                                             cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                           >  2:904602/1‑140 (MQ=255)
                                                atatTTTGAGTTTTTCTTTTTTTTTAAAATATCTATATGGGTTTTTTAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                        <  2:2436312/140‑1 (MQ=255)
                                                  aaTTTGAGCTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                                      <  1:3042933/140‑1 (MQ=255)
                                                       gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                 >  1:2948543/1‑140 (MQ=255)
                                                       gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                 >  2:2146181/1‑140 (MQ=255)
                                                         gCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                               <  1:1652819/140‑1 (MQ=255)
                                                             ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGc                                                                           >  2:2557596/1‑140 (MQ=255)
                                                             ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGc                                                                           >  1:3133487/1‑140 (MQ=255)
                                                                    ttttttAAATATATATTTAGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAATGCTTGatat                                                                    <  2:2793008/139‑1 (MQ=255)
                                                                       tttAAATATCTCTATAGGTCTGTTAATGAATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGCTGAAAATATTTAGAACTTTGACAGAAATCAATAACGCTTTCTATCgg                                                                 >  2:2900310/1‑139 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTTCAACTTAGCAATCAACCATTAAAGCTTGGTATCGc                                                                 >  2:1122954/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  2:2715936/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:2434816/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  2:2169324/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:2120382/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:1455530/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:94494/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:724569/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:1435768/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:1016460/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACCTAACAATCAACCATTAATGCTTGATATCGc                                                                 >  1:2793008/1‑140 (MQ=255)
                                                                       tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                 >  1:683984/1‑140 (MQ=255)
                                                                                     taGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                   <  1:463279/140‑1 (MQ=255)
                                                                                        gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                >  1:2222832/1‑140 (MQ=255)
                                                                                        gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                >  1:1965827/1‑140 (MQ=255)
                                                                                        gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGCTTAAATATTTACAACTTATCAAAAAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGtttt                                                >  2:1991143/1‑140 (MQ=255)
                                                                                          ttGTTAATGGTTTTTTTTTTTACATTTTTTTGTGTTTAGCCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                              <  2:1435768/139‑1 (MQ=255)
                                                                                          ttGTTAATCGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                              <  1:2900310/139‑1 (MQ=255)
                                                                                           tGTTAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             <  1:510413/140‑1 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  1:808643/1‑140 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  1:1249129/1‑140 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGAAATCAACCATTAAAGCTTGATATCGCTTTTTAAAACGCGTTTTTCa                                             >  2:1106602/1‑140 (MQ=255)
                                                                                           tGTTAATGGATTTTATTTTTACAAATTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                             >  1:3111211/1‑140 (MQ=255)
                                                                                            tttatgtgATTTTATTTTTACATTTTTTTGTGTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:683984/133‑1 (MQ=255)
                                                                                            tttAATGTTTTTTTTTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:808643/139‑1 (MQ=255)
                                                                                            tttAATGTTTTTTATTTTAAAAATTTTTTGTTTTTAGTCAAAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:1057536/139‑1 (MQ=255)
                                                                                            tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  2:94494/139‑1 (MQ=255)
                                                                                            tttAATGGTTTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:1576474/139‑1 (MQ=255)
                                                                                            gTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:2557596/140‑1 (MQ=255)
                                                                                            gTTAATGGATTTTTTTTTTAAAATTTTTTGTGTTTAGCCAAAGAAAAACCACCCCGCCATATGAACGGCTGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                            <  1:620727/140‑1 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:1807031/1‑140 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:2556003/1‑140 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:2744926/1‑140 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  2:268745/1‑140 (MQ=255)
                                                                                             ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                           >  1:612160/1‑140 (MQ=255)
                                                                                               attGGATTTTTTTTTTACATTTTTTTGTGTTTGGGCATATAAAAATCAGCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  1:1991143/138‑1 (MQ=255)
                                                                                               aaTGGATTTTATTTTTACATTTTTTTGTTTTTAGAAATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  1:1122954/140‑1 (MQ=255)
                                                                                               aaTGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                         <  1:1122943/140‑1 (MQ=255)
                                                                                                  ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAAATTAACAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTg                                      >  2:1576474/1‑140 (MQ=255)
                                                                                                  ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCTTATCCtt                                      >  1:1790853/1‑139 (MQ=255)
                                                                                                       ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtataa                                 >  2:694060/1‑139 (MQ=255)
                                                                                                       ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                                 >  1:2631264/1‑140 (MQ=255)
                                                                                                       ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                                 >  2:1884543/1‑140 (MQ=255)
                                                                                                       ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                                 >  1:1057536/1‑140 (MQ=255)
                                                                                                        tattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACa                                >  2:1681626/1‑140 (MQ=255)
                                                                                                                 atttttttttttttAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                       <  2:1807031/138‑1 (MQ=255)
                                                                                                                   ttttttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa                     <  1:1617970/140‑1 (MQ=255)
                                                                                                                        tgttTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  1:3060592/140‑1 (MQ=255)
                                                                                                                        tgttTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  2:2120382/140‑1 (MQ=255)
                                                                                                                        tgttTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  2:1249129/140‑1 (MQ=255)
                                                                                                                        tgttTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  2:2222832/140‑1 (MQ=255)
                                                                                                                        tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  1:2715936/140‑1 (MQ=255)
                                                                                                                        tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc                <  2:724569/140‑1 (MQ=255)
                                                                                                                         ttttttAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa               <  2:943124/137‑1 (MQ=255)
                                                                                                                              tagaaataaAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt          <  2:2556003/135‑1 (MQ=255)
                                                                                                                              tAGGAATAAAAAAATCATCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt          <  2:1965827/140‑1 (MQ=255)
                                                                                                                               aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc         >  1:2000016/1‑140 (MQ=255)
                                                                                                                                ggCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa        <  2:1173020/140‑1 (MQ=255)
                                                                                                                                     ataAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCg   >  2:752901/1‑140 (MQ=255)
                                                                                                                                      taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt  >  1:230943/1‑140 (MQ=255)
                                                                                                                                      taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt  >  1:1004148/1‑140 (MQ=255)
                                                                                                                                      taAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAACAACCATTAACGCTTGGTATCGCTTTTAAAGTCGCGGTTTTCATATCCCGTATACAGCTTAAGCGGAGGGGCAATCTTTATACCGt  >  2:207904/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                           
ATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGT  >  NC_000913/752991‑753264

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: