Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,290,380 T→A 7.5% intergenic (‑91/+93) PP_4660 ← / ← mdaB LysR family transcriptional regulator/drug activity modulator B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,290,3800TA 7.5% 452.8 / 20.8 148intergenic (‑91/+93)PP_4660/mdaBLysR family transcriptional regulator/drug activity modulator B
Reads supporting (aligned to +/- strand):  ref base T (64/71);  new base A (3/8);  total (68/80)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 149 total aligned reads displayed.

ATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5290231‑5290528
                                                                                                                                                     |                                                                                                                                                      
aTCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATACCCGCTTGTCGAGGCCAAGGCCAGCGGCGCCTGTATGGGCCTCGTCGCGGGCTT                                                                                                                                                        <  2:450618‑M1/150‑1 (MQ=255)
      gaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGc                                                                                                                                                  >  1:1166036‑M1/128‑150 (MQ=255)
      gaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCGtgccccc                                                                                                                                                  <  1:2127008‑M1/23‑3 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                                <  2:2858362‑M1/150‑1 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                                <  1:2120360‑M1/150‑1 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGAGCAAGGATTCGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCga                                                                                                                                                <  2:4372530‑M1/150‑3 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGAGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCga                                                                                                                                                <  2:3583604‑M1/150‑3 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGAGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCga                                                                                                                                                <  2:4372706‑M1/150‑3 (MQ=255)
        ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                >  2:3903681‑M1/126‑150 (MQ=255)
        ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                <  2:784166‑M1/25‑1 (MQ=255)
        ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                >  2:1846592‑M1/126‑150 (MQ=255)
        ccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCCcgcga                                                                                                                                                >  2:853724‑M1/126‑148 (MQ=255)
         cggcccgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                               >  2:3303751‑M1/125‑147 (MQ=255)
         cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaa                                                                                                                                               <  2:1819449‑M1/150‑4 (MQ=255)
         cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaa                                                                                                                                               <  1:1333514‑M1/150‑4 (MQ=255)
           gcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                             >  2:423508‑M1/123‑150 (MQ=255)
           gcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                             >  2:2679201‑M1/123‑150 (MQ=255)
           gTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGCgaaga                                                                                                                                             <  2:1663580‑M1/150‑6 (MQ=255)
             ctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGGGTGGCCgcgcccccc                                                                                                                                           <  1:982049‑M1/30‑4 (MQ=255)
              gATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                          <  2:1387397/150‑1 (MQ=255)
               aTGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAg                                                                                                                                         <  2:2662911/150‑1 (MQ=255)
                 gACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTa                                                                                                                                       <  2:475225/150‑1 (MQ=255)
                  aCGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTAc                                                                                                                                      <  1:543537/150‑1 (MQ=255)
                   cGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACt                                                                                                                                     <  2:1999053/150‑1 (MQ=255)
                      cgacgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCa                                                                                                                                  <  1:434580/150‑1 (MQ=255)
                       gacgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCaa                                                                                                                                 <  2:2334820/150‑1 (MQ=255)
                                cTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCAACCACAACTACAGCAAAATCCATGa                                                                                                                        >  1:3366441/1‑150 (MQ=255)
                                cTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGa                                                                                                                        >  2:3815593/1‑150 (MQ=255)
                                cTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCATGCCCGCTCCCACAAGTACTGCAAAGTCCATGa                                                                                                                        >  2:290643/1‑150 (MQ=255)
                                  ggAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGaaa                                                                                                                      >  1:1793161/1‑150 (MQ=255)
                                     gTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCACCCACAAGTAATGCAAAGTCCTTGAAAccg                                                                                                                   >  1:1098658/1‑149 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                                 <  1:729141/150‑1 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                                 <  1:2165682/150‑1 (MQ=255)
                                         ttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGt                                                                                                               >  2:1784903/1‑150 (MQ=255)
                                         ttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACATGCGt                                                                                                               >  2:741092/1‑150 (MQ=255)
                                           cGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTaa                                                                                                             <  1:2149304/150‑1 (MQ=255)
                                                cGGGTTTTCACATCAGTGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCt                                                                                                        >  1:2913963‑M2/1‑149 (MQ=255)
                                                cGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCt                                                                                                        >  1:2466490‑M2/1‑149 (MQ=255)
                                                 gggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGGCCAGCGGCGCCTGTCTGGGCCTCTTCGCGGGCGTGCCCGCTCCCACAACTACATCCAAGTCCTGGTAACCGGCGTCATAccag                                                                                                       >  2:1235900/1‑148 (MQ=255)
                                                  ggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                      >  1:1082106‑M2/1‑147 (MQ=255)
                                                  ggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                      >  2:2687945‑M2/1‑147 (MQ=255)
                                                   gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  2:437899‑M2/1‑146 (MQ=255)
                                                   gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  2:1808852‑M2/1‑146 (MQ=255)
                                                   gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGAGAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  1:3988875‑M2/1‑146 (MQ=255)
                                                    ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                    >  2:2420696‑M2/1‑145 (MQ=255)
                                                    ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                    >  1:851191‑M2/1‑145 (MQ=255)
                                                    ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCagtgg                                                                                                    >  1:3591882‑M2/1‑145 (MQ=255)
                                                    ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCCCGCCCCCAACTACAGCAAAGTACGTGAAACCGGCGGAATAACagtgg                                                                                                    <  1:3989597‑M2/150‑6 (MQ=255)
                                                        cacaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGGTTCCCGCGCCCCCAAGGACGGCCAAGTCCGTGAAAACTGCGGAATACCggtgggagc                                                                                                <  1:294249‑M2/150‑10 (MQ=255)
                                                         acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcg                                                                                               <  2:2179076‑M1/150‑45 (MQ=255)
                                                         acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAAACTGCGTAATACCggtgggggcg                                                                                               <  2:2091487‑M2/150‑11 (MQ=255)
                                                          ccggcggcgttggtagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            >  1:4398308‑M1/76‑150 (MQ=255)
                                                          ccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            >  2:1683001‑M1/76‑150 (MQ=255)
                                                          ccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            >  1:830791‑M1/76‑150 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgg                                                                                              <  1:1779042‑M1/150‑46 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgg                                                                                              >  1:585564‑M2/1‑139 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgg                                                                                              >  1:1424757‑M2/1‑139 (MQ=255)
                                                           cggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCggg                                                                                           >  2:3596921‑M1/75‑150 (MQ=255)
                                                           aTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGc                                                                                                                                   >  1:80519/1‑112 (MQ=255)
                                                            tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGTTGCCCGCTACCACAAGTACTGCAAAGTCCTTTTAAACTGGGTAATACCtgtgggtgcgggc                                                                                            <  2:3480756‑M2/150‑14 (MQ=255)
                                                            tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCTCCCACAACTACTGCAAAGTCCTTGAAACCTGCGGAATACCtgtgggagcgggc                                                                                            <  2:3879770‑M2/150‑14 (MQ=255)
                                                            tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAAAAACtgtgggagcgggc                                                                                            <  2:2127008‑M2/150‑14 (MQ=255)
                                                            tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACAGCACAGTCCTTGAAACCGGCGGAATACCagtgggagcgggc                                                                                            <  2:4194116‑M2/150‑14 (MQ=255)
                                                            tCAGGGCCTGTGGTTGGTCTGTGCAAGGATGAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGTGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtactcgtgggagcgggc                                                                                            <  2:4085034‑M1/150‑48 (MQ=255)
                                                             cAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGGACTGCAAAGTCCATGAAACCTGCGTAATACCtgtgggagcgggcc                                                                                           <  2:253886‑M2/150‑15 (MQ=255)
                                                                 gCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgggcaagcc                                                                                       >  1:2338747‑M2/1‑132 (MQ=255)
                                                                            gtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGGTGCCCCCGCCCCCCCAAACCGGCAAGTGCTTGAAACCGGCGGAATAACGGTGGGGGCGGGGCAGCCCCCGACGGGGc                                                                          >  1:2816851/1‑150 (MQ=255)
                                                                            gtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCCCGCCCACAACAAAAGCCAAAAACGTGAAAAAAGCCGAATAACAGGGGGTGCGGGCAAGCCCCCGAACAGGc                                                                          >  1:4224629/1‑150 (MQ=255)
                                                                            gtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTAAGGGCCTCTTCGCGGGCGTGCCCCCTCCCACCAGAAAGGCCAAATGCCTGAAAACCCCGGAATAACAGTGGGGGCCGGCCAAGCCCCGAATAGGc                                                                          >  1:1404902/1‑150 (MQ=255)
                                                                                  gCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGTTGCCCGCTCGCACAATTAGTGGAAAGTGGTTTTAAACTGGGTAATTACTGGGGGAGGGGGGAAGGCCGGGAAGAGGGGGGTag                                                                    >  2:2100908/1‑149 (MQ=255)
                                                                                   cAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGGTGCCCGCGCCCCCCAGGACGGCAAAGTCCATGGAAACCGCGGAATAACAGTGGGAGCGGGCCACCCCCCGGAGACGCCGGGAcc                                                                   >  1:1346400/1‑149 (MQ=255)
                                                                                   cAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCAAAAATACTGCAAAGTCCTTGAAACCTGCGGAATACCTGTGGGAGCGGGCAAACCCGCGAAGAGGCCGGTACa                                                                   >  2:440587/1‑150 (MQ=255)
                                                                                   cAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCCCCAGGACGGCCAACTCCGTGGAACCGGCGTAAAAACTGTGTGGGCCGGCCAGCCCCCCAAGAGGCCggggca                                                                   >  1:2741862/1‑146 (MQ=255)
                                                                                     ctgATTAGATCACCGCTTGTCTCGGCCACGGCCAGCGGCGCCTGTTTGGGCCTCTCGGGGGGCTTCCCGGCTCCCAAACTTACTGAAAAGTCCTTGAAACCTGCTAATTCCGGTTGGTTGCGGGAAAGCCCGCGAAGAGGCCGGTACAGa                                                                 <  2:1662140/148‑1 (MQ=255)
                                                                                       gtttAGATAACCGCTTGTCGAGGCCAAGTCCTGGGGCGCCTGTAGGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAATTACTGCAACGTCCTTTAAACCTGCTTAATACCTGTGGGAGCGGGCAAGTCCGCGAAGAGGCCGGTACAGACa                                                               <  2:605570/148‑1 (MQ=38)
                                                                                               aaCCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCTCCCACAAGTAAAGCAAAATGCTTGAAAACAGCGTAATAACTGTGGGAGCGGGCAAGCCCGCGCAGAGGCCCGTAAAGACAAAAAATTa                                                       >  2:1243822/1‑150 (MQ=37)
                                                                                                 ccGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGGTTGCCCGCTCCCACATGTACTGCAACGTCCTTTAAACCTGCTTAATACCTGTTGGAGCGGGAAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACt                                                     <  2:4019465/150‑1 (MQ=38)
                                                                                                     ttGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGACCACAAGTAATGGAAAGACCTTGAAAACATCGTAATTCCTGTGGGTGCGGGCCAGCCCGCGAAGAGGCCGGTAAAGACAAAAAAGTACAGTGc                                                 >  2:1195526/1‑150 (MQ=25)
                                                                                                      tGTCGAGGCCACGACAAGCGGGGCCTTTTTGGGCCTCTGCGCGGGCTTGCCCGCTCCCACATGTACTGAAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGGGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGcc                                                <  2:3070531/150‑1 (MQ=37)
                                                                                                        tCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGTCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACCTAAAATTAATGTGCCtt                                              >  2:612384/1‑150 (MQ=255)
                                                                                                              ccAACTCCAGCGGCGCCTGAATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGAACTCCAAACTCCTTGAAACCTGCGTAAAACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCaaa                                        <  2:4124396/150‑1 (MQ=37)
                                                                                                                aaGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACATCAAAGTGCTTGTAAACCGCGTAAAACCTGTGGGCGCGGGGCAGCCCGCGAAGCGGCCGGTGCAGACATAAAAATACTGTGGCTTGCCCAaca                                      >  2:3887090/1‑150 (MQ=25)
                                                                                                                     cAGCGGCGCCTGTATGGGCATCTTCGCGGGCTTGCCCGCTCCCACCAGTACTGCAACGTCCTTGAAACCTGCGTAATACCTGTGGTAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTg                                 <  2:2174019/150‑1 (MQ=255)
                                                                                                                      aGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGGAAAACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAAATACTGTGCCGTGCCAAACACCTTGt                                >  2:1101353/1‑150 (MQ=37)
                                                                                                                       gcggcgCCTGTATGGGCCTCTTGGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTc                               <  2:3179174/150‑1 (MQ=255)
                                                                                                                       gcggcgCCTGTATGGGCCTCTTGGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTc                               <  2:3179066/150‑1 (MQ=255)
                                                                                                                              cTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGTAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATgg                        <  2:3698332/150‑1 (MQ=255)
                                                                                                                               tGTATGGGCCTCTTCGCGGGCTTGCCCCCTCCCACAAGTACTGCAAACTACTTGAAACCAGCGGAAAACCGGTGGGAGCGGGCAAGCCCCCGAAGAGGCCGGTACAGACATAAAAATAATGTGCCATTCCAAACACCTTGTCCAGATtgg                       >  2:261467/1‑150 (MQ=25)
                                                                                                                                    gggCCTCTTCGCGGGCGTGCCCCCACCCACAACAACAGCAAACTCCGAGAAAACGGCCTAAAACCAGTGGGAGCGGGCAACCCCCCGAAGAGGCCGGTACAGACATAAAAATACAGAGCCATGCCAAACACCATGTCCCGATGGGCCTCa                  >  2:243135/1‑150 (MQ=255)
                                                                                                                                       CCTCTTCGCGGGCTGGCCCGCTCCCACACGTACTGCAACGTCCTTGAAACCGGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATaa               >  2:3366441‑M1/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCTCCCACAggtattacgccgggttcaaggactttgcagtacttgtgggagcgggcaagcccgcggagaggcccatacaggcgccgctggacttggccgcgaccagcgggtctctaatccttgcacagaccaa               >  2:4248496‑M1/1‑26 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgggcaagcccgcgaagaggcccatacaggcgccgctggacttggcctcgacaagcggttatctaatccgtgcacagaccaa               <  2:4444493‑M2/150‑91 (MQ=255)
                                                                                                                                          CTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTcataatct            <  2:3057280‑M1/150‑4 (MQ=255)
                                                                                                                                           TTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGCGCCTTGCCAAACACCTTGTCCAGATGGGCCGCATAAGCTg           <  2:2509481‑M1/150‑1 (MQ=255)
                                                                                                                                            TCGCGGGCTTGCCCGCTCCAACAAGTACTCCAATGTCCTTGAACCCGGCTTAATACCTGTGGTAGCGGGCACGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTAGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGc          >  2:531599‑M1/1‑150 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  1:1634769‑M2/1‑54 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     >  1:1634764‑M2/1‑54 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                     <  2:1634769‑M2/58‑5 (MQ=255)
                                                                                                                                                  gCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagggggcaagcccgcgaagaggcccatacaggcgccgctggacttggcctcgacaagcggttatctaatccttgcacagaccaaccacaggcc      <  2:982049‑M1/150‑134 (MQ=255)
                                                                                                                                                    TTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcgccg  <  1:592017‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5290231‑5290528

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: