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,3690TA98.3% 181.9 / NA 59intergenic (+71/‑76)dnaJ/insL1chaperone protein DnaJ/IS186/IS421 family transposase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (15/42);  minor base C (1/0);  total (16/43)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTA  >  NC_000913/15130‑15615
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         
gCTGTGCCGCGTGGTCGTCGAAACCCCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATGAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCggg                                                                                                                                                                                                                                                <  2:456894/251‑1 (MQ=255)
     gCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAGAGTGGCGGGGATc                                                                                                                                                                                                                                            <  1:236502/250‑1 (MQ=255)
      ccGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCTCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCgggggtacc                                                                                                                                                                                                                                          >  2:112392/1‑246 (MQ=255)
        gcgTTGTCGTCGAAACACCGGTAGGCCGGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTc                                                                                                                                                                                                                                        <  2:28400/251‑1 (MQ=255)
           ttGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATccc                                                                                                                                                                                                                                      >  1:364168/1‑250 (MQ=255)
             gtcgtcGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                                                                                                                                                                                                                    <  2:237356/250‑1 (MQ=255)
                  cgAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACcccccttagagc                                                                                                                                                                                                                              >  2:428753/1‑243 (MQ=255)
                    aaaCACCGGTAGGCCTGAACGAGCAGCAGAAACAGCTGCTGCAAAAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCggaggtcaatccccaaacagcc                                                                                                                                                                                                                             >  2:212975/1‑230 (MQ=255)
                                 ccTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTAAACGGCAAACTT‑AA‑GGGttt                                                                                                                                                                                                             >  1:349855/1‑250 (MQ=255)
                                 ccTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTtgt                                                                                                                                                                                                            <  1:175811/251‑1 (MQ=255)
                                  cTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCAATCCCCAAAACG‑CTAACTTTAA‑GGGtgtg                                                                                                                                                                                                            >  1:450319/1‑248 (MQ=255)
                                    gAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCGCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGt                                                                                                                                                                                                         <  2:47701/251‑1 (MQ=255)
                                       cgagaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTAtt                                                                                                                                                                                                      <  1:152779/245‑1 (MQ=255)
                                       cgagaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCAAAAAGCCTGCCCGTGGGCAGGCTTGGGTAACATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCTCTCCCATCAACT‑CTAACAT‑AA‑GGGGTGTgggttt                                                                                                                                                                                                      >  2:498974/7‑247 (MQ=255)
                                         agaagCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTAc                                                                                                                                                                                                    <  1:159834/247‑1 (MQ=255)
                                          gaagCAGAAACATCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTGGACGACCTGACCCGCTAACCGCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATAGGCGAGCACCTGAAAAGTGGCGGGGATCACGCCCATAAGCG‑CTAACTG‑AA‑GGGTAGTGGTATTACg                                                                                                                                                                                                   <  2:422747/248‑1 (MQ=255)
                                            agCAGAAACAGCTGCTGCCAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGACGATATGCGAGCACCTGAAAAGTGGCGGGGGTACCACCCCTAAGCG‑CTTACTT‑TA‑GGGTTGGGGTATTAAGcc                                                                                                                                                                                                 >  1:189465/2‑251 (MQ=255)
                                               aGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGa                                                                                                                                                                                              <  1:10040/251‑1 (MQ=255)
                                                 aaaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACACCCCTAAACC‑CCAACAT‑AA‑AGGTTGTggggttacaccccgt                                                                                                                                                                                             >  1:380401/1‑237 (MQ=255)
                                                  aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGatat                                                                                                                                                                                           <  2:198576/251‑1 (MQ=255)
                                                  aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGACACTCCCACAAACG‑CAAACTT‑AC‑AGGTTGGGGTAtttatccatgta                                                                                                                                                                                            >  1:454305/1‑240 (MQ=255)
                                                  aaCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGatat                                                                                                                                                                                           <  1:120610/251‑1 (MQ=255)
                                                      gctgctGCAAGAGCTGCAAGAAAGCTGCGGTTGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGATCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCGTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GTGTTGTGTTATTACGCCTGATATGAtt                                                                                                                                                                                       <  2:555209/251‑1 (MQ=255)
                                                         gctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCTCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTGGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                                                      <  1:586745/249‑1 (MQ=255)
                                                         gctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                                                      <  2:13169/249‑1 (MQ=255)
                                                          ctgcaagagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGGGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTGACACCCCTAAACG‑CTAACTT‑AA‑GGGTTGTGGTTTTTCGACCCATATGATTTatt                                                                                                                                                                                   >  1:125552/1‑249 (MQ=255)
                                                             caagagctgcaagaAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCGTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCGGACCCGCTAACCTCCCCAAAAGCCGGCCCGTGGGCAGGCTTGGGTAAAAGTAGGGTGCGTTGAAGATAGGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAGAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTg                                                                                                                                                                                <  1:330087/251‑1 (MQ=255)
                                                                gagctgcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGGAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                                                                                                             <  1:566777/251‑1 (MQ=255)
                                                                     gcaagaAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                                                                                                          <  2:450319/249‑1 (MQ=255)
                                                                            aGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct                                                                                                                                                                 <  2:188528/251‑1 (MQ=255)
                                                                             gCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                                                                                                                <  2:503815/251‑1 (MQ=255)
                                                                                tCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTTGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                                                                                                              <  1:505929/250‑1 (MQ=255)
                                                                                 cGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACGCTCACg                                                                                                                                                             <  2:4021/250‑1 (MQ=255)
                                                                                       cccAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                                                                                                      <  1:249522/251‑1 (MQ=255)
                                                                                            ccttggttcacaACCGCCCGCGCTCAAAGAGCCTCTGTGATGGTGTGAAGAAGTTTTTTGACGACCTCACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCACAAGCG‑CTCACTT‑AA‑CGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                                                                                                                 <  2:125552/243‑1 (MQ=37)
                                                                                            ccGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACCTCCATTACCG‑CCAACTT‑AA‑AGGTTGTGGGATTTCGCCCGATTTTATTTTACGTGCCGATGGATTTCTTTCACATTAACTgttc                                                                                                                                                   >  2:14455/1‑246 (MQ=255)
                                                                                                  aGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                                                                                                           <  1:624504/251‑1 (MQ=255)
                                                                                                     acaacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                                                                                                        <  1:119982/251‑1 (MQ=255)
                                                                                                      caacaGCCCGCGCTAAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAAGTAGGGTGCGTTGAAGATATGCGAGCTCCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                                                                                                                       <  2:4485/251‑1 (MQ=255)
                                                                                                                cgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTCCCATAAGCG‑CTAAATT‑AA‑AGGGTGTGGGATTACGCCTGATATGATTTTAAGTGCCGGTGAATTATTCTCACGATAACTGGTCAGCAAATCTGGCCCCTATTGGt                                                                                                                             >  1:244270/1‑251 (MQ=255)
                                                                                                                                ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCT‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTg                                                                                                              <  1:71049/250‑1 (MQ=255)
                                                                                                                                           aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCa                                                                                                  <  2:523059/251‑1 (MQ=255)
                                                                                                                                           aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGCTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTAAATCCCAAAAATG‑CCAACTT‑AAAGGCTTG‑GGGGTTACCCCTGCTATGATTTTACGCTCCCAAGAAATATTCTCACGATAACAGGGTACTAATTCTTGCCCAAATTGGTAAGCACgtagaaataggtaattctgcc                                                                                                  >  1:95824/1‑236 (MQ=255)
                                                                                                                                                  tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt                                                                                            <  1:603409/250‑1 (MQ=255)
                                                                                                                                                               gACCCTCTAACCTCACCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTGACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGAGTTAGCGTGCCGATGAATTTCTCTCACGATAACTGGTCAGCAATTCTGGCCCATAGTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTa                                                                                <  1:620925/249‑1 (MQ=18)
                                                                                                                                                                                    gCCTGCCCGTGTGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGTGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGGGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGa                                                         <  1:611921/251‑1 (MQ=255)
                                                                                                                                                                                             tGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGCAAAGTGGCGGGGATAACTCCCATAAACG‑CTCACTT‑TA‑GGGTTGGGGGATTAAGCCCGATATGATTTAAAGTGCCGCTGAATTAATCTCACGATAACTGGTCAGCAATTCTGGCCCATATTTGTAAACCCGCAGAACTGGATAATTCCGCACCTAATGACGGGGCTCTAACCCGCCGCCGCGAAAATCGGGAtgctgccac                                                >  1:200007/1‑248 (MQ=17)
                                                                                                                                                                                             ggggCAGTCTTGGGTAAAAGTAGGGTGCGTTGAAGATATGCTAGCATCTGAAAAGTGGCGTGGATCACGCCCATAAGCG‑CGAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAAc                                                <  2:511982/250‑1 (MQ=14)
                                                                                                                                                                                             ggggCAGGCGTGGGGAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGAATGACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCTTGATGCTGCAAc                                                <  2:301891/250‑1 (MQ=17)
                                                                                                                                                                                                        gggTAAAATTAGGGGGCGTGGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTc                                     <  2:349855/251‑1 (MQ=255)
                                                                                                                                                                                                          gTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                                                                                                                                                                               <  1:95999/111‑1 (MQ=255)
                                                                                                                                                                                                               attAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGggct                              <  1:498974/249‑1 (MQ=255)
                                                                                                                                                                                                                  aGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctct                                                                                 <  1:467886/197‑1 (MQ=255)
                                                                                                                                                                                                                     gTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCCTAAACG‑CTAACTT‑AA‑GGGGTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCCCGATAACTGGTCAACAATTCTGGCCCATATTGGTAAACCCGAAGAAATGGATACTTCGGCACGTAAAGCCGGGGCGCTAAACCGCCGCCGCGGAATTAGTGATGCTGCAACTCTGCTACGTCTGGGGCCGGCttt                        >  1:202502/1‑250 (MQ=14)
                                                                                                                                                                                                                            gAAGATATGCTAGCACCTGTAAAGGGGCGGGGATCACCCCCAGAAGCG‑CATACTT‑AA‑GGGTGGTGGTATGACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGTTCAGCAATTCTGGCCCAGATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGTGCGCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTTGCTGTCGGccc                  <  2:64607/250‑1 (MQ=14)
                                                                                                                                                                                                                                      gAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGt       <  1:428753/251‑1 (MQ=39)
                                                                                                                                                                                                                                       aGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTc      <  1:578418/251‑1 (MQ=38)
                                                                                                                                                                                                                                        gctccTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGGCAGCAAGTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCGGGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCa     <  1:323527/248‑1 (MQ=255)
                                                                                                                                                                                                                                           ccTGAAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTTGTATTACGCCTTATATGATTTAACGTGCCGATGAATTTCTCTCACGATAACTGGTCAGCAATTCTTGCCCATATTTGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTa  <  1:554926/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                         
GCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCG‑CTAACTT‑AA‑GGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTA  >  NC_000913/15130‑15615

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: ‑