Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 15,369 T→A intergenic (+71/‑76) dnaJ → / → insL1 chaperone protein DnaJ/IS186/IS421 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,3690TA96.1% 148.1 / ‑1.8 52intergenic (+71/‑76)dnaJ/insL1chaperone protein DnaJ/IS186/IS421 family transposase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (10/39);  total (12/40)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATT  >  NC_000913/15126‑15614
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
aTTTGCTGGGCCGTGTTGTCGTCGAAACTCCGGTCGGCCTGAACGAGAATCAGAAACAGCTGCTGCAAGAGCTGCAAGAAATCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTGTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCATGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTgg                                                                                                                                                                                                                                                   <  2:429318/251‑1 (MQ=255)
   tGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCg                                                                                                                                                                                                                                                 <  1:582369/250‑1 (MQ=255)
     cTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGCAAAGTGGCgggg                                                                                                                                                                                                                                              >  1:516477/1‑251 (MQ=255)
          ccGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGagat                                                                                                                                                                                                                                          >  1:59642/1‑247 (MQ=255)
             ggTTGCCGTCGAAACACCGGTAGGCCTGAATGAGAAGCAGAAACAGCTGTTGCAAGATCTGCAAGAAAGCTTCGGTGGACCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTGTGACGACCTGACTCGCTAACCTCCCCAAAAGCCTGCCCGTTGGCAGGATTGGGTAAAATTAGGGTTCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTcc                                                                                                                                                                                                                                      <  1:455746/250‑1 (MQ=255)
                 gtcgtcGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                                                                                                                                                                                                                   <  1:345227/250‑1 (MQ=255)
                   cgtcgAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAg                                                                                                                                                                                                                                <  1:218810/251‑1 (MQ=255)
                       gAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGaacaaacccctcagacgc                                                                                                                                                                                                                            >  1:229314/1‑234 (MQ=255)
                         aaccccGGTAGGCCTGTACAAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCATCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCGGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAATTGGCGGGGATCACTCTCATAAGCGCTaa                                                                                                                                                                                                                          <  2:17864/247‑1 (MQ=255)
                         aaCACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAcaacccacaaggaca                                                                                                                                                                                                                           >  2:574735/1‑236 (MQ=255)
                          acacCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc                                                                                                                                                                                                                         <  1:451476/251‑1 (MQ=255)
                           cacCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                                                                                                                                                                                        <  1:469343/251‑1 (MQ=255)
                                gTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑gg                                                                                                                                                                                                                 <  1:236775/251‑1 (MQ=255)
                                             agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTAc                                                                                                                                                                                                    <  1:538321/247‑1 (MQ=255)
                                                  cAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTACCCCAAGCGCTAAATATAA‑GGGtgtgggatttgtcccg                                                                                                                                                                                                >  1:172517/1‑236 (MQ=255)
                                                    gAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACTCCCCTAAGCGCTAACT‑Taaagggtggggggttatgcctgtt                                                                                                                                                                                              >  1:232317/1‑229 (MQ=255)
                                                             gctgcaagagcggcaagaAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGAGGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                                                                                                     <  2:376911/250‑1 (MQ=255)
                                                                      gctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                                                                                                           <  1:53315/251‑1 (MQ=255)
                                                                          caagaAAGCTGCGTTGGCCCAACCGGTGAGCACACCAGCCCGCGATGTAAGAGCTTCTTTGATGGTGAGAAGTAGTTTTTTGACGGTCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTGGGGTAAAAGTAGGGTGCGTTGAAGATATGCGAGAACCTGAAAAGTGGCGGGGATCTCTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                                                                                                       <  2:152633/251‑1 (MQ=37)
                                                                                    tcgctGGCCCAAACGGGGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGAACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAGGAATTACTCTCAc                                                                                                                                                              <  2:129217/246‑1 (MQ=38)
                                                                                      ggTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGCTGGTGTGACGACGTTTTTTGACGACCTGACCCGCTAAACTCCCCAAAGGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTACTCCCAAAAACGCATACA‑TAA‑AGGGTGTGGCATTACGCCCCATCTGAATTTATGTGCCGATGAATTAATCTCACgct                                                                                                                                                           >  1:286245/1‑249 (MQ=255)
                                                                                               aCCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                                                                                                                                  <  1:668826/251‑1 (MQ=255)
                                                                                                ccGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                                                                                                 <  1:432761/251‑1 (MQ=255)
                                                                                                    tgAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGGGACTCCCCTAAGCGCTAACT‑TAAAGGGTT‑TGGTGTTTCGCCTGATATGATTTTACGTGGCGGTGGATTAATCTCACGATAAAAGGTTAGCAAt                                                                                                                                             >  1:214277/2‑251 (MQ=255)
                                                                                                            acaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                                                                                                                     <  1:639051/251‑1 (MQ=255)
                                                                                                             caGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCa                                                                                                                                    <  1:415108/251‑1 (MQ=255)
                                                                                                                             agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                                                                                                                    <  1:189796/251‑1 (MQ=255)
                                                                                                                              gctgcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGACTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCTgaag                                                                                                                   <  2:420768/247‑1 (MQ=255)
                                                                                                                                    ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGAtt                                                                                                                                                                                       <  2:327524/177‑1 (MQ=255)
                                                                                                                                       aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCTTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTTGCGTGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTTATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                                                                                                          <  1:517157/251‑1 (MQ=25)
                                                                                                                                       aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                                                                                                          <  1:463549/251‑1 (MQ=255)
                                                                                                                                       aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                                                                                                          <  2:214277/251‑1 (MQ=255)
                                                                                                                                            gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                                                                                                      <  1:392563/250‑1 (MQ=255)
                                                                                                                                                agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                                                                                                 <  1:146516/251‑1 (MQ=255)
                                                                                                                                                     ttttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCACAACCGCTAACT‑TTA‑GGGGTGGGGGATTAAGCCTGCTATGATTTAAAGTGCCGCTGAATTACTTTCACGCTAACTGGTCAGCACTTTTGGCCCCTATTTGTTAGACCCAAGAACTGGATACTTCGGCACgcaa                                                                                             >  1:63674/1‑247 (MQ=255)
                                                                                                                                                           gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTAc                                                                                                                                                                                                    <  1:108537/141‑1 (MQ=255)
                                                                                                                                                                   gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGAGATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGtt                                                                                                                                                                                                              <  2:556284/123‑1 (MQ=255)
                                                                                                                                                                     cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTAAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAAccc                                                                            <  1:415268/251‑1 (MQ=255)
                                                                                                                                                                                     aaaGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCg                                                           <  1:346947/251‑1 (MQ=255)
                                                                                                                                                                                                gTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAGGAATTACTCTCACGATAACTGGTCa                                                                                                                                                  <  2:647227/154‑1 (MQ=21)
                                                                                                                                                                                                     cAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTAAACGCTAACT‑TAA‑AGGGTGTGGGATTAAGCCTGGTATGGTTTTACGTGGCGCTGGATTAATCTCTCGATAAATTGTCAGCAATTCTTGCCCATATTTGTAAACCCGAAGAACAGGATAATTCGGCACGCAATGCCGGGGCTCTAACCCCCCCCCGCGGAAATT‑GGGATGCTGCAAatctc                                           >  1:162205/1‑250 (MQ=255)
                                                                                                                                                                                                       ggcttGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCt                                         <  2:232317/247‑1 (MQ=255)
                                                                                                                                                                                                         cttGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATCCCCTTAAAGCTAACT‑TAA‑GGGGTGTGGTATTTAGCCCGGTATTATTTTAAATGCCGGTGGAATAATCTCACGATAACTGGTCAGCAATTCTGGCCCATATTTGTAAGCCCCAAGGAATGGATACTTTGGCACGTAAAGCCGGGGCACTAACCCCCCCCCGCG‑AAATTTGTGATGCTGCGACACTGCttc                                       >  1:665730/3‑249 (MQ=9)
                                                                                                                                                                                                                aaaatTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTgg                                 <  1:481169/250‑1 (MQ=255)
                                                                                                                                                                                                                   attAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTGGggct                             <  1:93929/249‑1 (MQ=255)
                                                                                                                                                                                                                   attAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTGCGGCACGTAATGCCGGGGCTCTAGCCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTGGggct                             <  1:40183/249‑1 (MQ=255)
                                                                                                                                                                                                                                aaaGATATGCGAGCACGTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTAAGCCTGATATGAGTTAACTGGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTGCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTTCGTCTGGTGCTGGCTTACGGcccc                <  2:63674/250‑1 (MQ=14)
                                                                                                                                                                                                                                    atgtGCGAGCACCTGAAAAGTGGCGTGGATCACTCCCAAAAGCGCTATCT‑TAA‑GGGTTTTGGTATTGCGGCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAATCCCGAAGAACTGGATACTTGGGCACGTAAGGCCGGTGCTCTAACCCGCCGCCGCG‑AAAGTCGTGATGCGGCAACTCTGCTACGTCTGGGGCTTGCTTACGGCCCCgggg            <  2:310970/248‑1 (MQ=14)
                                                                                                                                                                                                                                      atGCGAGCACCTGAAAAGTGGCGGGGTTCACTCCCATAAGCGCTAACT‑TAA‑GGGTGGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTATCCCGCCGCCGCG‑AAATTCGTGATGATGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCggggg           <  2:140686/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCACCTGAAAAGTGGCGGGGATCACTCCCATCAACGCTCACT‑TAA‑GGGGTGTGGGATTTCGCCTGATATGATTTTAAGGGCCGGTGAATTTATCTCACGATAAATGGGCCGCAATTCTGGCCCCTTTTTGTTAACCCCCAGAAATGGATACTTCGGCACGTAATGCCGGGGCTCTAAACCCCCCCCCCC‑AAATTCGTGGTGCTGCAACTCTGCTACCGCTGGGGCTGGCTTACGGCCCCcggggggtggg     >  1:105210/1‑246 (MQ=9)
                                                                                                                                                                                                                                            gCACCTGAAAAGTGGCGGGGAGCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATACCTGGTCAGCAATTCTGGCACATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCAGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACTGCCCCGGGGGGATGTCa    <  2:276008/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCAtt  <  1:492700/251‑1 (MQ=37)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
ATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACT‑TAA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCG‑AAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATT  >  NC_000913/15126‑15614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑