Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,051 A→T 22.2% intergenic (+30/‑166) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0510AT22.2% 203.1 / 14.8 96intergenic (+30/‑166)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (47/22);  new base T (2/20);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 9.69e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGA  >  CP009273/160902‑161198
                                                                                                                                                     |                                                                                                                                                    
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      >  1:4432283/1‑150 (MQ=255)
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      >  1:2101325/1‑150 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                     <  1:3413368/150‑1 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                     <  2:4483115/150‑1 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                     <  1:2564582/150‑1 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                                  <  1:4090384/150‑1 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                                  <  1:4090630/150‑1 (MQ=255)
     gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                 <  1:2045899/150‑1 (MQ=255)
      aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                >  1:3078803/1‑150 (MQ=255)
       aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                               >  2:122180/1‑150 (MQ=255)
         ggAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                             <  2:248052/150‑1 (MQ=255)
           aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                           >  1:1574259/1‑150 (MQ=255)
             cGTACCGTGAGGTGCACAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         >  2:4011894/1‑150 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         <  2:1066209/150‑1 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         >  2:3074465/1‑150 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         >  1:323212/1‑150 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         >  2:3074461/1‑150 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         <  1:3005002/150‑1 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        >  2:658367/1‑150 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        >  2:1849262/1‑150 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        >  1:4126294/1‑150 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        >  1:4126317/1‑150 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        >  1:4126449/1‑150 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        >  1:836916/1‑150 (MQ=255)
                   gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTc                                                                                                                                   >  1:1532806/1‑150 (MQ=255)
                   gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                   >  1:2808784/1‑150 (MQ=255)
                   gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                   >  1:1532781/1‑150 (MQ=255)
                     gAGTTGAAAAAAGTGATGAAAGGGCGTTATCCAAAACAGGTCTGCCCGGACGACCCGGCAAATAGTGCACAGACGCGACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCTCCTTGTGGGTGAGGGCCGGGGTTAGGGCATCAg                                                                                                                                 <  2:882629/150‑1 (MQ=255)
                     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAg                                                                                                                                 >  1:3984749/1‑150 (MQ=255)
                     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                 >  1:3984708/1‑150 (MQ=255)
                     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                 <  2:3923247/150‑1 (MQ=255)
                       tgTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                               >  2:3730097/2‑150 (MQ=255)
                        gTGCAATAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGTCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                              >  1:2180885/1‑150 (MQ=255)
                            aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                          >  1:704361/1‑150 (MQ=255)
                             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                         >  1:3048403/1‑150 (MQ=255)
                                  agaTGACAGGGGTTTAGCCAAACCATGTGGGGCCGGCCGCCCCGTCAATTACGGCCCCGACGCGCGGGACTAAAAAGTTTTCTTATGGTTGGTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                    <  2:3540135/150‑1 (MQ=255)
                                  agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggac                                                                                                                    >  1:2211188/1‑146 (MQ=255)
                                   gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTGAGGTGGGTTTTCTCACTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCGGCGCGCACggggtcg                                                                                                                   >  2:3227884/1‑144 (MQ=255)
                                   gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttggcc                                                                                                                   >  2:3227829/1‑146 (MQ=255)
                                    tttaacggggGTTATCCAAAACTTGTGGGGCCGGACGACCGGGAAAATACTGCCCCGCCGCGCGGGACTAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcaccc                                                                                                                  <  2:3660952/144‑1 (MQ=255)
                                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                  >  2:736806/1‑150 (MQ=255)
                                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccc                                                                                                                  >  2:2609658/1‑150 (MQ=255)
                                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaacc                                                                                                                  >  2:664538/1‑150 (MQ=255)
                                     tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccca                                                                                                                 >  2:1137097/1‑149 (MQ=255)
                                      gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccgc                                                                                                                >  2:3314716/1‑148 (MQ=255)
                                      gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtcacccac                                                                                                                >  2:2263109/1‑148 (MQ=255)
                                        atgggCTTTATCCCAAACAGGTCTGGCCGGACCACCCGGAAAATACGGCACCGACGCGACGGACTAAAAAGTATTCGTAAGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccctcac                                                                                                              <  2:1049168/148‑1 (MQ=37)
                                           ggCGTTAGCCCACACTTGTCTGGCCGGACGACCGGGAAAATACTGCACCGACGCTCCGTACTAAAATGTTTTCGTACGTTTGGTTTTCTCCCTCTCCTTGGGGGTGTGGGCCGGGGTTTGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                           <  2:4361021/150‑1 (MQ=255)
                                             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAGCCACACCCaaa                                                                                                         >  2:925400/1‑150 (MQ=38)
                                             cGTTACCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACTAAAAATTATTCGTAAGGTTGTTTTTCTCCCTCTCCCGGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                         <  2:92818/150‑1 (MQ=255)
                                               ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACCCCCCAAca                                                                                                       >  2:3206047/1‑149 (MQ=255)
                                                  tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                    <  2:2928498/150‑1 (MQ=255)
                                                  tCCAAAACTTGTGTGGCCGTACGACCCGGAAATTAGTCCACCGACGCGACGTACTAAAATTTATTGTTATGTTTGTTTTTCTCCCTCGCCCTGGGGTTTAGGGCCGGTGTGCGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                    <  2:303423/150‑1 (MQ=255)
                                                  ccccAAACTTTTCTGGCCGGACTACCCGGCAAATCCTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                    <  2:303428/149‑1 (MQ=38)
                                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCGAACCCAATccc                                                                                                >  2:320079/1‑150 (MQ=38)
                                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCACTccc                                                                                                >  2:2518406/1‑150 (MQ=255)
                                                       aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAAACCACAcccg                                                                                               >  2:948922/1‑149 (MQ=38)
                                                         cttGTCTGGCCGCACGCCCCGGCAAATACTGCACCTACGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                             <  2:4389777/148‑1 (MQ=38)
                                                         cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGACCCTCACCCTAAACCACAcccacc                                                                                             >  2:2884921/1‑147 (MQ=38)
                                                         cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCCCACCCCAAACCACAcccacc                                                                                             >  2:2466792/1‑147 (MQ=37)
                                                           tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGCGcccgccac                                                                                           >  2:3376635/1‑145 (MQ=38)
                                                                ggcggcacgacCCGCAAAATACTGCCCCGCCGCGACGCACTAAAATTTATTCTTATGTTGGTTTTTCTCCATCTCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                      <  2:836916/144‑1 (MQ=255)
                                                                ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                      <  2:4126317/150‑1 (MQ=255)
                                                                ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                      <  2:4126294/150‑1 (MQ=255)
                                                                ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                      <  2:4126449/150‑1 (MQ=255)
                                                                  ccGGACGACCCGGCAAATACTGCCCCGACGCGAGGGACTAAAAATTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGTGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                    <  2:1946230/150‑1 (MQ=38)
                                                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACACCCTAACCCAAAGCCGCAAGGGAGAgggg                                                                                  >  2:2846760/1‑150 (MQ=38)
                                                                      acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGACCCCCCCCCAAAACCACAACCCCAAGGGAgggggggc                                                                                >  1:1484618/1‑148 (MQ=25)
                                                                       cgacCCGGAAAATACTGCCCCGACGCGCCGGACTAAAATGTATTCGTATGGTGGTTTTTCCCCCTCCCCCGGTGGGATCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                               <  1:3853670/150‑1 (MQ=255)
                                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCGAACCCCCTCCCCCAAGGGGGCggggggcgctg                                                                           >  2:4278955/1‑143 (MQ=37)
                                                                            cgcaaaaTACTGCCCCGACGCGCCGGACTAAAAAGTATACTTCAGGTTGTTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                          <  2:3137300/146‑1 (MQ=37)
                                                                                aaaTACTCCCCCGCCGCGACGGACGAAAAATTATTCGTACGGTTGGTTTCCCCCCACTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                     <  2:1375861/150‑1 (MQ=255)
                                                                                     cTGCACCGACGCTCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCATCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                <  2:331003/150‑1 (MQ=255)
                                                                                      ggCACCGAGGCGACGGACGAAAAAGATTACGTCTGTTTGGTTTAACCACACGATGTGTGGTAGTGGGCTGGGGATTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                               <  2:1200551/149‑1 (MQ=255)
                                                                                        cACCGCCGCGAGGCACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTCTCGCGGTGGGAGGGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                             <  1:1765755/150‑1 (MQ=255)
                                                                                         aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                            <  2:4252329/150‑1 (MQ=255)
                                                                                              cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAGCCCCACCCCAAACCACAACCACCAGGGGGAGGGGGGCGGACGAGCCCAAAAGTTgaggga                                                       >  1:4549361/1‑150 (MQ=18)
                                                                                                cgacgGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                     <  2:2570418/150‑1 (MQ=39)
                                                                                                    ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTTGACCCACACCCCAAACCCCAACCACAAGGGGGAGGGGACCGCTCGAGCACAAAAATTTATAGATAtatt                                                 >  2:1767203/1‑150 (MQ=21)
                                                                                                     gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACACCCAAACCCACAACCACAATGGCGAGGGGACCGAACGAGCACAAAAATTGAGTTATAtatta                                                >  1:478596/1‑149 (MQ=21)
                                                                                                     gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACACGCACCCAAAACCACACCCACCAAGGCGAGGGGGCCCAGCGAGCACAAAATTTGAGAGATttcttc                                                >  2:4431998/1‑150 (MQ=21)
                                                                                                     gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCAAACCCAAAAACACACCCACAAATGAGAGGGGCCCCCACGAGCACACAAATTGATAGatatataa                                                >  1:478595/1‑146 (MQ=12)
                                                                                                              tttttCTTAGGGTTGTTTTTCCCCCCCTCCCTGTGGTTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                       <  1:1326392/147‑1 (MQ=255)
                                                                                                               tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGAACCCAATACCACAATGGGGGGGGGACCCATAGAGGCCCAAATTTTTGAGATATATTATTCTGTCtgg                                      >  2:3025940/1‑148 (MQ=16)
                                                                                                                aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                     <  2:3828904/150‑1 (MQ=255)
                                                                                                                    ttagggttggttTTCCCCCTCGCCGTGTGGGTGGGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                 <  1:4335182/146‑1 (MQ=255)
                                                                                                                    gTATGGTTGTTTTTCCCCCTCCCCCGGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                 <  1:4335180/150‑1 (MQ=255)
                                                                                                                         gttggttTCCTCCCGGCCCGTGTGGGTGCGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                            <  1:1801445/150‑1 (MQ=255)
                                                                                                                           tggttTTCCCCCCCTCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                          <  2:323212/150‑1 (MQ=255)
                                                                                                                                   tccctctccctGGGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                  <  2:1987483/150‑1 (MQ=255)
                                                                                                                                    ccctcgcccgGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                 <  2:51845/150‑1 (MQ=255)
                                                                                                                                     cagctccgtGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                <  1:658367/147‑1 (MQ=255)
                                                                                                                                        ctccctGTGGTCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt             <  2:3606276/150‑1 (MQ=255)
                                                                                                                                         tccctGTGGGCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt            <  1:1921504/150‑1 (MQ=255)
                                                                                                                                           cgtgtgGGAGTGGCCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt          <  1:3954712/148‑1 (MQ=255)
                                                                                                                                           cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt          <  1:2603525/150‑1 (MQ=255)
                                                                                                                                           cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt          <  2:3470301/150‑1 (MQ=255)
                                                                                                                                            ctGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg         <  2:621394/150‑1 (MQ=255)
                                                                                                                                                 ggAGAGGGCGGGGTTCTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt    <  2:1977252/150‑1 (MQ=255)
                                                                                                                                                   agagGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa  <  2:1452387/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGA  >  CP009273/160902‑161198

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