Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,393 T→G 8.2% intergenic (‑231/+153) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,3930TG 8.2% 336.4 / 2.1 120intergenic (‑231/+153)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (36/70);  new base G (11/2);  total (48/72)
Fisher's exact test for biased strand distribution p-value = 6.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGC  >  NC_002947/4081244‑4081542
                                                                                                                                                     |                                                                                                                                                     
ggggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCt                                                                                                                                                       <  1:1817714/150‑1 (MQ=255)
 tcgacagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCtt                                                                                                                                                      <  2:4951805‑M1/15‑1 (MQ=255)
  ggCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTg                                                                                                                                                     <  1:1591834/150‑1 (MQ=255)
  cgacagcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTg                                                                                                                                                     >  1:4449145‑M1/135‑150 (MQ=255)
   gcgcGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGc                                                                                                                                                    <  2:4830216/150‑1 (MQ=255)
      agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                 >  1:1944110‑M1/131‑150 (MQ=255)
      agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTc                                                                                                                                                 >  1:1944183‑M1/131‑150 (MQ=255)
      agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCGTgcggc                                                                                                                                                 <  2:636619‑M1/20‑3 (MQ=255)
      agcggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCGTGCGTc                                                                                                                                                 <  2:1558570‑M1/20‑1 (MQ=255)
        cggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaccgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                               >  1:3302406‑M1/129‑150 (MQ=255)
        cggccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                               >  1:4709358‑M1/129‑150 (MQ=255)
        aTTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                               <  1:237448/150‑1 (MQ=255)
        aTTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                                               <  2:2239789/150‑1 (MQ=255)
          tgtgCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGa                                                                                                                                             <  2:4683607/150‑1 (MQ=255)
           ccgtgtgtgcatgtccggcctgcgcccggccaacctgcaacgggttgccgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctccagtggggcgcggttggacttTGTGGGAGCGGGGTGgcgtcgggcg                                                                                                                                            <  1:3835455‑M1/25‑5 (MQ=255)
               aCGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAggg                                                                                                                                        >  2:3178841/1‑150 (MQ=255)
                cGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                       <  2:3890696/150‑1 (MQ=255)
                cGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                       <  2:3948947/150‑1 (MQ=255)
                 ggAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                      <  1:1193469/150‑1 (MQ=255)
                      ggTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                 <  1:4861378/150‑1 (MQ=255)
                       gTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGTCGGGAGGGGGCCGCAAAg                                                                                                                                >  1:319110/1‑150 (MQ=255)
                       gTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGAGAGGGCCGCaaac                                                                                                                                >  2:3777463/1‑149 (MQ=255)
                        ttGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                               >  1:1960907/1‑150 (MQ=255)
                        ttGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGACGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                               <  2:929520/150‑1 (MQ=255)
                           cACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                            <  2:1120823/150‑1 (MQ=255)
                             cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGccc                                                                                                                          <  1:5105748/150‑1 (MQ=255)
                             cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGccc                                                                                                                          <  1:5105783/150‑1 (MQ=255)
                              ttGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTGGCGGCGCGGTAGGGGCGCAAAGCGGgccc                                                                                                                         >  2:3569791/1‑150 (MQ=255)
                               tGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCa                                                                                                                        <  1:1435189/150‑1 (MQ=255)
                               tGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCCGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCa                                                                                                                        <  1:3479842/150‑1 (MQ=255)
                                gCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                       >  2:5131942/1‑134 (MQ=255)
                                gCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                       <  1:5131942/134‑1 (MQ=255)
                                gCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAg                                                                                                                       <  1:891248/150‑1 (MQ=255)
                                 cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGGCGCGGTAGGGCCGCCAAGCGGccccccc                                                                                                                      >  1:4693690/1‑147 (MQ=255)
                                 cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGccctatc                                                                                                                      <  1:3747962/150‑5 (MQ=255)
                                 cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                      <  1:4563170/150‑5 (MQ=255)
                                 cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                      <  1:4563219/150‑5 (MQ=255)
                                 cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                      <  2:3725712/150‑5 (MQ=255)
                                   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCaa                                                                                                                    <  1:5093647/150‑1 (MQ=255)
                                   aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCaa                                                                                                                    <  1:4887526/150‑1 (MQ=255)
                                      ggacaacctgcaacgggttggcgaagcgttcgctgccgttcagaagtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAttt                                                                                                                 <  2:2095323‑M1/52‑1 (MQ=255)
                                       ccTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTc                                                                                                                <  2:2539031/150‑1 (MQ=255)
                                         ttttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAg                                                                                                              <  2:4297187/150‑1 (MQ=255)
                                          tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                             <  1:4902720/150‑1 (MQ=255)
                                          tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                             <  2:3720456/150‑1 (MQ=255)
                                             ccGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGc                                                                                                          <  1:5138819/150‑1 (MQ=255)
                                               gaggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAAATCAGCGGCaa                                                                                                        >  2:2712324/1‑150 (MQ=255)
                                               gaggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAACCGGCCCCACCAAATTAAGCGGCaa                                                                                                        >  2:4379814/1‑150 (MQ=39)
                                                aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaaa                                                                                                       <  2:4574490/150‑1 (MQ=255)
                                                aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaaa                                                                                                       <  1:3097282/150‑1 (MQ=255)
                                                aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaaa                                                                                                       <  1:4005692/150‑1 (MQ=255)
                                                   agCGGTTGCACTTGTGTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCt                                                                                                    <  1:3081156/150‑1 (MQ=255)
                                                     cGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGc                                                                                                  <  1:2630739/150‑1 (MQ=255)
                                                       gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGc                                                                                                          <  1:675874/140‑1 (MQ=255)
                                                         tGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATc                                                                                              <  2:776883/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGGGTGGGGGGGGGTTGGGGCGCAAAGCGGGGCCCGCAATTTGAGCGGCAAAGGTGGAAtgttt                                                                                           >  1:3622505/1‑146 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  1:3480912/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  2:4024644/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                           <  1:242492/150‑1 (MQ=255)
                                                            aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCCCACAGCGCCCCCAGCAAATACAGCCGCCAAGCGGCAAACgga                                                                                           >  1:3532844/1‑148 (MQ=255)
                                                              ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                         <  1:3494773/150‑1 (MQ=255)
                                                               tGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAggg                                                                                        <  2:2202320/150‑1 (MQ=255)
                                                               tGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAggg                                                                                        <  1:28289/150‑1 (MQ=255)
                                                                gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATa                                                                                                                                           <  1:2112212/98‑1 (MQ=255)
                                                                gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgggg                                                                                       <  2:2869477/150‑1 (MQ=255)
                                                                gTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTAGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCCATCGTAgggg                                                                                       >  2:5364129/1‑150 (MQ=255)
                                                                 tCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTCAGCGGCAAAGCTGCAATCGTAGGGgc                                                                                      >  2:4971197/1‑150 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGGGGGGCGGGATGGGCGCAGAGCGGGCCCCGCAATAAAAGCGGCAAAGCGGCAATAATAGGGgcgc                                                                                    >  1:1558570/1‑150 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCAGCAATTTAAGCGGCAAAGCTGCAAACGGAGGGgcgc                                                                                    >  2:2152216/1‑150 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                               <  2:2345583/107‑1 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcc                                                                                                                           >  2:3802403/1‑111 (MQ=255)
                                                                   tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGgcgc                                                                                    <  1:927047/150‑1 (MQ=255)
                                                                       cAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGGGCGTGGGGCTAGGGGCGCCACGCGGCCCCCCCCAATTACGCGGCAAAGCAGCCAACGTAGGGGGGCGccc                                                                                >  1:4360446/1‑149 (MQ=255)
                                                                          gCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGCGAGGGCCGCCAACCCGCCCCACCAAAATCAGCGGCAAAGCTGCAAGCGTAGGGGCGCTCCGCAc                                                                             >  2:620755/1‑150 (MQ=255)
                                                                            ggAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                           <  2:4882577/150‑1 (MQ=255)
                                                                            ggAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                           <  2:4882686/150‑1 (MQ=255)
                                                                                  cATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTCCTAGGCAAATGGCCGCCAAACGGCCCCAGCAAATTAAGAGGCAAAGCTGCAAAAGTAGGGGCCCTCCGCACCCctttag                                                                     >  2:1996635/1‑148 (MQ=255)
                                                                                  agtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcg                                                                     >  1:4851968‑M1/55‑150 (MQ=255)
                                                                                   gtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcgc                                                                    >  2:480695‑M1/54‑150 (MQ=255)
                                                                                   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCtt                                                                         <  2:2754244/145‑1 (MQ=255)
                                                                                   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCtt                                                                         >  1:2754244/1‑145 (MQ=255)
                                                                                   aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCACCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCACCGCACCCATTTcgc                                                                    >  2:4407793/1‑150 (MQ=255)
                                                                                      gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCATGCGTCGCGATAGGGCCGCAAACCGGCCCCACCAAATTAATCGGCAAAGCTGCAAAAGTAGGGGCGCACCGCACCCTTTTAGCGAc                                                                 >  2:2172539/1‑150 (MQ=37)
                                                                                      actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                 >  1:2424032‑M1/51‑150 (MQ=255)
                                                                                        tGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                           >  1:3273146/1‑138 (MQ=255)
                                                                                        tGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                           <  2:3273146/138‑1 (MQ=255)
                                                                                            gtCGATCGGTGGAATGGAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGTTAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                           <  2:1398532/150‑1 (MQ=255)
                                                                                            caacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                           >  1:2095323‑M1/45‑150 (MQ=255)
                                                                                                        aaTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCCAGAGGGCCCCAAACCCGCCCCACCCAATTCAGCGGCAAAGCGGCAAACGGAGGGGCGCGCCCCAACCGGTTCGCGCCGCCACGCCGCGCCCACaa                                               >  2:982200/1‑150 (MQ=255)
                                                                                                               cggcgggggCGGCTGGGGGGGGGCGGGGGGGCGCGGCGTGGCGTCGCGATCGGGCCGCAACGCGGCCCCCGCATTTTCCGCGGCAAAGCTGCTTTGGTTGGGGCGCTCGGCCCCCTTTTCGCGACGCACGGCCGCTCCCACAAAGTCCaa                                        <  2:4642097/144‑1 (MQ=255)
                                                                                                                           cTGTGGGGCGGCCGGTGGGCGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCAGCAATTTCCGCGGAAACGCTGCAATCGTCGGGGCGCTCCGCACCCTTTCCGCGCCGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCcaca                            <  2:407205/150‑1 (MQ=255)
                                                                                                                                 ggcggcCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGAAACGCGGCCCCAGCAATTTCCGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAggg                      <  2:4394320/150‑1 (MQ=255)
                                                                                                                                    ggcCTGGGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGGGGCCCCAGCATTTTCCGCGGCAAAGCGGCAATCGTAGGGGCGCGCCGCCCCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGcc                   <  2:793216/150‑1 (MQ=37)
                                                                                                                                          tgGGCGCGGCCTTGCGGCGCGAACGGGCCGCAACGCGGCCCCAGCAATTTCCGCGGCAAAGCGGAAAAGGTCGGGGCGCGCCGCACCCTTTGCGCCCCGCAAGCCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt             <  2:54129/150‑1 (MQ=255)
                                                                                                                                           GGGCGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCCGCATTGTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTGCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTg            >  2:1584353‑M1/1‑150 (MQ=255)
                                                                                                                                               GCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCAGCAATTTCAGCGGCAATGCTGCAATCGTAGGGGCGCACCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGccc        >  2:3466401‑M1/1‑150 (MQ=255)
                                                                                                                                                 GGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGt      <  1:3109105‑M1/150‑1 (MQ=255)
                                                                                                                                                   CCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTg    <  2:3628246‑M1/150‑1 (MQ=255)
                                                                                                                                                    CTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg   <  1:38877‑M1/150‑1 (MQ=255)
                                                                                                                                                     TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTCCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc  >  2:4064833‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGC  >  NC_002947/4081244‑4081542

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: