Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,925 G→A G52G (GGG→GGA mprA → DNA‑binding transcriptional regulator MprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,9250GA90.6% 174.1 / 6.9 64G52G (GGG→GGAmprADNA‑binding transcriptional regulator MprA
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (24/34);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

ACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAA  >  NC_000913/2810679‑2811170
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGTTGAAATTAATGAGGTCATACTCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTa                                                                                                                                                                                                                                                   <  1:841850/251‑1 (MQ=255)
      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCACTAAAATGCTGAAGGCTCAGGGAATTAACg                                                                                                                                                                                                                                                >  2:167560/1‑248 (MQ=255)
           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACgttg                                                                                                                                                                                                                                        >  1:32665/1‑251 (MQ=255)
                ttgcTAAAACAATCTATCTACTGTATTTGTCCATGTCGTTACTATATCGTATGAAATTAATGAGGAAATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTGCGCTCCAGCCGCCACGAATATTTTCCTTATCAGGAGCTCCTTCTGACTCGTATTTGCATGCACATTCAATGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGTAATTAACGAGTCgttttcta                                                                                                                                                                                                                                    <  2:369561/247‑6 (MQ=255)
                       aaCAATACATGTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTCCGTTTACGCCCATTGAACAAATGCTAAAATTGCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTg                                                                                                                                                                                                                            <  2:523537/251‑1 (MQ=255)
                        aCAATACATTTATTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGATGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAGTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGa                                                                                                                                                                                                                           <  2:659342/251‑1 (MQ=255)
                         caattaattgaCTTTATTTTTCACTGTCGTTACTATATCGGTTGAAATTAATGAGGTCATTCCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAGCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGAt                                                                                                                                                                                                                          <  2:387118/241‑1 (MQ=255)
                               aTTTACTTTATTTTGCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGATCAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGATAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCt                                                                                                                                                                                                                    <  2:416822/251‑1 (MQ=255)
                                   aCTTTTTTTGTCACTGTCGTTACTATATCGTCTGAAATTAATGAGGACATACCCAAATGGATATTTCGTTTCCGCCCATTGAACAAATGCTAAAATTCCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                                                                                                                                                                                                 <  2:729977/250‑1 (MQ=255)
                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtctc                                                                                                                                                                                                            >  1:703044/1‑249 (MQ=255)
                                          ttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg                                                                                                                                                                                                         >  1:499744/1‑251 (MQ=255)
                                              cACTGTCGTTACTTTATCGGCTGAAATTAATGAGGTCATACCCATATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCAAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaa                                                                                                                                                                                                      >  2:155986/1‑250 (MQ=255)
                                                cTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCAGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTACCGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAcc                                                                                                                                                                                                   >  2:62122/1‑251 (MQ=255)
                                                cTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTACTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGGGTCTAAGGAAAAcc                                                                                                                                                                                                   >  2:62154/1‑251 (MQ=255)
                                                            tatCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTAAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGc                                                                                                                                                                                       >  1:118630/1‑251 (MQ=255)
                                                               cggttgaatagaatgagtgaaTACCCAAATGGATATTTTGTTTCCGCCCATTGACAATATGTTAAAATTTCGCGCCAGCCGCCACGAAGAGTTTCCTTATCATGAGATCCGTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTAGGGCGTTGATTACGCTGGAGTCTCAGGAAAACCTCAGTATTCAGCCt                                                                                                                                                                                     <  2:527040/230‑1 (MQ=255)
                                                                 gcggAAATTAATGAGTTCATCCCCAATTGGATAGTTCGTTTACGCCCAGTGAAAAAATGATAAAATTTCGCGCCATCCGCAACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCTATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCt                                                                                                                                                                                  <  2:746284/248‑1 (MQ=255)
                                                                     aaaTTAATGAGGGCATACCCAAATGAATAATTCGTTTACGCCCCTTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCATTATCAGTAGATCCTTATGACTCGTCTTTGCATGCACATGCAAAGCAATCTGCTGGAGAACCGCAATAAAAGGCTGAAGGCTCAGGGAATTAAAGAGACGTTGTTTATGGCGTTGAGTAAACTGGAGTCTCAGACAAACCACAGTATTCAGCATTCTGact                                                                                                                                                                              >  2:757649/1‑249 (MQ=255)
                                                                      aaTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAAtt                                                                                                                                                                             >  1:15317/1‑251 (MQ=255)
                                                                                   attcccataggggtAGTCCGTTTACGCCCATTTACCAAATACTAAAATGGCGCGACAGCCGCCACGAGGATTTTCCTTATCAGGAGATCCTTCTGTCTCGTCTATGCATTCGCATGTAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTTCAGGCTCAGGTAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCGCAGGAAAACCTCAGTATTCAGCCTTTTGAATTAAGGTGGGCTCtt                                                                                                                                                                <  2:533232/238‑1 (MQ=255)
                                                                                   aTACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGAAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCtt                                                                                                                                                                <  1:529066/251‑1 (MQ=255)
                                                                                    tACCCAAATGGTTAGTTCGTTTCCGCAAATTGGACACATGATAAAATCTCGCGCCAGCCGCCACGAAGATTTTCCTTATACGAAGATCCTTCTGCCTCGTCTTTGCATGCACATGCAATGCAAGATGGTGGCAAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAAAGAGACGTTGTTTATGGCGTTGATTACGCTAGAGTCTAAGGAAAACCCAAGTATTCAGCCTTCTGAATTAAGGTGTGCTCTTg                                                                                                                                                               >  2:77022/1‑251 (MQ=255)
                                                                                             ggATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCAAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATcc                                                                                                                                                       <  1:155986/250‑1 (MQ=255)
                                                                                              gATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCg                                                                                                                                                     <  1:319143/251‑1 (MQ=255)
                                                                                                  gTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTAcc                                                                                                                                                 >  1:454523/1‑251 (MQ=255)
                                                                                                   ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGACGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAACTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGACTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTAcc                                                                                                                                                 >  2:556082/1‑250 (MQ=255)
                                                                                                   ttCGTTTAAGACCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGACAATCCTTCTGACTCGTCTTTGCAGGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAACTGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTAAGGAAAACCACAGTATTCAGCCTTCTGAATTACGTTGTGCTCTTGGATCATCCTGTAcc                                                                                                                                                 >  2:125019/1‑250 (MQ=255)
                                                                                                      gTTTACGACCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGCTTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGCCGTTGTTTCTGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTATTGGATCATCCCGTACAAacg                                                                                                                                             >  2:805658/1‑251 (MQ=255)
                                                                                                           cGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAacgcgacg                                                                                                                                        >  2:560678/1‑251 (MQ=255)
                                                                                                              ccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAGCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGt                                                                                                                                     <  2:15317/251‑1 (MQ=255)
                                                                                                                aTTGAACAAATGCTAAAATATCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGCGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTAt                                                                                                                                   >  2:43443/1‑251 (MQ=255)
                                                                                                                   gAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTg                                                                                                                                 <  1:498371/250‑1 (MQ=255)
                                                                                                                       aaaTGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGAt                                                                                                                            <  1:580707/251‑1 (MQ=255)
                                                                                                                          attcaaaactttcgcccatgccgccACGAAGATTCTCCTTATCAGGAGATCGTTCTGAATCGTCGTTGCATGCACATGCAAAGCAAGCTGCTGGACAACCGCAATAAAATTCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCCTGTGTCTCAGGAAAACCACATTATCCCGCCTGCTGAATTAAGTTGTGCTCTTGGAACAGCACGTACCAACTTGACGCCGATTGCCGATGa                                                                                                                          <  2:116078/232‑1 (MQ=255)
                                                                                                                                    tccgcgCCCGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGaaaa                                                                                                                 <  1:42571/248‑1 (MQ=255)
                                                                                                                                         gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTAGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCgg                                                                                                           <  1:425730/251‑1 (MQ=255)
                                                                                                                                            aGCCGCCACGAAGATTTTCCTTATCAGGAGAGCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCTGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTg                                                                                                        <  1:125019/251‑1 (MQ=255)
                                                                                                                                                                 tATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGc                                                                                     <  1:845795/249‑1 (MQ=255)
                                                                                                                                                                    cAGGAGATCCTGCTGAATCGTCTTTGCATGCACATGCAAATCAAGCTGCTGGAGAACAGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCGCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATaa                                                                                <  2:515451/251‑1 (MQ=255)
                                                                                                                                                                      ggAGATCCTTCTGACTCGTCTTTGCCTGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACCGTATTCAGCCTTCTGAATTAAGTTGAGCTCTTGGATAATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTGGGATCGAACGTCGTGAAAGCgctaccg                                                                              >  2:459184/1‑245 (MQ=255)
                                                                                                                                                                          aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATcg                                                                          <  2:499744/251‑1 (MQ=255)
                                                                                                                                                                              ttCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATcgccgc                                                                      >  2:138864/1‑251 (MQ=255)
                                                                                                                                                                                      cGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGctgcctgc                                                               <  1:544240/250‑1 (MQ=255)
                                                                                                                                                                                          tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATg                                                                                                                           <  2:485309/186‑1 (MQ=255)
                                                                                                                                                                                          tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATg                                                                                                                           >  1:485309/1‑186 (MQ=255)
                                                                                                                                                                                            tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGc                                                         <  1:138864/250‑1 (MQ=255)
                                                                                                                                                                                                   acatgcaAAGCAAGCTGCGTGAGAACCGCAATAAAATGCTGAAGGCTTAGGGAAGTAACGAGACGTTGCTTATGACGGTGATTACGCTTGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGTATTAAGTTGTGCTCTTGGATCATCCCGTACCAACTCTATGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGTACGTCGTGAAAGCGATAACGAgcgcagcggtctgcagctgtaatgaaca                                                 <  2:529765/251‑25 (MQ=255)
                                                                                                                                                                                                    catgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTAAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACgg                                                <  2:118630/251‑1 (MQ=255)
                                                                                                                                                                                                              aagcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCa                                      <  1:713269/251‑1 (MQ=255)
                                                                                                                                                                                                               agcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCat                                     <  1:62154/251‑2 (MQ=255)
                                                                                                                                                                                                               agcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCAc                                     <  1:62122/251‑1 (MQ=255)
                                                                                                                                                                                                               agcTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTGGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGGCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCAc                                     <  1:872939/251‑1 (MQ=255)
                                                                                                                                                                                                                   gctgGAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGGTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCGTCTGAATTAAGTTGGGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAATCGATAACTATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTTACGAg                                  <  1:502671/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggAGAACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGtt                                <  1:43443/249‑1 (MQ=255)
                                                                                                                                                                                                                         gaACCGCAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTgc                           >  2:715251/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAATAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGACTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAgg                     >  1:3683/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGtt                   >  1:186672/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTAAAATGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGtt                   >  1:141512/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tAAAATGCTGAAAGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTAGAGTCTCAGGAAAACCACAGTATTCAGACTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTCACGGAAACAGGTCACGAGATTTcgcgcgcggatt                  >  2:49687/1‑245 (MQ=255)
                                                                                                                                                                                                                                   aaaaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGtttt                 >  2:816888/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTAGGGCGTTGATTACGCTGGAGTCGCTGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTTCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATGAACGGAAAAAGGTCACGAGTTTTTGCGCGCGGTTTTac               <  2:32665/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTGCTGAAGGCTCAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGACAACAACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGGAAAACGCGGTTGGATCGACCGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCCATTAACGGAAAAAGGTCACGCGTTTTTGCGCTCGGTTTTacc              >  2:569842/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  cAGGGAATTAACGAGTCGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCTATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATaa  <  1:805658/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cAGGGAATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCAt    <  1:556082/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
ACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAA  >  NC_000913/2810679‑2811170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: