Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,407 T→A 5.1% 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 5.1% 603.0 / 17.3 197intergenic (‑91/+93)PP_4660/mdaBLysR family transcriptional regulator/drug activity modulator B
Reads supporting (aligned to +/- strand):  ref base T (95/88);  new base A (8/2);  total (103/94)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 198 total aligned reads displayed.

ATCGAGCCGCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAACGCCG  >  NC_002947/5293258‑5293555
                                                                                                                                                      |                                                                                                                                                     
gggcatgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCtt                                                                                                                                                       >  2:811642‑M1/134‑150 (MQ=255)
  gcatgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattgtatgtctgtaccGGCCTCTTCGCGGGCTTGc                                                                                                                                                     <  2:5116608‑M1/19‑1 (MQ=255)
     tgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCg                                                                                                                                                  <  2:2401051‑M1/22‑1 (MQ=255)
     tgaccggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCGTGCCCg                                                                                                                                                  <  1:2746562‑M1/22‑1 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                               >  2:4679121‑M1/1‑150 (MQ=255)
        gCAGTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                                                                               <  1:4269736‑M1/150‑1 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                              >  2:5150812‑M1/125‑147 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                              >  1:4424015‑M1/125‑147 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                              >  1:3029329‑M1/125‑147 (MQ=255)
         cggcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagcaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                              >  2:3144526‑M1/125‑147 (MQ=255)
         cggcctgccgaccgtcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAGCccgcgaa                                                                                                                                              >  2:102563‑M1/125‑147 (MQ=255)
           gcctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCAAgcccgcgaaga                                                                                                                                            >  2:1300892‑M1/123‑145 (MQ=255)
           gTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCA                                                                                                                                            <  1:1139329‑M1/150‑1 (MQ=255)
           gTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCA                                                                                                                                            <  1:5300410‑M1/150‑1 (MQ=255)
           gTCGATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCA                                                                                                                                            <  1:3526526‑M1/150‑1 (MQ=255)
            cctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCcac                                                                                                                                           >  1:979672‑M1/122‑150 (MQ=255)
             ctgccgaccttcctggctgtggatgtgatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                          >  1:252562‑M1/121‑150 (MQ=255)
              gATGACGGCGACGAATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACaa                                                                                                                                         <  1:2257364/150‑1 (MQ=255)
                        acgaATACCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCGCTCCCACAAGTACTGCaaa                                                                                                                               >  2:1101787/1‑150 (MQ=255)
                              aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                         <  1:2670535/150‑1 (MQ=255)
                              aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                         <  1:28882/150‑1 (MQ=255)
                              aCCTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCtt                                                                                                                         <  1:3734865/150‑1 (MQ=255)
                                cTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGa                                                                                                                      >  2:3946687/1‑150 (MQ=255)
                                cTGGAGTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCCCGCCCCCCAGTACTGCAAACTccaggg                                                                                                                       >  1:4452055/1‑146 (MQ=255)
                                     gTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCt                                                                                                                 <  2:3380253/150‑1 (MQ=255)
                                     gTTCTTCGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCt                                                                                                                 <  1:1295748/150‑1 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                               >  2:2746562/1‑150 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCCCAACCACAACTAAAGCAAAGTCCGTTAAAccagg                                                                                                               >  1:2401051/1‑147 (MQ=255)
                                       tcttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGACCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                               <  1:3443176/150‑1 (MQ=255)
                                       gatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                               >  2:2255238‑M1/95‑150 (MQ=255)
                                       gatgaagcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGc                                                                                                               <  2:3450277‑M1/56‑1 (MQ=255)
                                        cttcGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCg                                                                                                              <  1:3177389/150‑1 (MQ=255)
                                           cGGATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTaa                                                                                                           <  2:2612830/150‑1 (MQ=255)
                                            agcgcccggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAAt                                                                                                          >  1:1725581‑M1/90‑150 (MQ=255)
                                             gATCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACAGCAAACTCCTTGAAAACGGCGGAATa                                                                                                         >  2:1310459/1‑150 (MQ=39)
                                              aTCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATAc                                                                                                        <  1:1338728/150‑1 (MQ=255)
                                               tCGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCCCGCCCACAAGGACGGCAAAGTCCGGGAAACCAGCCGAATAcc                                                                                                       >  1:5116616/1‑150 (MQ=255)
                                                cGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCGCTCCCACAAGAACAGCAAACTCCATGAAAACAGCGTaaaaacg                                                                                                      >  2:25601/1‑145 (MQ=255)
                                                cGGGTTTTCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGAGTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCCCGCCCACAAGGACTGCAAAGTCCGGGAAACCAGCGGAATAccg                                                                                                      >  1:1367781/1‑149 (MQ=255)
                                                   ggatgtgccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtg                                                                                                   >  1:2995633‑M1/83‑150 (MQ=255)
                                                    ttttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgg                                                                                                  >  1:3060302‑M2/1‑145 (MQ=255)
                                                      ttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtggga                                                                                                >  1:1030370‑M2/1‑143 (MQ=255)
                                                      ttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtggga                                                                                                >  1:1030144‑M2/1‑143 (MQ=255)
                                                      ttCACATCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtggga                                                                                                >  1:5085388‑M2/1‑143 (MQ=255)
                                                         acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCGTGAAACCAGCGTAATACCggtgggggcg                                                                                             <  2:283738‑M2/150‑11 (MQ=255)
                                                         acaTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCTCCCACAAGTACAGCAAAGTCCATGAAAACATCGTAATAACagtgggggcg                                                                                             <  1:1341421‑M2/150‑11 (MQ=255)
                                                         acaTCAGGGCCGGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCACAAGTACAGCAAAGTCCGTGAAACCAGCGTAAAACCtgggggagcg                                                                                             <  2:612870‑M2/150‑11 (MQ=255)
                                                          ccggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCgg                                                                                            >  2:4035859‑M1/76‑150 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTGCCCCACCCCACAAAACAGCAAAATACATGAAACCAGCGGAATAAACGTGGGGGCgg                                                                                            >  1:1111297/1‑150 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgg                                                                                             <  1:2844683‑M1/150‑46 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgg                                                                                             <  1:4163857‑M1/150‑46 (MQ=255)
                                                          caTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgg                                                                                            >  1:2242382‑M2/1‑139 (MQ=255)
                                                          caTCAGGGCCTGTGGCTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgg                                                                                            >  2:212273‑M2/1‑139 (MQ=255)
                                                           cggcggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCggg                                                                                           >  1:1203874‑M1/75‑150 (MQ=255)
                                                           aTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcggg                                                                                            <  1:1923315‑M1/150‑47 (MQ=255)
                                                           aTCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCGGCGGAATACCtgtgggagcggg                                                                                           <  2:5121004‑M2/150‑13 (MQ=255)
                                                            tCAGGGCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgggc                                                                                           <  2:5355384‑M1/150‑48 (MQ=255)
                                                              cggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCaa                                                                                        >  1:518‑M1/72‑150 (MQ=255)
                                                              cggcgttggcagcttatgaggcccatctggacaaggtgtttggcaaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCaa                                                                                        >  2:1919667‑M1/72‑150 (MQ=255)
                                                                ggCCTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgggcaagc                                                                                      >  2:2402801‑M2/1‑133 (MQ=255)
                                                                  ccTGTGGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgggcaagccc                                                                                    >  1:847713‑M2/1‑131 (MQ=255)
                                                                     gtgGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCCCACCCACAACTACTGCCAAGTCCGTGAAACCAGCGGAATAACggtgggagcgggccatcccgcg                                                                                 <  2:1526865‑M2/150‑23 (MQ=255)
                                                                     gtgGTTGGTCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCCCGCCCCAAAGGACTGCAAAGACCGTGAAAACAGCCGAATAACAGGGGGTGCGGGCAAccccccc                                                                                 >  1:1238206/1‑147 (MQ=255)
                                                                      tggttggtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATAc                                                                                                        <  1:4010389/126‑1 (MQ=255)
                                                                      tggttggtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATAc                                                                                                        >  2:4010389/1‑126 (MQ=255)
                                                                        cgcttccgaggcccatctggacacggtgtttggcaaggcacagtaattttatgactgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACACGTACTGCACAGTCCTTGAACCCTGCGTAAGACCTGTGGGAGCGGGCAAGCCCGCGAag                                                                              >  2:3954946‑M1/62‑150 (MQ=255)
                                                                         ttgttCTGTGCATGTATTACTTAAGCTCTTTTGGCGGCCAAGTCCTGCGGCGCCTGTTTGGGCATCTTGGGGGGTTTGCCCGCTCCCACAAGTACTGAAATGTCCTTTAAACCTGCTTAATACCTGTGGTTGCGGGAAAGCCCGCGAaga                                                                             <  2:220734/150‑1 (MQ=255)
                                                                         ttggtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTTCCCGCTCCCAAAAGTAATGCAAAGTACTTGAAACCTGCGTAATACCagtgggagcgggcaaccccgcggaga                                                                             <  2:5174935‑M2/150‑27 (MQ=255)
                                                                         gcttatgaggcccatctggacacggtgttgggccaggcacagtaattttatgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCCCAAGTACTGCAAAGTCCTTTAAACCTGCGTAATACCTGTGGGCGCGGGCACGCCCGCGAaga                                                                             >  2:1677592‑M1/61‑150 (MQ=255)
                                                                            tatgtggcccatctggacaagttttttggctaggcacagtaattttttgtctgtaccGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTTCAAAGTCATTTAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGc                                                                          >  2:1354738‑M1/58‑150 (MQ=255)
                                                                            gtCTGTGCATGTTTTATATACCGGGTTGTCGAGGCCATGTCCTGCGGGGCTTGTTTGGGCTTCTTCGCGGGCTTGCCCGCTCCCACACTTACTGCAATGTCTTTTAAACCTGCTTATTACGTGTGGGAGCGGGCAAGCCCGCGAAGAGGc                                                                          <  2:4019270/150‑1 (MQ=255)
                                                                            gtCTGTGCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtgggagcgggcaagcccgcgaagaggc                                                                          >  2:4153164‑M2/1‑121 (MQ=255)
                                                                                  gCAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggactttgcagtacttgtgggagcgggc                                                                                           >  1:3785896‑M1/1‑81 (MQ=255)
                                                                                  gCAAGGATTAGAGAACCGATTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCCTGCCCGCACCCCAAATTACTGCACACTCCATGCAACCGGCGGAATACCTGTGGGAGCGGGCAAGCCCCCGAAGAGGcccggac                                                                    >  2:3584950/1‑145 (MQ=255)
                                                                                   cAAGGATTAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCCCAAGTACTGCCAAGTCCGTGGAACCAGCGGAATACCGGTGGGGGCGGGCCAGCCCCCGCAGAGGCCGGTACa                                                                   >  2:3929423/1‑150 (MQ=255)
                                                                                      gtaTTACATACCGGCTGGTCTCGGCCAGGTCACGGGGCGCCTGTATGGGCCCGTTCGGGGGCTTGCCCGCCCCCAAACGTACTGAAACGTCGTTGAACCGTGCTTAATACCTGTGGGTGGGGGCACGCCCGCGAAGAGGCCGGTACAGAc                                                                <  2:4408377/148‑1 (MQ=255)
                                                                                         ttAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTTCCCGCTCCCCCAAGTACTGCAAAGTCCGTGAAACCAGCGTAATAACAGTGGGAGCCGGCCAGCCCCCCAAGAGGCCCGTACAGACaaa                                                             >  1:5609855/1‑148 (MQ=255)
                                                                                         ttAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCGCGCCCACAAGGACTGCAAAGTCCGTGAAAACGGGGTTAAACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATa                                                             >  1:5472499/1‑150 (MQ=37)
                                                                                          tAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCCAAGTCCTTGAAACCTGCGGAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATaa                                                            >  2:1853477/1‑150 (MQ=255)
                                                                                          tAGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATaa                                                            >  2:3657792/1‑150 (MQ=255)
                                                                                           aGATAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCGCCCACAATTACAGCAAAGTCCTTGAAACCTGCGTAATAACTGTGGGAGCGGGCAAGCCCCCGAAGATGCCGGTAAAGACAaaaa                                                           >  2:5428267/1‑150 (MQ=38)
                                                                                              tAACCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCATGCCCCCTCCCACAAATACTGCCAAGTCCGTGAAACCGGCGGAAAACCAGTGGGGGCGGGCAAGCCCCCGAAGAGGCCGGTAAACACAGaaaaat                                                        >  1:4659448/1‑148 (MQ=255)
                                                                                               cacccgttGTCTCGGCCACGCCCAGGGGCGCGTGTAGGGGCCTCGTCGGGGGCTGGCCCGCTCCCACAAGTACTGAAATGTCCTGGAAACCTGGTTAATACATGTGGGTGCGGGCAAGCCGGCGAAGAGGCCGGTACAGAAATAAAATTa                                                       <  2:2891007/144‑1 (MQ=255)
                                                                                               acccGCTTGTCGCGGCAACGCCCAGCGGCGCCTGTAGGGGCCTCTCCGCGGGCTTGCCCGCTCCCACAAGTACTCCAAAGTCCTTCAAACCTGCTTAACACGTGTGGGTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACATAAAATTa                                                       <  2:3848874/148‑1 (MQ=255)
                                                                                               aaCCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCCCAGCCACCAGTAATGCAAAGTCCTGTGAACCGGCGTGATACCGGTGGGCGCGGGCCACCCCGCGCAGAGGCCGGTACCGACATAAAAATa                                                       >  2:2507545/1‑150 (MQ=255)
                                                                                                aCCGCTTGTCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCGTGCCCCCGCCCAAAAGTACAGCAAAGTGCTTGTAAACATGGGTATAACAGTGGGGGCGGGGCAGCCCGCGAAGAGGCCGGTAAAGAAATAAAATTaa                                                      >  2:2382770/1‑149 (MQ=255)
                                                                                                 ccGCTTGTCGAGGCCAAGTCCTGGGGCGCCTGTATGGGCCTCTTGGCGGGCTTGCCCGCTCCCACATGTAGTGCAAAGTCATTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACt                                                     <  2:4101685/150‑1 (MQ=37)
                                                                                                        tCGAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACCAGGACAGCAAAGTCCTTGTAACCAGCGTAATACCAGTGGGAGCGGGCAAGCCCCCCAAGAGGCCCGTACAGACATAAAAATAAAGTGCCtt                                              >  1:3236543/1‑150 (MQ=255)
                                                                                                          gAGGCCACGTCCTGCGGCGCCTTTTTGGGCTTCTTGGGGGGCTTGCCCGCTCCAACACTTACTTCAACGGCCTTGAAACCTGGTTAATACGTGTGGGAGGGGCCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGc                                            <  2:1517056/150‑1 (MQ=255)
                                                                                                          gAGGCCAAGTCCAGCGGCGCCTGTATGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCtgtggga                                                                                                <  1:727791‑M2/98‑8 (MQ=255)
                                                                                                                               tGTACGGGCCTCGTCGCGGGCTTGCCCGCTCCCACCggtattccgcaggtttcaaggactttgctgtacttgggggagcgggcaagcccgcgaagaggcccatacaggcgccgctggacttggcctcgacaagcggttatctaatccttg                        <  2:3929367‑M1/150‑115 (MQ=255)
                                                                                                                                 tATGGGCCTCTTCGCGGGCGTGCCCCCTCCCACAAGGACTGCAAAGTCCAGGTAACCGGCGTCATACCTGGGGGAGCGGGCACGCCCCCGCAGAGGCCGGGACAGCCAGAAAATTACAGTGCCTTGCCAACCACCGTGTCCACATGGGcc                     >  2:3228105/1‑150 (MQ=255)
                                                                                                                                     ggCCTCTTCGCGGGCTTGCCCCCCCCCACAAATAATTCAAAGTACATGAAACAAGCCGAAAACAAGGGGGCGCGGGCAACCCCCCCAAGAAGCCGGTAAAGAAAAAAAATTACTGTGCCTTTCCAAAAACCGTTTACAGAAGGGccaaat                 >  2:1970508/1‑146 (MQ=255)
                                                                                                                                         TCTTCGCGGGCATGCCCGCTCCCACAAGTACTGCAAAGTCCTTGTAACCTGCGGAATACCTGTGGGAGCGGGCAAGCCCGCGAAGATGCCGGTACAGACATACAATTACTGTGCCTTGCCAAACACCGTGTCCAGATGGGCCGCATAAGc             <  2:2220048‑M1/150‑1 (MQ=255)
                                                                                                                                           ggaGCGGGCAAGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTg           <  1:5511070‑M1/147‑1 (MQ=255)
                                                                                                                                           TTCGCGGGCTTGCCCGCTCCCACAAGTACAGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTg           <  2:1368686‑M1/150‑1 (MQ=255)
                                                                                                                                             cgcgGGCTTGCCCGCTCCCACAggtattacgcaggtttcaaggacttttcagtacttgtgggagcgggcaagcccccgcagaggcccatacaggcgccgcaggacttggccgcgacaagcggttatctaatccttgcacagaccaaccac          >  2:2592522‑M1/1‑22 (MQ=255)
                                                                                                                                              gggggCGTGCCCGGACCCACAAGTACTGCAACGCCCTGGAACCCTGCTTTCTCCGGGGGGGCGCGGGCACGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCa        <  2:3780935/148‑1 (MQ=255)
                                                                                                                                                GGGCTTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAc      <  1:1551659‑M1/150‑1 (MQ=255)
                                                                                                                                                     TTGCCCGCTCCCACAAGTACTGCAAAGTCCTTGAAACCTGCGTAATACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGTACAGACATAAAATTACTGTGCCTTGCCAAACACCTTGTCCAGATGGGCCTCATAAGCTGCCAAcgccg  <  2:4289120‑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 < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: