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

ACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTGAACAAACTGGGCG  >  NC_000913/939079‑939665
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
aagCCCGGATTATGTCGTTCAGGTTCATTTGCGTAATGTGCCAACTAAAATAATGAAGGAAATGGGCTACGGGCAGGAATATCGATACGCTCATGATGAAGCAAACGGTTATGCTGCCGGTGATGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCg                                                                                                                                                                                                                                                                                                <  2:11318/299‑1 (MQ=255)
 cGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGTTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGAGATCgg                                                                                                                                                                                                                                                                                               >  2:242255/1‑300 (MQ=255)
                gggCCGGTTCATTTTCTTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTAGGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCCTGCGACcacag                                                                                                                                                                                                                                                                                <  1:246132/299‑2 (MQ=255)
                   ccGGTTCATTTGCGTAATGCGCCGTCGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCCTGCGACcacaggga                                                                                                                                                                                                                                                                             <  1:609017/301‑5 (MQ=255)
                    cGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGACGGCTGAACAGGATCAAAATAGCACCATAAAACGCTACCGTTAATGTTATCGTTGAGGTAATGTTGTTACTGTATCCCTTTGTCCGCAAACTGGGGCcacaccacc                                                                                                                                                                                                                                                                            >  1:621380/1‑296 (MQ=255)
                     ggttgtttgtgggagttcgtcgaagaATTTAAGGAAGGAAGTGGGTTACGGGAAGGAAGATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTGCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTccc                                                                                                                                                                                                                                                                           <  1:377456/281‑1 (MQ=255)
                         cATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAAAGCAACCGTTAAGGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCGGTGGCCACATAGCCCAttt                                                                                                                                                                                                                                                                       >  2:214795/1‑301 (MQ=255)
                          ttttGGGGAATGCGCCGACGAAATTAATTAAGGAAATGGGCTACGGGCAGGAATATGGTTACGCTCATGATTAAGCAAACGCTTATGCTGCCGGTGGGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCCTGCGACcacagggatacagta                                                                                                                                                                                                                                                                      <  1:376671/300‑12 (MQ=255)
                          aTTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCCTGCGACcaaagggatacagta                                                                                                                                                                                                                                                                      >  2:377167/1‑288 (MQ=255)
                                      cgcCGACGAAATTAATGAAGGAAATGGGTTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGca                                                                                                                                                                                                                                                          <  2:622239/301‑1 (MQ=255)
                                         tgctcttccgatctaaGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCGCAg                                                                                                                                                                                                                                                       <  2:300683/287‑1 (MQ=255)
                                            cgaAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGAt                                                                                                                                                                                                                                                    <  2:465792/301‑1 (MQ=255)
                                             gaAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATa                                                                                                                                                                                                                                                   >  2:594363/1‑301 (MQ=255)
                                                     tGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCCGTGAGGTTTACTTCCCGCCGGAAATAGCCCAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAACACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGAAGGCTGTGGCAACACCTCCCATTTAATTCGATAAGCCCAGGATAAgactcct                                                                                                                                                                                                                                           >  1:329994/1‑295 (MQ=255)
                                                      gAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCAAAGGATAAGCATGCTc                                                                                                                                                                                                                                          >  2:131187/1‑301 (MQ=255)
                                                       aaGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTCATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCCCCCATTTAATTCGATAAGCGCAgaaatcggaagagcg                                                                                                                                                                                                                                         >  1:300279/1‑287 (MQ=255)
                                                               gggTTACGGGCAGGAATGTCGTTACGCTCATGATGAAGCAACCTCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACTGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAAt                                                                                                                                                                                                                                 <  1:96260/301‑1 (MQ=255)
                                                                ggCTACGGTCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCCGTGAGGTTTACTTCCCGCCGGAAATAGCCCAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATc                                                                                                                                                                                                                                <  2:330399/301‑1 (MQ=255)
                                                                     cGGGCATGAATATCGTTACGCTCTTGATGAAGCAAACGCTTATGCTGCCTGTTAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAACCAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATctgctg                                                                                                                                                                                                                           <  1:214552/301‑1 (MQ=255)
                                                                         aaGGAATATCGTTACGCTAATGATGAAGCAAACTCTTAGGCGTCCGGTGAGGTTTACTTCCCTCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGATGTTACTGTATCCGTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTa                                                                                                                                                                                                                       <  1:328087/300‑1 (MQ=255)
                                                                               tatCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGc                                                                                                                                                                                                                 >  2:378260/1‑301 (MQ=255)
                                                                                          tCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCGACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTCATTCGATCAGCCAAGGAAAAGCATGCTCGATCCAAATCTGCTGCGTAATGAGACAAACGctagc                                                                                                                                                                                                      >  1:623101/1‑297 (MQ=255)
                                                                                              gatgaAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGACCGATACCAATCTGCTGCGTAATGAGCCAGAAGCAGTccccg                                                                                                                                                                                                  >  1:267787/1‑297 (MQ=255)
                                                                                                     cAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAACAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATACCAATCTGCTGCGTAATGAGCCAGACGCAGACGCTGAAAAACt                                                                                                                                                                                           >  1:439762/1‑301 (MQ=255)
                                                                                                      aaaCGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTAAATGTTATCGTTGCGGTAATGTTGTTACCGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGAGTAATGAGCCAGAAGCAATCGCTGAAAAAccg                                                                                                                                                                                          >  1:263496/1‑299 (MQ=255)
                                                                                                                  tgcCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGc                                                                                                                                                                              <  2:559318/301‑1 (MQ=255)
                                                                                                                             tttACTTCCCGACGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTGCCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGa                                                                                                                                                                   <  2:443789/301‑1 (MQ=255)
                                                                                                                                tcTTCCCGCCGGAAATAGCACAAACACGTTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGt                                                                                                                                                                <  2:440326/300‑1 (MQ=255)
                                                                                                                                 cctttccctacaggacgctcttccgatcTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTa                                                                                                                                                               <  1:533662/274‑1 (MQ=255)
                                                                                                                                  ttCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTACTGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATAAGCCAGACGCAGTCGACGAAAAACTGGCACGCCGGGGCTTTAACAGGGATGTAg                                                                                                                                                              >  1:548798/1‑301 (MQ=255)
                                                                                                                                        ccgGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGc                                                                                                                                                        <  2:561767/301‑1 (MQ=255)
                                                                                                                                                  aacaAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGctct                                                                                                                                              <  1:278590/300‑1 (MQ=255)
                                                                                                                                                              attattTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTaagatcggaagagcacacgtctgaactcca                                                                                                                                  >  2:534335/1‑272 (MQ=255)
                                                                                                                                                                    tCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCAAGCAGGGACTTTAAGATGGAAGTAGAAAAACAGGGCGATCTTGAAGAACGTCGAAAAGt                                                                                                                            >  1:360664/1‑301 (MQ=255)
                                                                                                                                                                          aaaaCAGGTGCCTTGAAGGTAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAACACGCTACCGTTAAGGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCa                                                                                                                      <  1:593602/300‑1 (MQ=255)
                                                                                                                                                                           aaaCAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCACTTGAAGAGCGTCGTAAAGTAATTCAg                                                                                                                     >  1:215183/1‑301 (MQ=255)
                                                                                                                                                                                          aGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATTGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAAc                                                                                                      >  1:42380/1‑301 (MQ=255)
                                                                                                                                                                                              aaGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGc                                                                                                  <  2:215426/301‑1 (MQ=255)
                                                                                                                                                                                               aGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATTGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCa                                                                                                 <  2:42433/301‑1 (MQ=255)
                                                                                                                                                                                                      cGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAAACTGCAAGCGGaa                                                                                          >  2:412686/1‑300 (MQ=255)
                                                                                                                                                                                                        aaaaGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATATAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGttttccgttttgacctgcaatactttac                                                                                        <  2:623960/301‑28 (MQ=255)
                                                                                                                                                                                                                       tggctgAACAGGTTCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGaa                                                                         <  2:42847/301‑1 (MQ=255)
                                                                                                                                                                                                                         gctgAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGaaaa                                                                       >  2:361067/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     aaCGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTAGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGa                                           <  2:439741/301‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggttccGTTAATGTTTTCGTTGCGGTAATGTTGTTCTTTGTGTCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGa                                        <  1:131062/297‑1 (MQ=255)
                                                                                                                                                                                                                                                          tACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCAGCAAGCGGAGCGGAACTCCCGAACGAAAACCATTGGCCAGGCGAAAGCGCGCGGGGaacaat                                      >  1:97373/1‑297 (MQ=255)
                                                                                                                                                                                                                                                             cGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCg                                   >  2:578246/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 aaTGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGAACGAAATCCACAGGCCAGGCGAAAGCGCGCGGGGAAGATATccacac                               >  1:78190/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                     ttATCGTTGCGGTAATGTTGTTATTGTATCCCTTTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTa                           <  2:268055/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ttATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTa                           <  2:97497/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ggTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAg                 <  2:623573/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 ttaTGTTGTTACTGTTTCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTg               <  1:305514/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aTGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCAAGCCTTTACGTCTGGAAGTGaa             >  2:305918/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          tACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGACGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGAGTCGTAAAGTATTGCAGGTCAAAACGGAAACACTGCAAGCGGAGCGTAACTCCCGATCGAAATACATTGGCCAGGCGAAAACGCGCGGGGAAGATCACGACCCTTACAGGATGGAAGTGAACAAAATg      >  1:429589/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               gttggccttgggTGGCAGGGTGTGGCCAAATGTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTGAACAAACTGGGCg  <  1:377764/291‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
ACGCCCGGATTATGACGTGCCGGTTCATTTGCGTAATGCGCCGACGAAATTAATGAAGGAAATGGGCTACGGGCAGGAATATCGTTACGCTCATGATGAAGCAAACGCTTATGCTGCCGGTGAGGTTTACTTCCCGCCGGAAATAGCACAAACACGCTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGAAGAGCGTCGTAAAGTATTGCAGGTCAAAACGGAAAACCTGCAAGCGGAGCGTAACTCCCGATCGAAATCCATTGGCCAGGCGAAAGCGCGCGGGGAAGATATCGAGCCTTTACGTCTGGAAGTGAACAAACTGGGCG  >  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: