Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.4% intergenic (+47/‑149) 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
*NC_000913164,5810AC17.4% 330.7 / 31.2 193intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (50/106);  new base C (2/35);  total (52/141)
Fisher's exact test for biased strand distribution p-value = 7.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728
                                                                                                                                                     |                                                                                                                                                   
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                     >  2:1231520/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  <  2:1395528/150‑1 (MQ=255)
     tggTGCAAAATGAGATTAATGGGCGTTATCCCAAACTTGTCTGGCGGGACGACCCGGAAAATAGTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAgc                                                                                                                                                <  2:3887465/149‑1 (MQ=38)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgc                                                                                                                                              >  2:8439651/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                             >  1:6210334/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgcc                                                                                                                                           >  1:3322407/1‑149 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                           >  1:4899231/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                         >  1:7907832/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                         >  1:8016109/1‑150 (MQ=255)
                gtgttGAACGGGCTTTATCCCAAACTTGTCGGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTTTTCGTAGGGTTGGTTTTCGCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                     <  2:8845667/146‑1 (MQ=37)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGtgggc                                                                                                                                    >  1:7323160/1‑146 (MQ=255)
                  taTGAAAGGGCTTAATCCCAAACAGGTCGGGCCGGACGACCGGGAAAATAATGCACCGCCGGGACGGAGGAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                   <  2:1711987/149‑1 (MQ=25)
                  gttGAAGGGGGGTTATCCAAAACTTTTGGGGCGGGCCGCCCGGGCAAATCCTGCACCGCGGGGACGGACGAAAATGTTTTGGTTTGGTTGGTTTTCGCCCTCTCCGTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTtcacc                                                                                                                                   <  2:166374/148‑1 (MQ=255)
                  gttGAACGGTCGTTATCCCAAACATGTCGGGCCGGACGACCCGGCAAATCCGCCACCGACGCGACGGACGAAAATGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTTTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                   <  2:632726/148‑1 (MQ=37)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTgcccccac                                                                                                                                >  1:5163512/1‑148 (MQ=255)
                            cGTTATCCCAACCATGTGTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                         <  2:1561830/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                       <  2:6102494/150‑1 (MQ=255)
                               tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCACCCACCCCCGAAccc                                                                                                                      >  2:9005520/1‑150 (MQ=38)
                               tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCCCCCAAAacc                                                                                                                      >  2:2837016/1‑150 (MQ=38)
                                  cccAACCATGTCGGGCGGGACGCCCCGGAAAATCCTGCACCGCCGCGACGGACTAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                   <  2:8716698/150‑1 (MQ=25)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCcact                                                                                                                   >  1:8335855/1‑150 (MQ=255)
                                   ccAAACATGTCTGGCCGGACGACCCGTCAAATACGGCACCGACGCTAGGCACTAAAATTTTTTCTTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTTTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                  <  2:436757/150‑1 (MQ=25)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCGCACCCcctc                                                                                                                  >  2:7304776/1‑150 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAaccacac                                                                                                                  >  2:5808252/1‑146 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGCACCCCCCCCCCAAcccacac                                                                                                                  >  1:6561472/1‑146 (MQ=38)
                                    cAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                 <  2:6687417/150‑1 (MQ=255)
                                    cAAACTTGTCGGGCCGGACGCCCCGGAAAATAGTGCACCGACGCGACGGACTAAAATGTTTTCGTTTGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                 <  2:4815402/150‑1 (MQ=37)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACAACCAAAcccacaac                                                                                                                 >  2:5626118/1‑145 (MQ=38)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCCAACCCACAccc                                                                                                                >  1:1734140/1‑150 (MQ=37)
                                         cggTCTGGCCGGACGCCCCGGCAAATCCTGCCCGGCCGCGCCGGACTAAAACGTTTTGGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTTCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:6484700/148‑1 (MQ=255)
                                         cggTCGGGCCGCCCGACCCGGAAAAAACTCCCCCGCCGCCACGGACTAAAAATTATCCTTAGGTTTGTTTTACTCCCTCACCGTGGGGGAGGGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                            <  2:3272192/148‑1 (MQ=255)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACACCCCAACCCACACCCCCCAggg                                                                                                         >  2:2705220/1‑150 (MQ=37)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAACCCCAACCCACAAggg                                                                                                         >  1:1256125/1‑150 (MQ=37)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAAACCACAACCCCAAGGGcgag                                                                                                     >  2:6611028/1‑150 (MQ=37)
                                                  cGGACGCCCCGCCAAATCCTGCCCCGCCGCGCGGGCCGAAAACTTTTTCGTAGGGTGGGTGTTCCCCCCCGCCCGGGGGGAGGGGCCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                   <  2:2493825/150‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCAAAACCACTCCCTCCAGGGCGTggg                                                                                                   >  2:5400297/1‑150 (MQ=37)
                                                    gacgccCGGGAAACTACGGCCCCGCCGCGACGGCCGAAAAAGTAGCCGTATGTTGGTTTTTCTCCCTCTCCCTGTGGGCGGGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                 <  2:282856/150‑1 (MQ=25)
                                                       gacCCGGCAAACACTGCCCCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTGTTCCCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                              <  2:8897837/150‑6 (MQ=37)
                                                        acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCCCCCACCCCCCAACCCAATCCCCCCAGGGGGGGGGGGcccc                                                                                             >  2:3532602/1‑148 (MQ=25)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAACCCCACCCACAACGGGGCGGGGGccccgg                                                                                           >  2:3991899/1‑146 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAACCCAAAACCCCCAGGGAGGGGGGGCCGagag                                                                                          >  1:3060676/1‑147 (MQ=25)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCAACCCCAAACCACACCCAAAACGGCGCGGGGGCCGAACg                                                                                          >  2:1064875/1‑150 (MQ=25)
                                                            ggaaaaTCCTGCCCCGCCGCGCCGGACTAAAATGTTTTTGTTTGGTTGTTTTTCCCCCTCCCCGGGTGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                         <  1:6583960/147‑1 (MQ=255)
                                                            ggCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                         <  2:1001298/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGTACGAAAATTTATTCGTAAGGTTGGTTTGCACCCTCTCCCTGTGGGTTTGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                        <  2:2874083/150‑1 (MQ=255)
                                                               aactACGGCCCCGACGCGCGGGACGAACAAGTATTCGTAAGGTTGGTTTGCGCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                      <  2:6668717/147‑1 (MQ=38)
                                                               aaaTACGGCCCGGCCGCGCCGGCCGAAACAGTATGCGAACGGGTGGTTTCCGCCCCCCCCCTGGGGGCGTGGGCCGGGGTGGGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                      <  2:4076559/150‑1 (MQ=255)
                                                                    cggCACCGACGCGACCGACGCAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCGGTGGGAGCGGGCCGGGGTGAGGCCATCAGCGCGCACGGTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                 <  1:9048272/148‑1 (MQ=255)
                                                                       ccccgccgCTACGTACTAAAATTTATTCGTCAGTTTGGTTTCCCCCCTCTCCGTGTGGGTGGGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:5507844/148‑1 (MQ=255)
                                                                       cACCGACGCGACGGACGAAAAAGTTTTCTTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                              <  1:8428779/150‑1 (MQ=38)
                                                                         ccGAGGCTACGTACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGGGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  2:4210080/150‑1 (MQ=25)
                                                                         ccGAGGCGAGGTACTAAAAATTTTTCGTATGGTGGGTTTTCCCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  1:7305718/150‑1 (MQ=255)
                                                                         ccGACGCGACGGCCTAAAAAGTATTCGTATGGTTGTTTTTCCCCCCCCCCCTGTGGGTGTGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  1:2940736/150‑1 (MQ=255)
                                                                         ccGACGCCCCGGACTAAAATTTATTGTTATGGTTGTTTTTCGCCCGCTCCGTGGGGGTTTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                            <  1:8624760/150‑1 (MQ=255)
                                                                                 ccggacgAAAAAGTTTTCGTACGGTTGGTTTCCCCCCTCGCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                    <  2:5875190/149‑1 (MQ=39)
                                                                                     actaaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                                                                <  1:6785602/147‑1 (MQ=255)
                                                                                           caGTTTTCGTACGGTTGGTTTCCCCCCCCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  1:234183/149‑1 (MQ=255)
                                                                                           aattATTGGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  1:8108305/147‑1 (MQ=255)
                                                                                           aaGTATTCGTCCGGTTGGTTTTCCCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                          <  1:6640881/150‑1 (MQ=255)
                                                                                            aGTATTCGTAGGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                         <  2:2171132/150‑1 (MQ=255)
                                                                                             gTATTCGTACGGTTGTTTTTCTCCCTCACCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                        <  1:5279333/150‑1 (MQ=255)
                                                                                              ttttCTTATGGTTGTTTTTCCCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  2:3060676/148‑1 (MQ=38)
                                                                                              tATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  2:3828994/150‑1 (MQ=255)
                                                                                               atccGTAAGGTTGGTTTTCTCCCTGTCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                      <  1:5530059/147‑1 (MQ=255)
                                                                                                  cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGACATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                   <  2:2626557/150‑1 (MQ=255)
                                                                                                   gTATGGTTGGTTTTCACCCTCTCCCTGTGGTTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                  <  2:3767033/150‑1 (MQ=255)
                                                                                                     tgtgttggttTTCCCCCGCGCCTTGGGGTTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  1:3141778/147‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                               <  2:930394/150‑1 (MQ=255)
                                                                                                        tttggttTTCTCCCACTCCGTGTGGGAGAGGGCCGGGGTTAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                             <  1:5572983/149‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                          <  2:5152948/150‑1 (MQ=255)
                                                                                                            gtttagcccctctccctGTGGTAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                         <  2:2512003/143‑1 (MQ=255)
                                                                                                            gttTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                         <  1:5292768/150‑1 (MQ=255)
                                                                                                             tttgccccctctccgtGTGGGATAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                        <  2:4316980/144‑1 (MQ=255)
                                                                                                               ttCTCCCTCTCCCTGTGGTTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                      <  1:2074088/150‑1 (MQ=255)
                                                                                                                 ctccctcgccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGTGACCGATCGAGCACAAATTTTGCGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                    <  2:8268144/150‑1 (MQ=255)
                                                                                                                  tccctctccctGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                   <  2:7184183/150‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGTAAATCg                                  <  2:6749318/150‑1 (MQ=255)
                                                                                                                   ccccctccgtGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                  <  2:5196058/150‑1 (MQ=255)
                                                                                                                    cctccccgtGGGGGTGTGGCCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                 <  1:7179590/150‑1 (MQ=255)
                                                                                                                    cctccccctGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                 <  1:771795/150‑1 (MQ=255)
                                                                                                                          ccggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                           <  1:7638737/147‑1 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                         <  1:5400297/150‑1 (MQ=255)
                                                                                                                              tgGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                       <  1:220623/150‑1 (MQ=255)
                                                                                                                                ggAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                     <  1:1090918/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                   <  1:9009034/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCACACCCAAAACCCCTCCCACCAGGGGGAGGGGGCCGATCCAGCCAAAAATTTGTGGGATATAATAATATGTCTTGTGTCACAAGAACAGAAAAAAGGGCGTTTTCGCCGGa                   >  2:2312963/1‑150 (MQ=255)
                                                                                                                                   gtgggCGGGGGTGCGGGCATCAGCGAGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                  <  1:8528213/148‑1 (MQ=255)
                                                                                                                                        cgggggggCGGGCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  1:7420216/148‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  1:6047897/150‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg             <  2:7634021/150‑1 (MQ=255)
                                                                                                                                         gggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  2:8070788/149‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  2:7680615/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  2:7548747/150‑1 (MQ=255)
                                                                                                                                          gggtTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  2:1910030/150‑1 (MQ=255)
                                                                                                                                          ggggggCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  2:1853853/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  1:3953798/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  1:9170036/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg           <  2:8133480/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga         <  1:7664352/150‑1 (MQ=255)
                                                                                                                                                   gCATCAGCGCGCACGTTCAACCACACCCAAAACCAAACCCACAAGGGAGAGGGGACCCAAAGAGCACAAATTTTTTGAGATAAATTAATAATTCTTTTAAAAAAAGAACACAAAAAACGGCGTTTTCGCCGGAAAAATTCGGAGAAAAAg  >  2:5258601/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  NC_000913/164432‑164728

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: