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.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  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)
                     gAATGGGCTTTATCCAAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGAGGAAAAAGTTTTGTTACGGTTGGTTTTCTCCCTCCCCCTTTGGGTTTGGGCCGGGGTTTGGGCATCAGCGCGCACGTtcaccctc                                                                                                                               <  2:4661601/150‑1 (MQ=255)
                       atggggGTTATCCAAAACTTGTCGGGCCGGACGACCCGGAAAATACGGCCCCGACGCGCCGCACGAAACTTTATTCGTACGGTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                             <  2:1312967/148‑1 (MQ=255)
                             gTTATCCAAAACTTGTCTGGCGGCCCGCCCCGGAAAATCCTGCACCGAGGCGACGGAGTAAAATTTATTGGTCAGGTTGTTTTTCTCCGTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                       <  2:2447039/150‑1 (MQ=255)
                              ttATCCCAACCTTGTCTGGCCGGACGCCCCGGCAAATACTCCACCGACGCGACGTACTAAAAAGTTTGCGTATGGTGGTTTTTCCCCACGGCCCTGTGGTAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:1950094/150‑1 (MQ=255)
                                  ccaaaaCTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGCCGCGACGTACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                  <  1:5524488/147‑1 (MQ=25)
                                  ccaaaaCATGTCTGGCCGGACCACCGGTAAAATACGGCACCGACGCGACGGACTAAAAATTTTTCTTATGGTTGTTTTTCCCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                  <  1:6015293/147‑1 (MQ=255)
                                    cAAACTTGTGGGGTCGGACGACCCGGAAAATCCTGCACCGACGCTACGTACTAAAAAGTATTCGTACGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                <  2:6268095/150‑1 (MQ=25)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAAcccacacc                                                                                                                >  1:5368993/1‑145 (MQ=37)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCGTACCCTCTccc                                                                                                               >  1:5458842/1‑150 (MQ=38)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCCAAaccacaacc                                                                                                               >  1:817801/1‑144 (MQ=38)
                                        cttgtgtgGCCGGACTACCCAGAAAATACTGCAGCGCCGCGAGGGACGAAAATTTTTTCTTATTGTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:3030649/148‑1 (MQ=25)
                                               ggCCGGACGACCCGGAAAATACTCCCCCCCCGCGCCGGACTAAAATGTTTCCTTACGGTTGTTTTTCCCCCACCCCCTGGGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  1:1615824/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:4833101/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCAAACCCCAAACCACACCCCCAAGGGggag                                                                                                    >  2:2231233/1‑150 (MQ=37)
                                                 ccGGACGCCCCGGCACAGCCGGCACCGCCGCGACGGACGAAAACGGATGCGCCCGGTGGGTTTTCGCGCTCACCCGGGGGGAGAGGGCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                   <  2:5174586/150‑1 (MQ=255)
                                                  cGGACGCCCCGGCAAAAACTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGGTTGGTTTGCTCCCTCGCCCTGGGGGTGAGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                  <  2:6456070/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCCCACCCCCACCCAAACCCGCAAGGGGGAgggg                                                                                                 >  2:3998738/1‑150 (MQ=25)
                                                     acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                               <  2:4966372/150‑1 (MQ=255)
                                                      cgccccGGCAAACACCGCCCCGACGCGACGCACGAAACCGTATTCGCCTGGTTGGTTGTCGCCCCCTCCCTGGGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                                              <  1:1384132/147‑5 (MQ=255)
                                                      cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCAACCCAAAACACACCCACAAGGGCGCGGGGGcc                                                                                              >  1:968795/1‑150 (MQ=25)
                                                       cacCCGGAAAAAACTCCACCCACCCGACCAACCAAAAACTATACTAAAGGTTGCTTTACCCCCTCCCCCGGTGGCACAGGGCCGGGTTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                             <  2:723345/149‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCCAAACCACACCCACCAGGGGGAGGGGGccccgc                                                                                          >  2:2974881/1‑146 (MQ=25)
                                                           cgccAAATACGGCACCGCCGCGACGGCCGAAAAAGTATGCGTACGGTGGGTTTTCGCCCCCCCCCGGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                         <  2:3557927/147‑1 (MQ=255)
                                                           cGGCAAATACTGCACCGACACGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCCCCCAAAACCGCTCCCCCACCGGAGGGGGGccccctcg                                                                                         >  2:6015797/1‑145 (MQ=25)
                                                            ggCAAATACTCCCCCGACGCGACGGACGAAAAATTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                        <  1:5955182/150‑1 (MQ=38)
                                                             gaaaaTCCGGCACCCACGCGAGGCAACAAAATTTATCCGTATGGTTGTTTTCCTCCCCCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  1:879323/148‑1 (MQ=255)
                                                             gCAAATCCTCCCCCGCCGCGCCGGACGAACACGTATGCGTCCGGTGGGTTGTCGCCCTCGCCCGGGGGGAGCGGGCGGGGGTGAGGGCCTCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  1:5941704/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  2:3784338/150‑1 (MQ=255)
                                                               aaaTCCCGCACCGACGCGACGTCCTAAAACTTTTGCGTTGGTTTGTTTTTCCCCCCCCCCGGGGGGTTGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                     <  2:3003774/150‑1 (MQ=255)
                                                                 aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAACCCACACCCACAACGGAGAGGGGACCGAGCCCGcaca                                                                                   >  2:4695174/1‑150 (MQ=25)
                                                                 aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACCCCCAAAACCACAACCACCACGGGGGGGGGTCCGATCGGGcacc                                                                                   >  2:1467884/1‑149 (MQ=21)
                                                                   aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCACCCCCACCCAAACCCCCACCCACAACGGGGGGGGGCCCCAACGAGGCCaaa                                                                                 >  1:4803146/1‑150 (MQ=21)
                                                                    cTGCACCGACGCGACGGCCGAAAAAGTATTCGTATGTTTGTTTTTCTCCCTCTCCGTGTGGTAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                <  1:1528060/150‑1 (MQ=38)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                <  2:2180870/150‑1 (MQ=255)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCCAACCCACTACCGCACGGGGGAGGGGGCCGCGCGAGCCCAAAt                                                                                >  1:1535442/1‑150 (MQ=25)
                                                                    cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTTGGCCCTCACCCAAAACCACTACCACAACGGGGGGGGGGCCGCTCGTGCCCaaaa                                                                                >  1:2188483/1‑149 (MQ=25)
                                                                       cACCTAGGCGAGGGACTAAAATTTTTTCTTTTGTTTGTTTTCCCCCCTCTCGTTGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  1:1901568/150‑1 (MQ=255)
                                                                       cACCTACGCGACGGACGAAAAATTATTCGTAAGTTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  2:5416555/150‑1 (MQ=255)
                                                                       cACCGACGCGCCGGCCGAACAGTTATGCGTCAGGTTGGTTTGCTCCCTCGCCCGGGGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  2:3500823/150‑1 (MQ=255)
                                                                        aCCGACGCGACGGACTAAAAAGTATTCTTACGGTTGTTTTTCTCCCTCACCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                            <  1:892678/150‑1 (MQ=38)
                                                                        aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCGCACCCCAACCCGCACCCGCCCGGGGGCGGGGACCGCTCGAGCCCAACTTTTg                                                                            >  2:5486002/1‑150 (MQ=21)
                                                                         ccGACGCGCCGTACTAAAACGTATTCGTATGTTTGGTTTTCCCCCCCCCCCGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:1443152/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGACATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:1989526/150‑1 (MQ=255)
                                                                              gcgccgGACGAAAATGTATTGTTACGGTTGGTTTGCTCCCTCGCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:5825965/150‑1 (MQ=255)
                                                                                gacggactAAAAAGTTTTCGTACGGGTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                    <  2:5678797/150‑1 (MQ=38)
                                                                                 acggactAAAAATTATTCGTATGTTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                   <  1:1431105/150‑1 (MQ=38)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACAACCCACACGGAGTGGGGGCCCCTCGAGCCCAAAAATTTGGAGatataat                                                                 >  1:4181174/1‑147 (MQ=16)
                                                                                        aaaatGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                            <  1:5488491/150‑1 (MQ=255)
                                                                                           aaGTATTCGTAAGGTTGGGTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  2:3728265/150‑1 (MQ=255)
                                                                                           aaGTAGTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  2:4138329/150‑1 (MQ=255)
                                                                                               ttgcGTATGGTTGTTTTTCTCCCCCGCCCTGGGGGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:4580986/147‑1 (MQ=255)
                                                                                               aTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  1:3198334/150‑1 (MQ=255)
                                                                                               aTTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  2:6466486/150‑1 (MQ=255)
                                                                                                 tCGTAAGGTTGGTTTTCGCCCCCGCCGGGTGGGAGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  2:1394487/150‑1 (MQ=255)
                                                                                                            gttTTCACCCTCTCCCTGTGGGTGAGGGCCGGTGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                        <  1:3586154/150‑1 (MQ=255)
                                                                                                                    cgtctccgtGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:4811693/148‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGTGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  2:6360588/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:3895175/150‑1 (MQ=255)
                                                                                                                     ctctcccgGGGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCGAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGCGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:918889/150‑1 (MQ=255)
                                                                                                                       ctccctGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGTTACCCCACACCCCAAACCACACCCACAAGGGCGCGGGGGCCGCGCGAGCACAAAATTTGAGAGATATCATCATATCTCATGTAACAGAAGAACAGAAAAACGGGcg                             >  2:4021923/1‑149 (MQ=255)
                                                                                                                        tccctGTGGGAGAGGGCCGGGGTGGGGGCATCAGCGCGCACGGGCACCCGCACCCAAACCCCCACCCCCCCGGGGGAGGGGGCCGAGCGAGCCCAAAATTTGAGAGATATATTCTTCTGGCTTTTAACCGAAGAACCGAAAATCGGGcat                            >  2:4397565/1‑148 (MQ=255)
                                                                                                                         ccctGGGGGTGGGGGCGGGGGTGCGGGCCTCCGCGCGCCCGTTCACCCGCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:2070001/150‑1 (MQ=255)
                                                                                                                            tgtgGGAGCGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                        <  1:843613/150‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACTACCACAAGGGGGAGGGGCCCGATCGAGCCCAAAATTTTAGAGATATCATCTTCTGTCTTGTAAAAGAAGGACACAAAAACGGGCTTTTGCGccg                    >  2:3240468/1‑149 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCGCACCCAAAAACACTCCCTACAGGGGGAGGGGCCCGAAAGAGCACAAAAATTGAGAGAAATAATAATCAGTATTGTAAAAGAATAACAGAAAAAAGGGCCTTTTCGCCAGAatat              >  2:4243546/1‑150 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTGCACCCCACCCCAACCCACAACCACAACGGGGGGGGGGCCCCTCGAGCCCAAAAATTTTTAGTTATATTATTATGTCTTTTTAAAAAAAAAAAAAAAAAAGGGCTTTTTCGCCGGAaaat              >  2:1330541/1‑150 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACACCCACCACGGCGAGGGGCCCCATCGACCACAAAAATTGAGAGATATATTAATAAGTCTTGTAACACAAGAAAACAAAAACAGGCGTTTTCGCCCGAaaat              >  2:5908138/1‑150 (MQ=255)
                                                                                                                                        ccGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCGCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  2:5139178/150‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  1:2875238/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:2689992/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:291181/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:397472/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  1:403777/150‑1 (MQ=255)
                                                                                                                                             gTGAGGGCATCAGCGCGCACGTTCACCCTCAACCTAAACCACTCCCTCAAGGGAGAGGGGACCGATCGAGAACAAAATTTGAGATATATCTTATTCTGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTTCGCCGGAATATTTCGgag       >  2:1573620/1‑150 (MQ=255)
                                                                                                                                              tGAGGGCCTCAGCGCGCACGCTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga      <  2:5349801/150‑1 (MQ=255)
                                                                                                                                                aGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:5150054/150‑1 (MQ=255)
                                                                                                                                                aGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGATGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:5514012/150‑1 (MQ=255)
                                                                                                                                                aGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa    <  1:3181593/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  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: