Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,277,507 A→G 17.1% E58G (GAG→GGG)  yhaV → toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,277,5070AG17.1% 159.7 / 3.5 67E58G (GAG→GGG) yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (11/42);  new base G (10/2);  total (21/46)
Fisher's exact test for biased strand distribution p-value = 8.68e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAG  >  NC_000913/3277269‑3277741
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
acaaGGAAAGAAACTGGTCGCTGGTGTGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTCCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGAAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgagga                                                                                                                                                                                                                                              <  2:454325/239‑1 (MQ=255)
  aaGGAAAAAAACTTGTCGGGGGCATGGACGCAAAAATTGAGGAGGAGGTGGGGCAAGAGAAATAATGGATTTTCCAAAAAGGGTTAATGGTTGGGTGCTTTATGCTCTCCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCATGAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGc                                                                                                                                                                                                                                            <  1:240794/239‑1 (MQ=255)
    ggAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAATGATTTAGGGGcct                                                                                                                                                                                                                                          >  2:67265/1‑237 (MQ=255)
    ggAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTACATAAAATGATTTAGGGGCat                                                                                                                                                                                                                                          >  2:752415/1‑239 (MQ=255)
            cctgggtggtggatgtgACGTCAACATTGGTTATGGAGTGGGCGACGACGAATAATTGATTTTCCAAAAAGGGTTAATGGTTGGGAGCTATATGCTCACCCCTGTTTTCAGGAAACCTACGACGCTTTGGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                                                                                  <  1:777067/223‑1 (MQ=255)
            atttttccttGGCATTGAGGTCAAATTTGATGTTGGGTTGGGCGCCAAAAAATAAGGGTTTTTCCACAAAGGGTAAAGGGTTGGGCCTTATATGCTCATCCCTGTTTGCAGAAAACCTACGACGCTTTAGTTGCCAAAGTCGAGACATTAAAGGAAAAAGATCCAGAAAATTATAAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                                                                                  <  1:545023/230‑1 (MQ=255)
               tGTCGCTGGCAAGGACTAAAACATTGATGTTGAGATTGGCGCCAGCAAATAATGGATTTTCCACAAGGGGTTAATGGTTGGGCGTTATATGCTCCTCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCaa                                                                                                                                                                                                                               <  2:115085/239‑1 (MQ=255)
                    cTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGATAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTATTCCCTAAAGTGATTGTGGATCacattccgtgcaaaccct                                                                                                                                                                                                                          >  2:216325/1‑225 (MQ=255)
                                 cAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCCGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGAGAATTAGGGGCCTAACACGGGCAAACCATCATCACCgccct                                                                                                                                                                                                              >  2:1006119/1‑234 (MQ=255)
                                 cAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAAGGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCAAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCAt                                                                                                                                                                                                              <  2:258093/238‑1 (MQ=255)
                                             tgggaTTGGCGACGACGAATAATGGATTTCCCAAAAGGGGTTAAGGGTTGGGCGCTATATGCCCATCCCTGTTTCCAGGAAACCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGACCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCa                                                                                                                                                                                                 <  2:83978/235‑1 (MQ=255)
                                              gagtttggccgcgtggaATAGTGGATTTTCCCCAAGGGGTTAATGGTGGGTCGTTATATGCTCATCCCGGTTTGCAGGAAACCTACGACGCTTTAGTTGCCGAAGGCGAGACATAAAAGGAAAAAGATCCTAAAAATTAACAGAGAAAGGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCaa                                                                                                                                                                                                <  1:573393/223‑1 (MQ=255)
                                              gaggtgggcggcGCCAAATAATGGATTGTCCACAAGGGTTAAAGGGTGGGGGGTTAAAGGACTCTCCTGGTTTTCGTGAATCTTCCGACTTTTTATTGGCGTAAGTCGAGACATTAAAGGGAAATGATCCAGAAAATTATCAGGAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCaa                                                                                                                                                                                                <  1:752415/232‑1 (MQ=255)
                                                   tGGCGAAGAAGAATAAGTGATTCTCCAAAAAGTGTAAATGTTGGGGGTTTATATTCTCCTCCAGTTTTTCAGAAATCCTACGCCGGTTGAGGTGCCAAAGCCGAGAATTTAAAGGGAAAAGTTCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGtt                                                                                                                                                                                          <  2:58098/239‑1 (MQ=255)
                                                              aataatGTATTTACCAAAAAGGGTTAAGGGTTGGGCGTTATATGCCCACCCCTGTTTTCAGAAACCCGACGACGTTTTAGTTGCCGAAGTCGAAACATTAAAGGGAAAAGATCCTAAAAATTACCGGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                                               <  2:449470/239‑1 (MQ=255)
                                                              aataatGGATTTTCCAAAAAGGGTTAATGGTTGGGCGTTATATGCTCATCCTTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                                                                                  <  2:740106/189‑1 (MQ=255)
                                                               atatgggATTTCCCACAAGGGGTTAAGGGTTGGGCGCTATATGCGCATCCCTGTTTCAAGGAAACCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                                               <  2:962341/233‑1 (MQ=255)
                                                               ataatGGATTTCCCACAAAGGGTTAATGGTTGGGCCGTATGGGCTGCTCCCTGTTTTCAGGAAACCTACGACGTTTTAGTTGCCGAAGTCAAGAAATAAAAGGGAAAAGACCCTAAAAATGATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                                               <  2:263335/238‑1 (MQ=255)
                                                                    tgtttttTCCCAAAGGGTTTAATGGTGGGCGTCATTATGCACTCCCCTGTTCTCAGAAAACCTACCACGGTTTAGTGCCCGAAGACGAGACATAAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAAt                                                                                                                                                                          <  1:667045/236‑1 (MQ=255)
                                                                               aaaGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCATAAGAGATTTAGGAGACAATACCGGGCAATCCATACTCACCCGCATTCCGTGCTGGACAGTAGGTAAACTCTGGGAaaaaaaaaaaagtgacc                                                                                                                                                              >  1:378232/1‑230 (MQ=255)
                                                                               aaaGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGCATGTCCAAAAAAGATTTAGGGGGAAATCACGGTCAATACATCATTAACGGCATTCCGTCCTGGGCAAGTGGTAGGGACTGGGCaaaaaaaaaaacgagga                                                                                                                                                              >  1:979477/1‑230 (MQ=255)
                                                                                      ttAAGGGTTGGGGGCTATATGCTAATCCCTGTTTTAAGGAAACCTTCGACGCTTTAGTTGCCGAAGTCGAGAAATTAAAGGAAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTa                                                                                                                                                       <  2:516769/239‑1 (MQ=255)
                                                                                       tAATGGTTGGGAGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTGTATTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTAGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTa                                                                                                                                                       <  2:163958/238‑1 (MQ=255)
                                                                                       tAAGGGTGGGGCGCTAATTGCCCATCCCTGTTTCAAGGAACCCTACGACGCTTTAGTTGCCGAAGTCGAGAATTTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTa                                                                                                                                                       <  2:949228/238‑1 (MQ=255)
                                                                                       aaagggTTGGGCGCTATAGTCTCATCTCTGTTTGCAGGAAACCTACGACGCTTTAGGTGCCGAAGTCGAGACTGAAAAGGGAAAAGATCCTGAAAATTTGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTa                                                                                                                                                       <  2:385797/234‑1 (MQ=255)
                                                                                              tGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTGGTCCACAAAAGTGTTGAGGGGACAATACCGGGACATCCACCCTCCCCCGCATTCCCTAAGTGGAAGACGTTAGGCTCTTGGAAAAAAAAAAACAGGACACGGGAaaaaaatgt                                                                                                                                               >  1:167854/1‑234 (MQ=255)
                                                                                              tGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAATGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTTCGTCATGGCAAGTCGTTTGGCTCTTGGGAAAATCAAGACTGGTCACGGGGCAAAATTgg                                                                                                                                               >  2:803006/1‑239 (MQ=255)
                                                                                                  gggtttgatGCCCACCCCTGTTTTAGGAAACCCTACGACGTTTTAGTTGCCGAAGACGAAAAATTAAAGGAAAAAGATCCTGAAAATTAAAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                                             <  2:351388/231‑1 (MQ=255)
                                                                                                  ggcgtttttgtgtcccccgtGTTTTCAGAAAACCTACTACGCTTTAGTGGACGAAGTCGAGAAAATAAAGGGAAAAAACTCTGAAAATTATGAGATAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                                             <  1:803006/220‑1 (MQ=255)
                                                                                                  gccttttttGCTCATTCTTTTTTTCAGAAAACCTGCGATGTTTTAGTTGCAGAAAGCGAGACATTAAAGGGAAAAGATCCGGAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                                             <  1:429919/230‑1 (MQ=255)
                                                                                                  ccgttataTTCTCATCCCTGTTTACAGGAAACCTACGACGTTTTAGTTGCCGAAGTCGGGATATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                                                                                             <  2:198836/234‑1 (MQ=255)
                                                                                                        tatGCTCATCCCTGTTTTCAGGAAACCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCAcc                                                                                                                                                                                                                   <  2:391027/164‑1 (MQ=255)
                                                                                                                aTCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGGGGCTATCACGGTCAATCCATCATCACCGGCATTCCGTCATTGCAAAGCGTTAGGGCTTGGGGAAAATAAAAACTGGCACCGGGGAAAAATTgtgggtgtccggtttg                                                                                                                                 >  2:260639/1‑222 (MQ=255)
                                                                                                                 tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATc                                                                                                                                                                                                                             <  1:688930/146‑1 (MQ=255)
                                                                                                                 tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAATGATTGAGGAGCATCTCCCGGTCAATc                                                                                                                                                                                                                             >  2:688930/1‑146 (MQ=255)
                                                                                                                 tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGGATTGAGGGGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAAAATGGTCACGGGGAAAAATTGGGGCTGGTGGTTATCg                                                                                                                                >  2:370340/1‑237 (MQ=255)
                                                                                                                   ccTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGGGGCTATCACGGTCAATCCATCATCACCGGCCTTCCGTCATGGCAAGTCGTTAGGGCCTGGGAAAAATAAAAACGGGGCCAGGGACAAAATTTGGGCTGGCGGTTTTCGtc                                                                                                                              >  2:548174/1‑237 (MQ=255)
                                                                                                                     tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTGCATAAAAGGATTTAAGGACATATAACCGCACATACATCATCCACAGCAATCCGTCGTGGGAAGGAGGTCGGGTCGGGGGAAAATAAAAATAGGTCACGGGGAAAAattttggcgggtggtatttttacgtt                                                                                                                          >  1:106103/1‑216 (MQ=255)
                                                                                                                     tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGGGTTTTGGAGCATCACCCGGGTAATCCCTCCTCAACGGCATTCCGTCCTGGGAAGTCGTTAAGCGTTGGGGGAAATAAAGACAGGACACCGGGAAAATATTGTGTTGTTCggattagtccct                                                                                                                           >  1:565562/1‑227 (MQ=255)
                                                                                                                            aGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGGTCTGGGGAAAATAAAAACTTGGCACGGGGAAAAATTTGTGCTTGGCGTTTTCGGTCCTTCTTTCGt                                                                                                                   >  2:140652/1‑239 (MQ=255)
                                                                                                                                 aacctcagtagttttGGTGGCCGAATGGGAGAAATAAAAGGAAAAAGATCCTGAAAATAAACGGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTtata                                                                                                               <  1:140652/227‑1 (MQ=255)
                                                                                                                                       cgccgCTTGAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGCTCCTAAAAATTACCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGaaa                                                                                                         <  2:461046/236‑1 (MQ=255)
                                                                                                                                             tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGGGTTGTGGAGGATATAACGGTCAATCCATCCTCAACGGCATTCCGTACAGGCCGGTCGTTAGGCTCCTGGGAAAATAAAGACTGGTGCCGGGGAAAATTTTGTGCGGGGTGGTATGGTCTCTTCTTTCTTTAAAAAGaaaaaaaaa                                                                                                   >  1:11759/1‑235 (MQ=255)
                                                                                                                                             tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCAATAAAATGATTTTGGGGGAAATAACGGTCAATCCACCATCACCGGCATTCCGTCATGGCAAGTACTTAGGCTCCGGGGAAAAAAAAAAATGGTCACGGGGCAAAATTGTGGCTGGTCGTTTTATGCTCTTTTTTTGTTTTAAGGaaaaaaaaa                                                                                                   >  1:696866/1‑235 (MQ=255)
                                                                                                                                                   taataagagacagaCATAAAAGGAAAAGGATCAAAAAAATTAACAAAAAAAAGCCGACACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATAGCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                                                                             <  1:170519/228‑1 (MQ=255)
                                                                                                                                                   tGCCGAAGTAGAGAAATTAAAGGGAAAAGATCCTAAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                                                                             <  2:754814/239‑1 (MQ=255)
                                                                                                                                                    cccGAAGGCGAGAAATTAAAGGGAAAAGTCCCGAAAAATTATGAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtcat                                                                                            <  2:979477/238‑1 (MQ=255)
                                                                                                                                                              ctaCAATAAGGGAAAAAGAGCTTGAAAATTAACAGAGAAAAGCCGGCCCAAAGTTATTGGCGGTAGCCCAAAAAGTTATTTAGGAACATATCACGGTCAAGCCATCATCACCGGCAGTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGgga                                                                                  <  1:491648/237‑1 (MQ=255)
                                                                                                                                                               aaCATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGgga                                                                                  <  2:11759/237‑1 (MQ=255)
                                                                                                                                                                       aGGGAAAAGATCCCTAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAAc                                                                         <  2:106103/239‑1 (MQ=255)
                                                                                                                                                                           aaaaGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGAGATTGAGGAGCATATCACGGTCAATCCATCATCCCCCGCATTCCGGCCTGGCAAGTCGTTAGGCTCTGGGGAAAATAAAGACTTGTCACCGGGAAAATTTGGTGCTGGTCGTTATCGCCTCTTCTTTCGTTTTAAGGAAAAAAAGAAAGtcatccatttggg                                                                                   >  2:170519/1‑217 (MQ=255)
                                                                                                                                                                                 tCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAaca                                                               <  2:696866/239‑1 (MQ=255)
                                                                                                                                                                                        atttATAAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTgcg                                                        <  2:249122/237‑1 (MQ=255)
                                                                                                                                                                                               aggaaaaaGCCGGCCAAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTAc                                                 <  2:576971/235‑1 (MQ=255)
                                                                                                                                                                                               aaagaAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGGCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTAc                                                 <  2:382170/237‑1 (MQ=255)
                                                                                                                                                                                                       gccgccACAAAGTTATTGGCGGTAATCCATAAAAAGATTGAAGAGGATATCACGGTCAAACCATCATCCACGGCATTCCGGCATAGGAAGACGGTAAGGTCTGGGGAAAAAAAAAACTGGGTACCGGGAAAAATTTGTGCTGGTCGTTATCGTCTCTTCTTTCTGTTTAGTGAAAAAAAAAAAATGAACATTTTTGGATGGGTGAACGATGAAAaaacaacacgacacaacggtaaaaa                                         >  1:135203/1‑218 (MQ=255)
                                                                                                                                                                                                           ccACAAATTTTTTGGCGGGTGTCCAGAAATGGATTTAGGGGAATATCACGGACAATTCAGCATCACAGAGATTCCGTCATGGCAAGTCGTTAGGGTGCGGGAAAAATAAAGACTGGTCACGGGTAAAAGTTGGTGCTGGTCGTTATCGTCTCTGCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                                     <  1:985956/239‑1 (MQ=255)
                                                                                                                                                                                                           ccACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                                     <  2:135203/239‑1 (MQ=255)
                                                                                                                                                                                                            cacaAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                                     <  2:565562/238‑1 (MQ=255)
                                                                                                                                                                                                             cccaaGTTTTTTGCGGTAGTCCATAAATTGATTGGGGTGCATATACCGGGCAATCCATCATAGCCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa                                   <  1:370340/236‑1 (MQ=255)
                                                                                                                                                                                                             acaATGTTATTGGCGTTAGTCCATAAAGTGAGTGAGGCGCATAGAAAGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTGTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa                                   <  2:391563/239‑1 (MQ=255)
                                                                                                                                                                                                             acaAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa                                   <  2:47166/239‑1 (MQ=255)
                                                                                                                                                                                                             acaAAGTTATTGGCGGTAATCCATAAAGTTATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa                                   <  2:14429/239‑1 (MQ=255)
                                                                                                                                                                                                                  tttATTCGCGGTCGGCCATAAAGATATTGTGGGGTATCTCACGTTCATTCCATCATCACCGGCATTCCGTAATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTTGTGCTGGTCGTTATCGTCTCTTCTGTCGTTATAGTGAAAAAGACAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACGCTGCGCACCTACGGTAAAAAAACAGATGCCt                              <  1:469019/238‑1 (MQ=255)
                                                                                                                                                                                                                                  ccATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAAGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCa               <  2:378232/239‑1 (MQ=255)
                                                                                                                                                                                                                                           tGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGACAAATAAAGACTGATCACGGGTAAAATTTGGTGCTGGGCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTTTGGGATGGCTGAACGTAGCAAAAACTCAGCGAACCTACGGATAAAAACAACATCCCAATACAGTAttaccaaaacagtaa       >  1:107123/1‑226 (MQ=255)
                                                                                                                                                                                                                                              ttGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAAAACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGAATTCCTCAAAATTTTaaaa    >  2:527402/1‑239 (MQ=255)
                                                                                                                                                                                                                                               tGAGGAGCATATCACGGTCAAACCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaaaaa   <  2:878619/239‑1 (MQ=255)
                                                                                                                                                                                                                                                gaggagCAAATCACGGGCAATCCATCCTCCCCGGCATTCCGTCATGGCAAGGCGTTAGGCTCTGGGAAAAATAAAGGCTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAag  <  2:177277/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
ACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAG  >  NC_000913/3277269‑3277741

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