Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,562 A→T 18.4% intergenic (+28/‑168) 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,5620AT18.4% 369.7 / 24.8 163intergenic (+28/‑168)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (56/70);  new base T (2/32);  total (59/104)
Fisher's exact test for biased strand distribution p-value = 1.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

ACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCT  >  NC_000913/164413‑164709
                                                                                                                                                     |                                                                                                                                                   
aCTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGga                                                                                                                                                     >  1:4741905/1‑150 (MQ=255)
 cTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgag                                                                                                                                                    >  2:4195077/1‑150 (MQ=255)
 cTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgag                                                                                                                                                    >  2:162384/1‑150 (MQ=255)
  ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                   >  2:3073138/1‑150 (MQ=255)
  ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                   <  1:1973376/150‑1 (MQ=255)
  gtCTGGAAAGGAGCGTACCGTGAGGTGCAAAACGCGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                   <  2:3143797/149‑1 (MQ=255)
    cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                                 <  2:3013296/150‑1 (MQ=255)
    cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                                 <  2:2559928/150‑1 (MQ=255)
    cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                                 <  1:1835497/150‑1 (MQ=255)
     tGGAAAGGAGCGTACCGGGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGGTTTCTCCCTCTCCCTGTGGGCGAggg                                                                                                                                                >  2:4811693/1‑150 (MQ=255)
      ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                               <  1:3872390/150‑1 (MQ=255)
      ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                               <  1:2211003/150‑1 (MQ=255)
       gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                              >  2:2551092/1‑150 (MQ=255)
       gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                              <  2:5476369/150‑1 (MQ=255)
        aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                             >  2:549632/1‑150 (MQ=255)
         aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                            <  2:5230565/150‑1 (MQ=255)
          aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                           <  2:680485/150‑1 (MQ=255)
            gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                         >  2:573783/1‑150 (MQ=255)
               cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                      <  2:497974/150‑1 (MQ=255)
                  aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                   >  2:2187300/1‑150 (MQ=255)
                   ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                  >  2:5980460/1‑150 (MQ=255)
                   ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                  >  1:2942905/1‑150 (MQ=255)
                    cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                 <  2:2070683/150‑1 (MQ=255)
                    cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                 <  2:3515431/150‑1 (MQ=255)
                    cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                 >  2:4098291/1‑150 (MQ=255)
                      tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                               >  2:3316024/1‑150 (MQ=255)
                         ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCcgcg                                                                                                                            >  1:3011309/1‑150 (MQ=255)
                         ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCcgcg                                                                                                                            >  1:6360588/1‑150 (MQ=255)
                          gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                           >  2:4666589/1‑150 (MQ=255)
                             cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                        >  2:737013/1‑150 (MQ=255)
                              aaaatGAGTTTAAAGGGCGTTATCCCAAACAGTTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                       <  2:1617199/150‑1 (MQ=37)
                               aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                      >  1:5678797/1‑150 (MQ=255)
                               aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                      >  2:5210867/1‑150 (MQ=255)
                                aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                     >  1:1259459/1‑150 (MQ=255)
                                  agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGATGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                   >  1:3471981/1‑148 (MQ=255)
                                  agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGtgc                                                                                                                   >  1:5825965/1‑148 (MQ=255)
                                  agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttg                                                                                                                   >  2:403777/1‑149 (MQ=255)
                                  agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                   >  1:2540762/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)
                                                gTTAACCCAAACATGTGTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                     <  2:6000734/150‑1 (MQ=255)
                                                 ttATCCCAACCTTGTCTGGCCGGACGCCCCGGCAAATACTCCACCGACGCGACGTACTAAAAAGTTTGCGTATGGTGGTTTTTCCCCACGGCCCTGTGGTAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                    <  2:1950094/150‑1 (MQ=255)
                                                      ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTTTTCGTACGGTGGGTTTTCTCCCTCTCCGTGTGGTTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                               <  2:1597188/150‑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)
                                                            aggTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCCGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                         <  2:1004568/148‑1 (MQ=255)
                                                              gTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCttgagg                                                                                       <  2:448838/150‑6 (MQ=38)
                                                                  ggCCGGACGACCCGGAAAATACTCCCCCCCCGCGCCGGACTAAAATGTTTCCTTACGGTTGTTTTTCCCCCACCCCCTGGGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                   <  1:1615824/150‑1 (MQ=255)
                                                                     cGGACGCCCCGGCAAAAACTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGGTTGGTTTGCTCCCTCGCCCTGGGGGTGAGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                <  2:6456070/150‑1 (MQ=255)
                                                                        acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                             <  2:4966372/150‑1 (MQ=255)
                                                                         cgccccGGCAAACACCGCCCCGACGCGACGCACGAAACCGTATTCGCCTGGTTGGTTGTCGCCCCCTCCCTGGGGGAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttag                                                                            <  1:1384132/147‑5 (MQ=255)
                                                                         cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCCAACCCACGCCCACAACGGGGAGGGGGcc                                                                            >  1:968344/1‑150 (MQ=25)
                                                                              cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCACTCCCCCGAGGGGGCGGGGGCCGATCg                                                                       >  2:1545810/1‑150 (MQ=37)
                                                                              cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCGCCCCCCAACCCACACCCCCCAGGGGGGGGGGGCCGCTCg                                                                       >  1:6068727/1‑150 (MQ=25)
                                                                               ggCAAATACTCCCCCGACGCGACGGACGAAAAATTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                      <  1:5955182/150‑1 (MQ=38)
                                                                                gaaaaTCCGGCACCCACGCGAGGCAACAAAATTTATCCGTATGGTTGTTTTCCTCCCCCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                     <  1:879323/148‑1 (MQ=255)
                                                                                  aaaTCCCGCACCGACGCGACGTCCTAAAACTTTTGCGTTGGTTTGTTTTTCCCCCCCCCCGGGGGGTTGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                   <  2:3003774/150‑1 (MQ=255)
                                                                                  aaaTACTGCACCGACGCGACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                   <  2:3824532/150‑1 (MQ=255)
                                                                                      aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCACCCCCACCCAAACCCCCACCCACAACGGGGGGGGGCCCCAACGAGGCCaaa                                                               >  1:4803146/1‑150 (MQ=21)
                                                                                       cTGCCCCGACGCGCCGGACTAAAAAGTATTCGTATGTTTGTTTTTCCCCATCTCCCGGGGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                              <  2:2415146/150‑1 (MQ=255)
                                                                                       cTGCACCGACGCGACGGCCGAAAAAGTATTCGTATGTTTGTTTTTCTCCCTCTCCGTGTGGTAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                              <  1:1528060/150‑1 (MQ=38)
                                                                                       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)
                                                                                          cACCGAGGCGAGGGACTAAAATTTTTTGTTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                           <  2:2075328/150‑1 (MQ=255)
                                                                                          cACCGACGCGCCGGCCGAACAGTTATGCGTCAGGTTGGTTTGCTCCCTCGCCCGGGGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                           <  2:3500823/150‑1 (MQ=255)
                                                                                            ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGACATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                         <  1:1512996/150‑1 (MQ=255)
                                                                                            ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGACATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                         <  1:1989526/150‑1 (MQ=255)
                                                                                                   gactgactAAAATTTTTTCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                  <  1:4032308/150‑1 (MQ=255)
                                                                                                   gacggactAAAAAGTTTTCGTACGGGTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                  <  2:5678797/150‑1 (MQ=38)
                                                                                                    acggactAAAAATTATTCGTATGTTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                 <  1:1431105/150‑1 (MQ=38)
                                                                                                           aaaatGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                          <  1:5488491/150‑1 (MQ=255)
                                                                                                              aaGTATTCGTAAGGTTGGGTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                       <  2:3728265/150‑1 (MQ=255)
                                                                                                                  ttgcGTATGGTTGTTTTTCTCCCCCGCCCTGGGGGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                   <  1:4580986/147‑1 (MQ=255)
                                                                                                                  cgtCGTAAGGTGGGTTTGCTCCCTCGCCCTGGGGGAGAGGGCCGGGGGGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                   <  2:5972941/148‑1 (MQ=255)
                                                                                                                  aTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                   <  1:3198334/150‑1 (MQ=255)
                                                                                                                  aTTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                   <  2:6466486/150‑1 (MQ=255)
                                                                                                                  aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                   <  1:2698516/150‑1 (MQ=255)
                                                                                                                    tggTATGGTTGGTTTGCGCCCTCTCCCTGTGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                 <  2:3668543/148‑1 (MQ=255)
                                                                                                                    tCGTAAGGTTGGTTTTCGCCCCCGCCGGGTGGGAGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                 <  2:1394487/150‑1 (MQ=255)
                                                                                                                         tggttggttTTCTCCCGCTCCCGGTGGGTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                            <  2:5864427/149‑1 (MQ=255)
                                                                                                                               gttTTCACCCTCTCCCTGTGGGTGAGGGCCGGTGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                      <  1:3586154/150‑1 (MQ=255)
                                                                                                                                      ccctccccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg               <  2:4238233/150‑1 (MQ=255)
                                                                                                                                       cctctccctGTGGGTGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg              <  2:6360588/150‑1 (MQ=255)
                                                                                                                                        ctctccctGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg             <  1:5980460/150‑1 (MQ=255)
                                                                                                                                        ctctcccgGTGGGCGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTGGAGAGATATCGTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg             <  2:2770392/150‑1 (MQ=255)
                                                                                                                                        ctctcccgGGGGGTGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCGAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGCGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg             <  2:918889/150‑1 (MQ=255)
                                                                                                                                          ctccctGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGTTACCCCACACCCCAAACCACACCCACAAGGGCGCGGGGGCCGCGCGAGCACAAAATTTGAGAGATATCATCATATCTCATGTAACAGAAGAACAGAAAAACGGGcg           >  2:4021923/1‑149 (MQ=255)
                                                                                                                                            ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt         <  2:4913413/150‑1 (MQ=255)
                                                                                                                                            cccgGTGGGAGTGGGCCGGGTTGTGGTCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt         <  1:2070305/150‑1 (MQ=255)
                                                                                                                                             cctGTGGGTGAGGGCCGGGGTGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt        <  1:5273343/150‑1 (MQ=255)
                                                                                                                                              cggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg       <  1:3792263/148‑1 (MQ=255)
                                                                                                                                               tgtgGGAGCGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc      <  1:843613/150‑1 (MQ=255)
                                                                                                                                               tgtgGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGGGCACCCACCCCCAAACCCACGCCCACAACGGAGGGGGGGCCGCTCGCCCCCCCATTTTTTGAGAGATATTATTGAGTCGTGTTACACAAGAACAGAAAAACggggattttg      >  2:5369979/1‑143 (MQ=255)
                                                                                                                                                   ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACTACCACAAGGGGGAGGGGCCCGATCGAGCCCAAAATTTTAGAGATATCATCTTCTGTCTTGTAAAAGAAGGACACAAAAACGGGCTTTTGCGccg  >  2:3240468/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCT  >  NC_000913/164413‑164709

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: