Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 939,377 T→C 8.2% intergenic (+40/‑51) rarA → / → serS recombination intermediate processing DNA‑dependent ATPase/seryl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913939,3770TC 8.2% 150.2 / 2.2 55intergenic (+40/‑51)rarA/serSrecombination intermediate processing DNA‑dependent ATPase/seryl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base C (2/3);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTGAACAAACTGGGCG  >  NC_000913/939079‑939665
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
aagCCCGGATTATGTCGTTCAGGTTCATTTGCGTAATGTGCCAACTAAAATAATGAAGGAAATGGGCTACGGGCAGGAATATCGATACGCTCATGATGAAGCAAACGGTTATGCTGCCGGTGATGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCg                                                                                                                                                                                                                                                                                                <  2:11318/299‑1 (MQ=255)
 cGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGTTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGAGATCgg                                                                                                                                                                                                                                                                                               >  2:242255/1‑300 (MQ=255)
                gggCCGGTTCATTTTCTTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAGGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGCCACAGCCTGCGACcacag                                                                                                                                                                                                                                                                                <  1:246132/299‑2 (MQ=255)
                   ccGGTTCATTTGCGTAATGCGCCGTCGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGCCACAGCCTGCGACcacaggga                                                                                                                                                                                                                                                                             <  1:609017/301‑5 (MQ=255)
                    cGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGACGGCTGAACAGGATCAAAATAGCACCATAAAACGCTACCGTTAATGTTATCGTTGAGGTAATGTTGTTACTGTAT‑CCCTTTGTCCGCAAACTGGGGCcacaccacc                                                                                                                                                                                                                                                                            >  1:621380/1‑296 (MQ=255)
                     ggttgtttgtgggagttcgtcgaagaATTTAAGGAAGGAAGTGGGTTACGGGAAGGAAGATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTGCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTccc                                                                                                                                                                                                                                                                           <  1:377456/281‑1 (MQ=255)
                         cATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAAAGCAACCGTTAAGGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCGGTGGCCACATAGCCCAttt                                                                                                                                                                                                                                                                       >  2:214795/1‑301 (MQ=255)
                          ttttGGGGAATGCGCCGACGAAATTAATTAAGGAAATGGGCTACGGGCAGGAATATGGTTACGCTCATGATTAAGCAAACGCTTATGCTGCCGGTGGGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGCCACAGCCTGCGACcacagggatacagta                                                                                                                                                                                                                                                                      <  1:376671/300‑12 (MQ=255)
                          aTTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGCCACAGCCTGCGACcaaagggatacagta                                                                                                                                                                                                                                                                      >  2:377167/1‑288 (MQ=255)
                                      cgcCGACGAAATTAATGAAGGAAATGGGTTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGca                                                                                                                                                                                                                                                          <  2:622239/301‑1 (MQ=255)
                                         tgctcttccgatctaaGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCGCAg                                                                                                                                                                                                                                                       <  2:300683/287‑1 (MQ=255)
                                            cgaAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGAt                                                                                                                                                                                                                                                    <  2:465792/301‑1 (MQ=255)
                                             gaAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATa                                                                                                                                                                                                                                                   >  2:594363/1‑301 (MQ=255)
                                                     tGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCCGTGAGGTTTACTTCCCGCCGGAAATAGCCCAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAACACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGAAGGCTGTGGCAACACCTCCCATTTAATTCGATAAGCCCAGGATAAgactcct                                                                                                                                                                                                                                           >  1:329994/1‑295 (MQ=255)
                                                      gAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCAAAGGATAAGCATGCTc                                                                                                                                                                                                                                          >  2:131187/1‑301 (MQ=255)
                                                       aaGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTCATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCCCCCATTTAATTCGATAAGCGCAgaaatcggaagagcg                                                                                                                                                                                                                                         >  1:300279/1‑287 (MQ=255)
                                                               gggTTACGGGCAGGAATGTCGTTACGCTCATGATGAAGCAACCTCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACTGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAAt                                                                                                                                                                                                                                 <  1:96260/301‑1 (MQ=255)
                                                                ggCTACGGTCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCCGTGAGGTTTACTTCCCGCCGGAAATAGCCCAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATc                                                                                                                                                                                                                                <  2:330399/301‑1 (MQ=255)
                                                                     cGGGCATGAATATCGTTACGCTCTTGATGAAGCAAACGCTTATGCTGCCTGTTAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAACCAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATctgctg                                                                                                                                                                                                                           <  1:214552/301‑1 (MQ=255)
                                                                         aaGGAATATCGTTACGCTAATGATGAAGCAAACTCTTAGGCGTCCGGTGAGGTTTACTTCCCTCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGATGTTACTGTAT‑CCGTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTa                                                                                                                                                                                                                       <  1:328087/300‑1 (MQ=255)
                                                                               tatCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGc                                                                                                                                                                                                                 >  2:378260/1‑301 (MQ=255)
                                                                                          tCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCGACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTCATTCGATCAGCCAAGGAAAAGCATGCTCGATCCAAATCTGCTGCGTAATGAGACAAACGctagc                                                                                                                                                                                                      >  1:623101/1‑297 (MQ=255)
                                                                                              gatgaAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGACCGATACCAATCTGCTGCGTAATGAGCCAGAAGCAGTccccg                                                                                                                                                                                                  >  1:267787/1‑297 (MQ=255)
                                                                                                     cAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAACAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATACCAATCTGCTGCGTAATGAGCCAGACGCAGACGCTGAAAAACt                                                                                                                                                                                           >  1:439762/1‑301 (MQ=255)
                                                                                                      aaaCGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTAAATGTTATCGTTGCGGTAATGTTGTTACCGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGAGTAATGAGCCAGAAGCAATCGCTGAAAAAccg                                                                                                                                                                                          >  1:263496/1‑299 (MQ=255)
                                                                                                                  tgcCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGc                                                                                                                                                                              <  2:559318/301‑1 (MQ=255)
                                                                                                                             tttACTTCCCGACGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTGCCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGa                                                                                                                                                                   <  2:443789/301‑1 (MQ=255)
                                                                                                                                tcTTCCCGCCGGAAATAGCACAAACACGTTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGt                                                                                                                                                                <  2:440326/300‑1 (MQ=255)
                                                                                                                                 cctttccctacaggacgctcttccgatcTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTa                                                                                                                                                               <  1:533662/274‑1 (MQ=255)
                                                                                                                                  ttCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTACTGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATAAGCCAGACGCAGTCGACGAAAAACTGGCACGCCGGGGCTTTAACAGGGATGTAg                                                                                                                                                              >  1:548798/1‑301 (MQ=255)
                                                                                                                                        ccgGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGc                                                                                                                                                        <  2:561767/301‑1 (MQ=255)
                                                                                                                                                  aacaAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGctct                                                                                                                                              <  1:278590/300‑1 (MQ=255)
                                                                                                                                                              attattTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTaagatcggaagagcacacgtctgaactcca                                                                                                                                  >  2:534335/1‑272 (MQ=255)
                                                                                                                                                                    tCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCAAGCAGGGACTTTAAGATGGAAGTAGAAAAACAGGGCGATCTTGAAGAACGTCGAAAAGt                                                                                                                            >  1:360664/1‑301 (MQ=255)
                                                                                                                                                                          aaaaCAGGTGCCTTGAAGGTAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAACACGCTACCGTTAAGGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCa                                                                                                                      <  1:593602/300‑1 (MQ=255)
                                                                                                                                                                           aaaCAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCACTTGAAGAGCGTCGTAAAGTAATTCAg                                                                                                                     >  1:215183/1‑301 (MQ=255)
                                                                                                                                                                                          aGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATTGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAAc                                                                                                      >  1:42380/1‑301 (MQ=255)
                                                                                                                                                                                              aaGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGc                                                                                                  <  2:215426/301‑1 (MQ=255)
                                                                                                                                                                                               aGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATTGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCa                                                                                                 <  2:42433/301‑1 (MQ=255)
                                                                                                                                                                                                      cGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAAACTGCAAGCGGaa                                                                                          >  2:412686/1‑300 (MQ=255)
                                                                                                                                                                                                        aaaaGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATATAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGttttccgttttgacctgcaatactttac                                                                                        <  2:623960/301‑28 (MQ=255)
                                                                                                                                                                                                                       tggctgAACAGGTTCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGaa                                                                         <  2:42847/301‑1 (MQ=255)
                                                                                                                                                                                                                         gctgAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGaaaa                                                                       >  2:361067/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     aaCGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTAGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGa                                           <  2:439741/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggttccGTTAATGTTTTCGTTGCGGTAATGTTGTT‑CTTTGTGTCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGa                                        <  1:131062/297‑1 (MQ=255)
                                                                                                                                                                                                                                                          tACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCAGCAAGCGGAGCGGAACTCCCGAACGAAAACCATTGGCCAGGCGAAAGCGCGCGGGGaacaat                                      >  1:97373/1‑297 (MQ=255)
                                                                                                                                                                                                                                                             cGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCg                                   >  2:578246/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 aaTGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGAACGAAATCCACAGGCCAGGCGAAAGCGCGCGGGGAAGATATccacac                               >  1:78190/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                     ttATCGTTGCGGTAATGTTGTTATTGTAT‑CCCTTTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTa                           <  2:268055/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ttATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTa                           <  2:97497/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ggTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAg                 <  2:623573/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ttaTGTTGTTACTGTTT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTg               <  1:305514/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aTGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCAAGCCTTTACGTCTGGAAGTGaa             >  2:305918/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          tACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGACGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGAGTCGTAAAGTATTGCAGGTCAAAACGGAAACACTGCAAGCGGAGCGTAACTCCCGATCGAAATACATTGGCCAGGCGAAAACGCGCGGGGAAGATCACGACCCTTACAGGATGGAAGTGAACAAAATg      >  1:429589/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               gttggccttgggTGGCAGGGTGTGGCCAAATGTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTGAACAAACTGGGCg  <  1:377764/291‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
ACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTAT‑CCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTGAACAAACTGGGCG  >  NC_000913/939079‑939665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑