Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,407 T→A 6.0% 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,293,4070TA 6.0% 665.5 / 21.9 214intergenic (‑91/+93)PP_4660/mdaBLysR family transcriptional regulator/drug activity modulator B
Reads supporting (aligned to +/- strand):  ref base T (83/113);  new base A (6/7);  total (89/125)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 216 total aligned reads displayed.

ATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5293258‑5293555
                                                                                                                                                     |                                                                                                                                                    
  gcatgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGc                                                                                                                                                    >  1:1359085‑M1/132‑150 (MQ=255)
  gcatgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgacgcccatctggacaaggtgtttggcaaggcacagtacttttatgtctgtaccGGCCTCTTCGCGGGCTTGc                                                                                                                                                    >  1:2371485‑M1/132‑150 (MQ=255)
  cGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGc                                                                                                                                                    <  1:750365/150‑1 (MQ=255)
  cGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGc                                                                                                                                                    <  2:5402882/150‑1 (MQ=255)
   gAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGcc                                                                                                                                                   >  2:1014750/1‑150 (MQ=255)
   gAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTgac                                                                                                                                                   >  1:3523743/1‑148 (MQ=255)
   gAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGGTGGTCTGTGCAAGGAGGAGATAACCGCTGGTCGAGGCCAAGTCCAGCGGCGCCTGGATGGGCCTCTTCGCGGGCGTGcc                                                                                                                                                   >  1:3885596/1‑150 (MQ=255)
    atgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  >  2:6053762‑M1/130‑150 (MQ=255)
    aGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  <  1:521664/150‑1 (MQ=255)
      gaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGc                                                                                                                                                >  2:2282751‑M1/128‑150 (MQ=255)
       cGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCt                                                                                                                                               >  2:516972/1‑150 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGcga                                                                                                                                              <  2:5889390/150‑3 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtggtgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                             >  2:3702701‑M1/125‑147 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTcc                                                                                                                                             >  1:2101229‑M1/125‑150 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                             >  2:864613‑M1/125‑147 (MQ=255)
         cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTcc                                                                                                                                             <  2:5895100/150‑1 (MQ=255)
         cAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCAAGCCCGcgaa                                                                                                                                             <  2:1949156/150‑4 (MQ=255)
            tCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          <  2:3677638/150‑1 (MQ=255)
            tCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          <  1:6231961/150‑1 (MQ=255)
            tCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          <  2:5778443/150‑1 (MQ=255)
            cctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                          >  2:525684‑M1/122‑150 (MQ=255)
             cGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                         >  2:4915396/1‑150 (MQ=255)
              tgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCGTCGCGGGGGTGGCcgcgccccccg                                                                                                                                        <  1:5995999‑M1/31‑5 (MQ=255)
              gATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                        <  1:256958/150‑1 (MQ=255)
              gATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                        <  2:178895/150‑1 (MQ=255)
              gATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                        <  2:1292055/150‑1 (MQ=255)
                   cGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACt                                                                                                                                   <  1:587085/150‑1 (MQ=255)
                   cGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCACCCCCAAGTACt                                                                                                                                   >  2:3760767/1‑150 (MQ=255)
                    ggCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTg                                                                                                                                  <  2:2452861/150‑1 (MQ=255)
                    ggCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTg                                                                                                                                  >  2:291206/1‑150 (MQ=255)
                    ggCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTg                                                                                                                                  <  2:2036831/150‑1 (MQ=255)
                        acgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGTattacgca                                                                                                                              <  1:4013937/150‑8 (MQ=255)
                        acgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCaaa                                                                                                                              <  1:6148275/150‑1 (MQ=255)
                        acgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCaaa                                                                                                                              <  2:6724721/150‑1 (MQ=255)
                        acgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCaaa                                                                                                                              <  1:5166408/150‑1 (MQ=255)
                          gaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAga                                                                                                                            >  2:4536243/1‑149 (MQ=255)
                              aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                        <  1:5486198/150‑1 (MQ=255)
                              aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                        <  1:4674421/150‑1 (MQ=255)
                              aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                        <  1:3479996/150‑1 (MQ=255)
                               ccTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTg                                                                                                                       <  1:660233/150‑1 (MQ=255)
                                cTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTACCACAAGTACTGCAAAGTCCTTGa                                                                                                                      >  1:2734133/1‑150 (MQ=255)
                                cTGGAGTTCGTCGGATCGGGTTTTCACATCAGGGCCGGGGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGa                                                                                                                      <  2:212259/150‑1 (MQ=255)
                                 tGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGaa                                                                                                                     <  1:20254/150‑1 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                               <  2:2806338/150‑1 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                               <  1:5234612/150‑1 (MQ=255)
                                        cttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCATCGGCGCCTGTATGGGCCTCTTCGCGGGCTAGCCCGCTCCCACAACTACAGCAAAGTCCATGAAAcaagcg                                                                                                              >  2:302121/1‑145 (MQ=255)
                                        cttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCg                                                                                                              >  2:6828880/1‑150 (MQ=255)
                                         ttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCAGCGt                                                                                                             >  2:1268418/1‑150 (MQ=255)
                                          tcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTa                                                                                                            <  2:6308620/150‑1 (MQ=255)
                                           cGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTaa                                                                                                           <  2:4594845/150‑1 (MQ=255)
                                           cGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTaa                                                                                                           <  2:1657787/150‑1 (MQ=255)
                                           cGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTaa                                                                                                           <  2:2171616/150‑1 (MQ=255)
                                              aTCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCGCCCCCCAATACAGCAAAGTCCTGTTAACCAGCGTAAAAc                                                                                                        >  1:1433035/1‑150 (MQ=255)
                                               tCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAACTCCTTGAAACCTGCGTAATAcc                                                                                                       >  2:2390969/1‑150 (MQ=255)
                                               tCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCGCTCCCCCAAGTACAGCAAAGTCCGTGAAACCGGCGTAATaac                                                                                                       >  2:1389318/1‑148 (MQ=255)
                                                cGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCt                                                                                                      <  1:6387716/150‑1 (MQ=255)
                                                 gggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtt                                                                                                     >  2:6419384/1‑149 (MQ=255)
                                                  ggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCGTGAAACCAGCGTAATACCagt                                                                                                    >  2:107383/1‑150 (MQ=255)
                                                  ggTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGATTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCACCAGTAAAGCCAAGTCCTTGAAAAATGCGGaaaaacagt                                                                                                    >  1:5643380/1‑143 (MQ=255)
                                                  cggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgt                                                                                                    >  2:404903‑M1/84‑150 (MQ=255)
                                                   ggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                   >  1:1739086‑M1/83‑150 (MQ=255)
                                                   gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCTGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                   <  1:4155683/150‑1 (MQ=255)
                                                   gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                   <  1:3326145/150‑1 (MQ=255)
                                                   gTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCGTGAAACCTGCGTAATACCggtg                                                                                                   >  2:4279651/1‑150 (MQ=255)
                                                    gatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTgg                                                                                                  >  2:1115102‑M1/82‑150 (MQ=255)
                                                     tttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTggg                                                                                                 >  2:5438016/1‑150 (MQ=255)
                                                       tCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCCAAGTCCTTGAAACCTGCGTAATACCTGTGGGAg                                                                                               >  2:4927165/1‑150 (MQ=255)
                                                        tgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtactGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGc                                                                                              >  1:6654971‑M1/78‑150 (MQ=255)
                                                         gccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCg                                                                                             >  1:4782951‑M1/77‑150 (MQ=255)
                                                         acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCg                                                                                             <  1:1014750/150‑1 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            <  1:2194199/150‑1 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            <  1:267900/150‑1 (MQ=255)
                                                              aGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCaa                                                                                        <  2:6211648/150‑1 (MQ=255)
                                                              aGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCaa                                                                                        <  2:5947602/150‑1 (MQ=255)
                                                              aGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCaa                                                                                        <  1:6124616/150‑1 (MQ=255)
                                                              aGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCaa                                                                                        <  1:4229404/150‑1 (MQ=255)
                                                                ggCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGc                                                                                      <  1:1832750/150‑1 (MQ=255)
                                                                ggCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGc                                                                                      <  1:3898247/150‑1 (MQ=255)
                                                                 gCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCCCGACCACAAGTACATCCAAGGCCATGAAAACAGCGTAATAACTGTGGGGGCGGGCAAGcc                                                                                     >  2:179892/1‑150 (MQ=255)
                                                                    tgtgGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCcgc                                                                                  >  2:460558/1‑150 (MQ=255)
                                                                    tgtgGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCcgc                                                                                  <  2:914831/150‑1 (MQ=255)
                                                                       ggttggtCTGTGCAAGGATTAGATAACCGTTTGTCGAGGCCAAGTCCTGGGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCGGCTCCCACAAGTACTGCAATGTCCTTGAAACCTGCGTAATACCTGTGGGTGCGGGCAAGCCCGCGaa                                                                               <  2:3240086/150‑1 (MQ=39)
                                                                             tCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCAGCGTAATAACTGTGGGAGCGGGCCAGCCCCCGAAGAGGcc                                                                         >  2:3591078/1‑150 (MQ=255)
                                                                                        aTTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCACAAGGCCGGCACAGGCCGTGAACCCGGCGGCAGACCGGGGGGCGCGGGCCAGCCCGCGACGAGGCCCGTCCAGAccg                                                              >  2:5341730/1‑148 (MQ=255)
                                                                                            gATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCACCCACAAGTACTGCAAAGTCCATGAAAACAGCGTAATACCAGGGGGAGCCGGCAAGCCCGCGAAGAGGCCCGTACAGACATaaaa                                                          >  2:4394101/1‑150 (MQ=255)
                                                                                                  cGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCGCCCACAAGGACAGCCCAGGCCTAGAAACCGGGGTAAAAACGGGGGGGGCCGGCCAGCCCCCCCCCAGGCCGGGTCCGACATAAAAATAcgg                                                    >  2:258602/1‑148 (MQ=255)
                                                                                                   gCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGGGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACtgt                                                   <  2:510713/150‑1 (MQ=255)
                                                                                                          gcggCCAAGGCCAGCGGCGCCTTTATGGGCCTCTTCGGGGGCTTGCCCGCTCCCACAAGTACTGCAACGTCATTGAAACGTGCGTAATACGTGTGGGAGGGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGc                                            <  2:4041098/148‑1 (MQ=38)
                                                                                                           aGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCACAAGTACTGCAAAGTCCTTTAAACCTGCGTAATAACGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGcc                                           >  2:4615265/1‑150 (MQ=255)
                                                                                                               cAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAac                                       >  2:6155994/1‑150 (MQ=255)
                                                                                                                           cgcCTTTATGGGCCTCTCCGGGGGCTTGCCCGCTCCCAAAAGTACTGAAAAGTCCTTTAAACCTGCTTACTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGa                           <  2:3530318/150‑1 (MQ=38)
                                                                                                                           cgcCGGTAGGGGCCCCTCCGCGGGCTGGCCCGCCCCCCCAAGAACTCCAACCCCCTTAAACCCTGCGAAATACCGGTGGGCGCGGGAACGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGa                           <  1:3188169/150‑1 (MQ=255)
                                                                                                                              cTGTAGGGGCCTCTGCGGGGGCTGGCCCGCTCCCACATGTACCGAAACGGCCTTCAACCCTGGGTAATACCTGTGGGAGGGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATgg                        <  2:98605/150‑1 (MQ=37)
                                                                                                                                    gggCCTCTTCGGGGGCTTGCCCGCTCCCACACTTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGAAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCa                  <  2:3523105/150‑1 (MQ=255)
                                                                                                                                     ggCCTCTGGGGGGGCGTGCCGGCCCCCACCAGTACTGCAACGCCGTGCAAACCTGGTTATTCCCGGTGGGTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTCGCCAAACACCTTGTCCAGATGGGCCTCAt                 <  2:1995149/150‑1 (MQ=255)
                                                                                                                                      GCCTCTTCGCGGGCTTGCCCGATCCCAAAATTACTGAAAAGTCCTTTAAACCTGCGAAATACCTGTGGTAGCGGGAAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATa                >  2:6405901‑M1/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCTCCCAAAAGTACTTCAAATTGCTTTAAAACAGCGTAATAACTGTTGGATCGGGGAAGCCCGCGAAGAGGCCGGTAAAGAAATAAAATTAATTTGCCTTTACAAAAAACTTTTACAGATTGGGCTAATAAGCTg           >  2:790416/1‑150 (MQ=255)
                                                                                                                                              GCGGGCTTGCCCGCTCCCACAACTACAGCAAACTCCTTGAAACCAGCCTAATAACTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTAAAGACAGAAAATTACAGTGCCAAGCCAAACACCTTGTCCACATGGGCCTCAAAAGCTGCCa        <  1:1473153‑M1/150‑1 (MQ=255)
                                                                                                                                               CGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCaa       <  1:4460793‑M1/150‑1 (MQ=255)
                                                                                                                                                    TTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcgccg  <  2:5486402‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5293258‑5293555

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: