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

ATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATT  >  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)
                                gTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAgg                                                                                                                                                                                                                 <  1:236775/251‑1 (MQ=255)
                                             agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                                                                                                    <  1:538321/247‑1 (MQ=255)
                                                  cAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTACCCCAAGCGCTAAATATAAGGGtgtgggatttgtcccg                                                                                                                                                                                                >  1:172517/1‑236 (MQ=255)
                                                    gAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACTCCCCTAAGCGCTAACTTaaagggtggggggttatgcctgtt                                                                                                                                                                                              >  1:232317/1‑229 (MQ=255)
                                                             gctgcaagagcggcaagaAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGAGGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                                                                                                     <  2:376911/250‑1 (MQ=255)
                                                                      gctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                                                                                                           <  1:53315/251‑1 (MQ=255)
                                                                          caagaAAGCTGCGTTGGCCCAACCGGTGAGCACACCAGCCCGCGATGTAAGAGCTTCTTTGATGGTGAGAAGTAGTTTTTTGACGGTCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTGGGGTAAAAGTAGGGTGCGTTGAAGATATGCGAGAACCTGAAAAGTGGCGGGGATCTCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                                                                                                       <  2:152633/251‑1 (MQ=37)
                                                                                    tcgctGGCCCAAACGGGGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGAACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAGGAATTACTCTCAc                                                                                                                                                              <  2:129217/246‑1 (MQ=38)
                                                                                      ggTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGCTGGTGTGACGACGTTTTTTGACGACCTGACCCGCTAAACTCCCCAAAGGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTACTCCCAAAAACGCATACATAAAGGGTGTGGCATTACGCCCCATCTGAATTTATGTGCCGATGAATTAATCTCACgct                                                                                                                                                           >  1:286245/1‑249 (MQ=255)
                                                                                               aCCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                                                                                                                                  <  1:668826/251‑1 (MQ=255)
                                                                                                ccGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                                                                                                 <  1:432761/251‑1 (MQ=255)
                                                                                                    tgAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGGGACTCCCCTAAGCGCTAACTTAAAGGGTTTGGTGTTTCGCCTGATATGATTTTACGTGGCGGTGGATTAATCTCACGATAAAAGGTTAGCAAt                                                                                                                                             >  1:214277/2‑251 (MQ=255)
                                                                                                            acaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                                                                                                                     <  1:639051/251‑1 (MQ=255)
                                                                                                             caGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCa                                                                                                                                    <  1:415108/251‑1 (MQ=255)
                                                                                                                             agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                                                                                                                    <  1:189796/251‑1 (MQ=255)
                                                                                                                              gctgcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGACTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCTgaag                                                                                                                   <  2:420768/247‑1 (MQ=255)
                                                                                                                                    ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                                                                                                                                                                                       <  2:327524/177‑1 (MQ=255)
                                                                                                                                       aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCTTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTTGCGTGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTTATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                                                                                                          <  1:517157/251‑1 (MQ=25)
                                                                                                                                       aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                                                                                                          <  1:463549/251‑1 (MQ=255)
                                                                                                                                       aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATa                                                                                                          <  2:214277/251‑1 (MQ=255)
                                                                                                                                            gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                                                                                                      <  1:392563/250‑1 (MQ=255)
                                                                                                                                                agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                                                                                                 <  1:146516/251‑1 (MQ=255)
                                                                                                                                                     ttttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCACAACCGCTAACTTTAGGGGTGGGGGATTAAGCCTGCTATGATTTAAAGTGCCGCTGAATTACTTTCACGCTAACTGGTCAGCACTTTTGGCCCCTATTTGTTAGACCCAAGAACTGGATACTTCGGCACgcaa                                                                                             >  1:63674/1‑247 (MQ=255)
                                                                                                                                                           gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                                                                                                    <  1:108537/141‑1 (MQ=255)
                                                                                                                                                                   gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGAGATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                                                                                                                                                                                              <  2:556284/123‑1 (MQ=255)
                                                                                                                                                                     cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTAAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAAccc                                                                            <  1:415268/251‑1 (MQ=255)
                                                                                                                                                                                     aaaGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCg                                                           <  1:346947/251‑1 (MQ=255)
                                                                                                                                                                                                gTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAGGAATTACTCTCACGATAACTGGTCa                                                                                                                                                  <  2:647227/154‑1 (MQ=21)
                                                                                                                                                                                                     cAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTAAACGCTAACTTAAAGGGTGTGGGATTAAGCCTGGTATGGTTTTACGTGGCGCTGGATTAATCTCTCGATAAATTGTCAGCAATTCTTGCCCATATTTGTAAACCCGAAGAACAGGATAATTCGGCACGCAATGCCGGGGCTCTAACCCCCCCCCGCGGAAATTGGGATGCTGCAAatctc                                           >  1:162205/1‑250 (MQ=255)
                                                                                                                                                                                                       ggcttGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCt                                         <  2:232317/247‑1 (MQ=255)
                                                                                                                                                                                                         cttGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATCCCCTTAAAGCTAACTTAAGGGGTGTGGTATTTAGCCCGGTATTATTTTAAATGCCGGTGGAATAATCTCACGATAACTGGTCAGCAATTCTGGCCCATATTTGTAAGCCCCAAGGAATGGATACTTTGGCACGTAAAGCCGGGGCACTAACCCCCCCCCGCGAAATTTGTGATGCTGCGACACTGCttc                                       >  1:665730/3‑249 (MQ=9)
                                                                                                                                                                                                                aaaatTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTgg                                 <  1:481169/250‑1 (MQ=255)
                                                                                                                                                                                                                   attAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGggct                             <  1:93929/249‑1 (MQ=255)
                                                                                                                                                                                                                   attAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTGCGGCACGTAATGCCGGGGCTCTAGCCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGggct                             <  1:40183/249‑1 (MQ=255)
                                                                                                                                                                                                                                aaaGATATGCGAGCACGTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAAGCCTGATATGAGTTAACTGGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTGCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTTCGTCTGGTGCTGGCTTACGGcccc                <  2:63674/250‑1 (MQ=14)
                                                                                                                                                                                                                                    atgtGCGAGCACCTGAAAAGTGGCGTGGATCACTCCCAAAAGCGCTATCTTAAGGGTTTTGGTATTGCGGCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAATCCCGAAGAACTGGATACTTGGGCACGTAAGGCCGGTGCTCTAACCCGCCGCCGCGAAAGTCGTGATGCGGCAACTCTGCTACGTCTGGGGCTTGCTTACGGCCCCgggg            <  2:310970/248‑1 (MQ=14)
                                                                                                                                                                                                                                      atGCGAGCACCTGAAAAGTGGCGGGGTTCACTCCCATAAGCGCTAACTTAAGGGTGGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTATCCCGCCGCCGCGAAATTCGTGATGATGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCggggg           <  2:140686/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCACCTGAAAAGTGGCGGGGATCACTCCCATCAACGCTCACTTAAGGGGTGTGGGATTTCGCCTGATATGATTTTAAGGGCCGGTGAATTTATCTCACGATAAATGGGCCGCAATTCTGGCCCCTTTTTGTTAACCCCCAGAAATGGATACTTCGGCACGTAATGCCGGGGCTCTAAACCCCCCCCCCCAAATTCGTGGTGCTGCAACTCTGCTACCGCTGGGGCTGGCTTACGGCCCCcggggggtggg     >  1:105210/1‑246 (MQ=9)
                                                                                                                                                                                                                                            gCACCTGAAAAGTGGCGGGGAGCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATACCTGGTCAGCAATTCTGGCACATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCAGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACTGCCCCGGGGGGATGTCa    <  2:276008/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCAtt  <  1:492700/251‑1 (MQ=37)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
ATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATT  >  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: