Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.5% 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.5% 256.7 / 16.8 156intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (48/77);  new base C (5/26);  total (53/103)
Fisher's exact test for biased strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164432‑164727
                                                                                                                                                     |                                                                                                                                                   
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                     >  2:5980460/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                    >  2:4098291/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                    <  2:2070683/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCag                                                                                                                                                    >  2:3196010/1‑149 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  <  2:5979978/150‑1 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  <  2:341394/150‑1 (MQ=255)
      ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCcgcg                                                                                                                                               >  1:3011309/1‑150 (MQ=255)
      ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCcgcg                                                                                                                                               >  1:6360588/1‑150 (MQ=255)
      ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                               >  1:6360501/1‑150 (MQ=255)
         gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                            >  2:2698516/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                           >  2:737013/1‑150 (MQ=255)
           aaaatGAGTTTAAAGGGCGTTATCCCAAACAGTTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                          <  2:1617199/150‑1 (MQ=37)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                          <  2:1617303/150‑1 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                         >  1:5678797/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                         >  1:1882543/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                         >  2:5210867/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                         >  1:5324647/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                        >  1:1259459/1‑150 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttg                                                                                                                                      >  2:403777/1‑149 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                      >  1:2540762/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                      >  1:3472026/1‑148 (MQ=255)
                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACggggcc                                                                                                                                    >  2:2052997/1‑145 (MQ=255)
                     gAATGGGCTTTATCCAAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGAGGAAAAAGTTTTGTTACGGTTGGTTTTCTCCCTCCCCCTTTGGGTTTGGGCCGGGGTTTGGGCATCAGCGCGCACGTt‑caccctc                                                                                                                               <  2:4661601/150‑1 (MQ=255)
                       atggggGTTATCCAAAACTTGTCGGGCCGGACGACCCGGAAAATACGGCCCCGACGCGCCGCACGAAACTTTATTCGTACGGTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTt‑caccctcac                                                                                                                             <  2:1312967/148‑1 (MQ=255)
                             gTTATCCAAAACTTGTCTGGCGGCCCGCCCCGGAAAATCCTGCACCGAGGCGACGGAGTAAAATTTATTGGTCAGGTTGTTTTTCTCCGTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAAc                                                                                                                       <  2:2447039/150‑1 (MQ=255)
                              ttATCCCAACCTTGTCTGGCCGGACGCCCCGGCAAATACTCCACCGACGCGACGTACTAAAAAGTTTGCGTATGGTGGTTTTTCCCCACGGCCCTGTGGTAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAAcc                                                                                                                      <  2:1950094/150‑1 (MQ=255)
                                  ccaaaaCTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGCCGCGACGTACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCctct                                                                                                                  <  1:5524488/147‑1 (MQ=25)
                                  ccaaaaCATGTCTGGCCGGACCACCGGTAAAATACGGCACCGACGCGACGGACTAAAAATTTTTCTTATGGTTGTTTTTCCCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCctct                                                                                                                  <  1:6015293/147‑1 (MQ=255)
                                    cAAACTTGTGGGGTCGGACGACCCGGAAAATCCTGCACCGACGCTACGTACTAAAAAGTATTCGTACGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTcc                                                                                                                <  2:6268095/150‑1 (MQ=25)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑GGCCCCCCCCCCAAcccacacc                                                                                                                >  1:5368993/1‑145 (MQ=37)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑GGCCCGCACCCGTACCCTCTccc                                                                                                               >  1:5458842/1‑150 (MQ=38)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑GACCCACACCCCAAaccacaacc                                                                                                               >  1:817801/1‑144 (MQ=38)
                                        cttgtgtgGCCGGACTACCCAGAAAATACTGCAGCGCCGCGAGGGACGAAAATTTTTTCTTATTGTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:3030649/148‑1 (MQ=25)
                                               ggCCGGACGACCCGGAAAATACTCCCCCCCCGCGCCGGACTAAAATGTTTCCTTACGGTTGTTTTTCCCCCACCCCCTGGGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  1:1615824/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:4833101/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑CACCCAAACCCCAAACCACACCCCCAAGGGggag                                                                                                    >  2:2231233/1‑150 (MQ=37)
                                                 ccGGACGCCCCGGCACAGCCGGCACCGCCGCGACGGACGAAAACGGATGCGCCCGGTGGGTTTTCGCGCTCACCCGGGGGGAGAGGGCGGGGGTGCGGGCCTCCGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                   <  2:5174586/150‑1 (MQ=255)
                                                  cGGACGCCCCGGCAAAAACTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGGTTGGTTTGCTCCCTCGCCCTGGGGGTGAGGGCCGGGGTGAGGCCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                  <  2:6456070/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTG‑GGCCCCCACCCCCACCCAAACCCGCAAGGGGGAgggg                                                                                                 >  2:3998738/1‑150 (MQ=25)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                               <  2:4966372/150‑1 (MQ=255)
                                                      cgccccGGCAAACACCGCCCCGACGCGACGCACGAAACCGTATTCGCCTGGTTGGTTGTCGCCCCCTCCCTGGGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                                              <  1:1384132/147‑5 (MQ=255)
                                                      cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑GGCCCCCAACCCAAAACACACCCACAAGGGCGCGGGGGcc                                                                                              >  1:968795/1‑150 (MQ=25)
                                                       cacCCGGAAAAAACTCCACCCACCCGACCAACCAAAAACTATACTAAAGGTTGCTTTACCCCCTCCCCCGGTGGCACAGGGCCGGGTTGCGGCCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                             <  2:723345/149‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑GACCCGCACCCCAAACCACACCCACCAGGGGGAGGGGGccccgc                                                                                          >  2:2974881/1‑146 (MQ=25)
                                                           cgccAAATACGGCACCGCCGCGACGGCCGAAAAAGTATGCGTACGGTGGGTTTTCGCCCCCCCCCGGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                         <  2:3557927/147‑1 (MQ=255)
                                                           cGGCAAATACTGCACCGACACGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑CACCCCCCCCCAAAACCGCTCCCCCACCGGAGGGGGGccccctcg                                                                                         >  2:6015797/1‑145 (MQ=25)
                                                            ggCAAATACTCCCCCGACGCGACGGACGAAAAATTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                        <  1:5955182/150‑1 (MQ=38)
                                                             gaaaaTCCGGCACCCACGCGAGGCAACAAAATTTATCCGTATGGTTGTTTTCCTCCCCCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  1:879323/148‑1 (MQ=255)
                                                             gCAAATCCTCCCCCGCCGCGCCGGACGAACACGTATGCGTCCGGTGGGTTGTCGCCCTCGCCCGGGGGGAGCGGGCGGGGGTGAGGGCCTCCGCGCGCCCGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  1:5941704/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  2:3784338/150‑1 (MQ=255)
                                                               aaaTCCCGCACCGACGCGACGTCCTAAAACTTTTGCGTTGGTTTGTTTTTCCCCCCCCCCGGGGGGTTGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                     <  2:3003774/150‑1 (MQ=255)
                                                                 aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑GACCCCCACCCCAACCCACACCCACAACGGAGAGGGGACCGAGCCCGcaca                                                                                   >  2:4695174/1‑150 (MQ=25)
                                                                 aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑CCCCCACCCCCAAAACCACAACCACCACGGGGGGGGGTCCGATCGGGcacc                                                                                   >  2:1467884/1‑149 (MQ=21)
                                                                   aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGG‑CACCCCCACCCAAACCCCCACCCACAACGGGGGGGGGCCCCAACGAGGCCaaa                                                                                 >  1:4803146/1‑150 (MQ=21)
                                                                    cTGCACCGACGCGACGGCCGAAAAAGTATTCGTATGTTTGTTTTTCTCCCTCTCCGTGTGGTAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                <  1:1528060/150‑1 (MQ=38)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                <  2:2180870/150‑1 (MQ=255)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑GCCCCGCACCCCAACCCACTACCGCACGGGGGAGGGGGCCGCGCGAGCCCAAAt                                                                                >  1:1535442/1‑150 (MQ=25)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTT‑GGCCCTCACCCAAAACCACTACCACAACGGGGGGGGGGCCGCTCGTGCCCaaaa                                                                                >  1:2188483/1‑149 (MQ=25)
                                                                       cACCTAGGCGAGGGACTAAAATTTTTTCTTTTGTTTGTTTTCCCCCCTCTCGTTGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  1:1901568/150‑1 (MQ=255)
                                                                       cACCTACGCGACGGACGAAAAATTATTCGTAAGTTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  2:5416555/150‑1 (MQ=255)
                                                                       cACCGACGCGCCGGCCGAACAGTTATGCGTCAGGTTGGTTTGCTCCCTCGCCCGGGGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGGT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  2:3500823/150‑1 (MQ=255)
                                                                        aCCGACGCGACGGACTAAAAAGTATTCTTACGGTTGTTTTTCTCCCTCACCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                            <  1:892678/150‑1 (MQ=38)
                                                                        aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGG‑CCCCCGCACCCCAACCCGCACCCGCCCGGGGGCGGGGACCGCTCGAGCCCAACTTTTg                                                                            >  2:5486002/1‑150 (MQ=21)
                                                                         ccGACGCGCCGTACTAAAACGTATTCGTATGTTTGGTTTTCCCCCCCCCCCGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:1443152/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGACATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:1989526/150‑1 (MQ=255)
                                                                              gcgccgGACGAAAATGTATTGTTACGGTTGGTTTGCTCCCTCGCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:5825965/150‑1 (MQ=255)
                                                                                gacggactAAAAAGTTTTCGTACGGGTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                    <  2:5678797/150‑1 (MQ=38)
                                                                                 acggactAAAAATTATTCGTATGTTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                   <  1:1431105/150‑1 (MQ=38)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑GACCCACACCCAAAACCACAACCCACACGGAGTGGGGGCCCCTCGAGCCCAAAAATTTGGAGatataat                                                                 >  1:4181174/1‑147 (MQ=16)
                                                                                        aaaatGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                            <  1:5488491/150‑1 (MQ=255)
                                                                                           aaGTATTCGTAAGGTTGGGTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  2:3728265/150‑1 (MQ=255)
                                                                                           aaGTAGTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  2:4138329/150‑1 (MQ=255)
                                                                                               ttgcGTATGGTTGTTTTTCTCCCCCGCCCTGGGGGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:4580986/147‑1 (MQ=255)
                                                                                               aTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:3198334/150‑1 (MQ=255)
                                                                                               aTTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  2:6466486/150‑1 (MQ=255)
                                                                                                 tCGTAAGGTTGGTTTTCGCCCCCGCCGGGTGGGAGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  2:1394487/150‑1 (MQ=255)
                                                                                                            gttTTCACCCTCTCCCTGTGGGTGAGGGCCGGTGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                        <  1:3586154/150‑1 (MQ=255)
                                                                                                                    cgtctccgtGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:4811693/148‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGTGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  2:6360588/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:3895175/150‑1 (MQ=255)
                                                                                                                     ctctcccgGGGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTT‑CACCCTCACCCGAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGCGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:918889/150‑1 (MQ=255)
                                                                                                                       ctccctGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGTT‑ACCCCACACCCCAAACCACACCCACAAGGGCGCGGGGGCCGCGCGAGCACAAAATTTGAGAGATATCATCATATCTCATGTAACAGAAGAACAGAAAAACGGGcg                             >  2:4021923/1‑149 (MQ=255)
                                                                                                                        tccctGTGGGAGAGGGCCGGGGTGGGGGCATCAGCGCGCACGGG‑CACCCGCACCCAAACCCCCACCCCCCCGGGGGAGGGGGCCGAGCGAGCCCAAAATTTGAGAGATATATTCTTCTGGCTTTTAACCGAAGAACCGAAAATCGGGcat                            >  2:4397565/1‑148 (MQ=255)
                                                                                                                         ccctGGGGGTGGGGGCGGGGGTGCGGGCCTCCGCGCGCCCGTT‑CACCCGCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:2070001/150‑1 (MQ=255)
                                                                                                                            tgtgGGAGCGGGCGGGGGTGAGGCCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                        <  1:843613/150‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTG‑CACCCACACCCAAAACCACTACCACAAGGGGGAGGGGCCCGATCGAGCCCAAAATTTTAGAGATATCATCTTCTGTCTTGTAAAAGAAGGACACAAAAACGGGCTTTTGCGccg                    >  2:3240468/1‑149 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCCCGTT‑CACCCGCACCCAAAAACACTCCCTACAGGGGGAGGGGCCCGAAAGAGCACAAAAATTGAGAGAAATAATAATCAGTATTGTAAAAGAATAACAGAAAAAAGGGCCTTTTCGCCAGAatat              >  2:4243546/1‑150 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTGCACCC‑CACCCCAACCCACAACCACAACGGGGGGGGGGCCCCTCGAGCCCAAAAATTTTTAGTTATATTATTATGTCTTTTTAAAAAAAAAAAAAAAAAAGGGCTTTTTCGCCGGAaaat              >  2:1330541/1‑150 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCGCACCCAAAACCACACCCACCACGGCGAGGGGCCCCATCGACCACAAAAATTGAGAGATATATTAATAAGTCTTGTAACACAAGAAAACAAAAACAGGCGTTTTCGCCCGAaaat              >  2:5908138/1‑150 (MQ=255)
                                                                                                                                        ccGGGGTGCGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCGCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  2:5139178/150‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGCCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  1:2875238/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:2689992/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGCCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:291181/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:397472/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  1:403777/150‑1 (MQ=255)
                                                                                                                                             gTGAGGGCATCAGCGCGCACGTT‑CACCCTCAACCTAAACCACTCCCTCAAGGGAGAGGGGACCGATCGAGAACAAAATTTGAGATATATCTTATTCTGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTTCGCCGGAATATTTCGgag       >  2:1573620/1‑150 (MQ=255)
                                                                                                                                              tGAGGGCCTCAGCGCGCACGCT‑CACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga      <  2:5349801/150‑1 (MQ=255)
                                                                                                                                                aGGGCCTCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:5150054/150‑1 (MQ=255)
                                                                                                                                                aGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGATGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:5514012/150‑1 (MQ=255)
                                                                                                                                                aGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:3181593/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTT‑CACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164432‑164727

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: ‑