Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,277,454 T→A 10.6% Y40* (TAT→TAA 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,4540TA10.6% 165.3 / 2.3 81Y40* (TAT→TAAyhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base A (2/7);  total (32/49)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGA  >  NC_000913/3277226‑3277665
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                         
atatCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAATTCGAGACATTAAAGGGAAAAAATCCGGAAaataat                                                                                                                                                                                                                           >  2:227779/1‑229 (MQ=255)
 tttCCAGAAAAACCCGAAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATc                                                                                                                                                                                                                          <  2:624579/227‑1 (MQ=255)
 tgtccagacacaccagtaaaaaaCTCGTCCATTCAACATTCAAAAAGGAAAGAAACTTGTTGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATc                                                                                                                                                                                                                          <  1:500266/212‑1 (MQ=255)
        aacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTTATGATGAAATTGGCGACGACGAATAATGGATTTTCCACAAAGGTTTAATGGTTGGGCGCTATATACTCATCCCCGTTTTCAGGAAACCTACGACCCTTTATCTGCCGAATTCGAGACATTAAAGGGAAAAAAACCTAAAAAATAACAaaaaaa                                                                                                                                                                                                                   >  2:350867/1‑229 (MQ=255)
        aacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGTAAACCTACGACGCTTTAGTTGCCGAAGTCGGGACATTAAAGGGAAAAGATCCTGAAAATTATAAGAaaaa                                                                                                                                                                                                                   >  2:728125/1‑229 (MQ=255)
        aacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGACACCTACGACGCTTTAGTGGCCGAAGTCGAGACATTAAAGGGAAAAGATCATGAAAATTATCaggcaaa                                                                                                                                                                                                                   >  1:47639/1‑224 (MQ=255)
        aacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGCTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGGAAATTATCAGGaaaa                                                                                                                                                                                                                  >  1:62319/1‑229 (MQ=255)
          caacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGGTAATGGTTGGGGGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAGATCGAGACATTAAAGGGAAAAAATCCTGAAAATTATAAAAGaaaaa                                                                                                                                                                                                                 >  2:725051/1‑228 (MQ=255)
             cccGCACTAAACTGGTGCATTCAACGTTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGATATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGGCGAAGGCGGGAAATTAAAGGTGTAAGATCCAGAAAATCATCAGCAAAAAgcac                                                                                                                                                                                                              >  2:466420/1‑227 (MQ=255)
             cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGCGACATTAAAGGGAAAAGATCCTGAAAATTATCCGAGAAAAgcgg                                                                                                                                                                                                              >  1:298962/1‑227 (MQ=255)
             cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGaaaaacgg                                                                                                                                                                                                              >  1:405859/1‑225 (MQ=255)
             cccGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACATACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGCTCCTGAAAATTATCAGAGAAAAgccg                                                                                                                                                                                                              >  1:740379/1‑229 (MQ=255)
             cccGCAAAAAACTCGTCCATTAAACATTCAACAAGGCAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGGGCCGAAGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAACCCTACGATGCTTTAGTAGCCGCAGCCGAGACGTCAAAGGGAACAATTCCGGAAATATATCAGATCAAAgacc                                                                                                                                                                                                              >  2:445417/1‑226 (MQ=255)
             cccGAAAAAAACTGGTCCTTTTAACTTGAAAAAAGGAAAAAAACTTGTCGCAGGCTTGGACGTCAACATGGATGAGGGATTGGGCGACGACGAATAAGGGATTTTCCACAAAGGGTTAAGGGTTGGGCGTTATACGCTCATCCTTGTTTTCAGGAAACCTACGACGCTTTAGTGGCCAAAGTCGAGACATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAgccg                                                                                                                                                                                                              <  2:163477/229‑1 (MQ=255)
                 cAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCca                                                                                                                                                                                                           >  1:644764/1‑228 (MQ=255)
                  aaaaagCTCTGCCATTAAACATTCAAGAAGGAAAGAAATTTGCCGCGGGCATGGACGTAAAAATGGAGGATGAGATGGGCGCGGACGAATAAGGGATTTTCCACAAGGGGTTATTGGTTGGGCGCTATATGCTCATCCCTGTTTCAAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGAAAAAGATCTTGAAAATTATCAGAAAAAAGCCGCcaca                                                                                                                                                                                                         <  1:218036/229‑1 (MQ=255)
                  aaaaacCTCGCCCATAAAACATTAAACAAGAAAAGAAACTTGTCGCTGGCATGGAGGTCAAAATGGAGGATGGGTTGGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAAGGGTTGGGCGCTATATGCCCACCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATAAAAGGAAAAAGATCCTAAAAATTATAAGAGAAAAGCCGCcaca                                                                                                                                                                                                         <  2:644763/229‑1 (MQ=255)
                  aaaaaaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGCCGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAgcggcccc                                                                                                                                                                                                          >  1:116363/1‑226 (MQ=255)
                  aaaaaaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGCCGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCccca                                                                                                                                                                                                         >  1:576611/1‑229 (MQ=255)
                  aaaaaaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGGGAAAAGCCGCcaca                                                                                                                                                                                                         >  1:506308/1‑229 (MQ=255)
                            tCCATTCAACCTTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGGGAGTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTAATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCTTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTgg                                                                                                                                                                                               <  1:752511/229‑1 (MQ=255)
                            tCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAGTAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCCCAAAGGTATTgg                                                                                                                                                                                               >  2:521007/1‑229 (MQ=255)
                                ttcaacattcaacaAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGATTAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGCGACATTAAAGGGAAAAAATTCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTAATGgggga                                                                                                                                                                                           >  2:741620/1‑228 (MQ=255)
                                ttcaacattcaacaAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCgga                                                                                                                                                                                           >  1:540082/1‑228 (MQ=255)
                                               ggAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGATAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGAGGGAGACCATaaaaggat                                                                                                                                                                            >  1:75508/1‑224 (MQ=255)
                                                 aaagaacttggtgtggggaaTGGACGTAAACTTTGAGTTGGATATTGGCGCCGACTATTCATGGATTTCCCACAGAGGTTAAAGGGTTGGGCGTTATATTCTCATCCCTTGTTGCAGGAAACCAACGACTCTTTATTTGCCAAAGTCGATACATTAAAGGGAAAATATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTg                                                                                                                                                                         <  1:741621/210‑1 (MQ=255)
                                                   agaaaCTTTGCGCTGGCAGGGACGTAAACAGTGATGATGGGATGGGACACAGCGAATAAGGGATTTTCCACAAGGGGTTAATGGTTGGGGGTTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATAATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgag                                                                                                                                                                       <  2:412583/229‑1 (MQ=255)
                                                    gaaaCTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTGGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCATAGCAAAACCGCCACAAAGTTATggggggaaatcagaaaaacagagtagg                                                                                                                                                                       >  2:651520/1‑208 (MQ=255)
                                                    gaaaCTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACATACGACGCTTTCGTTGCCGAAGTGCCGCCATTACAGGGAAAAGATCCTGAAAATTAACAGCGAAAAGCCCCCACAAAGTAATTGGAGGGTGTcccctacaaggattggg                                                                                                                                                                       >  2:366521/1‑214 (MQ=255)
                                                     aatcTTGTCGCTGGCAAGGATGTCAACATTGATGATGAGGTTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGATATAATCTCATCCTTGTTTACAGGAAACCTACGACGTTTTAGTTGCCGAAGCAGAAACATTAAAGGAAAAAGATCCTGAAAATTAACAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgagga                                                                                                                                                                     <  1:149928/226‑1 (MQ=255)
                                                             cGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTCCGACGATTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCa                                                                                                                                                             <  1:53544/229‑1 (MQ=255)
                                                                       tgcGTAAACATGGATGAGGAGATGGACGACGACGAATAATGGCTTCTGCACAAGGGGTTAATGGTTGGGCACTTGACGCTATTCCCGGTTTTAAGGAAACCTACGCACATTAAGTAAGCGAAGTCGATACATTAAAGGAAAAAGACCATGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATcc                                                                                                                                                   <  1:342149/227‑1 (MQ=255)
                                                                                         gagaTTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAGAACATGAGGGGCAAATCACGGTCAATACATCATCAACGGCATTCCg                                                                                                                               >  1:624599/1‑229 (MQ=255)
                                                                                              tGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCAGCCACAAAGTTATTGGCGGTAGACAATAAAGGGATTTAGGGGCATAACACGGGCAAACcacacataccgtcattccctcctg                                                                                                                            >  1:290456/1‑207 (MQ=255)
                                                                                                  gcaccacaaaactttgtttttCCCAAAATGGTTAAGGGTTGTGCTCTAGATCTACATCCCTTTTTTATAAAACCCGACGACATTTTATTTCCAAAGGTCGAAACTTAAAAGGCAAAGGACCCTAAAAATTCTCTGAAAAATGCCGCCCCAAAGTTATTGGCGGTAGCCCATAAAGTGATTTAGGAGCATATCGCGGTCAATCCATCATCACCGGCATTCCGTCATGGcg                                                                                                                       <  2:386779/212‑2 (MQ=255)
                                                                                                  ctaggaagaATAATGGATTTTCCACAAAGGGTTAATGGTTGGGTGATAAATGCTCACCCTTGTTTGCAGGAAACCTACGCCGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCa                                                                                                                       <  2:160865/222‑1 (MQ=255)
                                                                                                   tgagtctaataatGGGTTTCCCCAAAGGGGTTAATGGTGGGCGCGTAAGCGCTCATCCCTGTTTCCAAGAACCCGACGACGATTTAGTTGAAGAAGTCGAGAAATAAAGGGAAAAAGAGCCTGAAAAGTACAAGAGAAAAGCCGCCACAAAGGTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGAATTCCGTCATGGCa                                                                                                                       <  1:350867/221‑1 (MQ=255)
                                                                                                   cggcgaccaattagggttttCCCCAAAAGGTTGAATGGTTGCGCGTTATCGGGTCATCCCTCTTTCCAGGAAACTAACGGCGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCa                                                                                                                       <  1:366521/213‑1 (MQ=255)
                                                                                                        ccaataagGGATTTCCCAAAAGGGGTTAATGGTTGGGCGCTTAATGCTAACCCTTGTTTTCAGGAAACCTACGACGTTTAAGTACCCAAAGCCGAGAATTTAAAGGAAAAAGACCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGt                                                                                                                <  1:202446/227‑1 (MQ=255)
                                                                                                         aaataatGGATTTCCCACAAAGGGTTAAGGGTTGGGCGCTATATGCTCATCCGTTTTTTCAGGAAACTTACGACGTTTTAGTTGCCGAAGCCGAGAAATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGtt                                                                                                               <  1:521007/228‑1 (MQ=255)
                                                                                                               ttgATTTTCCAAAAAGGTTAAAGGGTTGGGCTCTATATGCTCATCCCTGTTTTCAGGAAACCTGCGACGTTTTAGTTGCCAAAGTCGAAAAATTAAAGGGAAAATACCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGct                                                                                                          <  1:647480/226‑1 (MQ=255)
                                                                                                                  aTTTCCCAAAAGGGGTTAAGGGTTGGGCGCTATATGCTCATCCCTGTTTTAAGGAAACCTCCGACTCTTTAGTTGCCGAAGTCGAGACATAAAAGGAAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTgg                                                                                                      <  1:476182/229‑1 (MQ=255)
                                                                                                                    tttCCCAAAGGGGTTAATGGTGGGGGTATTGATGCACACCCTTTTTTTCAGAAAACCTACGCTGTTTTAGTTGCAGAAGACAGAACATAAAAGGAAAAAGACCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                    <  2:523544/229‑1 (MQ=255)
                                                                                                                    tttCCAAAAAGGTTTAAGGGTTGGGTGTTTTATCTCTACCATTTGTTCAAGGAAACCTCCACGTATTTACGTCCAGAGGTGGAAACATTAAAGGAAAAAGACACTGAAAATTATCAGAAAAAGGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCTTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                    <  2:713590/229‑1 (MQ=255)
                                                                                                                    ttcccACAAAGGGTTGATGGTTGGGCATTATATCCATACACTGTTTTGCAGGAAACCTACGACGCTTAATTTCCCGAAGTCGAGACATTAAAGGGAAAAGTTCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                    <  1:466420/226‑1 (MQ=255)
                                                                                                                    ttcaccaaaggggTTAAGGGTGGGGCGTTAGAGGCACCTCACGTTTTACAGGAAAGCGAGGAGTATTTGTTTGCCAAAGTGGAGACATTAAAGGGAAAAGTTCCTAAAAATCATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCTCGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCTTTAGGCTCTGGGa                                                                                                    <  2:704017/219‑1 (MQ=255)
                                                                                                                    gtcaccaaaaGGGTTAATTGTGGGGCCTTTAGTCCCCAGCCCTTTTATAAGAAAACCAGCGACTGTTTGGTTGCCGAAGGCGACAAATAAAAGGAAAAAGATCGTAAAAATTATCAGAAAAAGTCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                    <  2:506308/222‑1 (MQ=255)
                                                                                                                     tccctcaAAGGTGTAGTGGTTGGGTGTTACATTCTCAGCCTTGTTTCAAGAAATCCGTCGACCTTTTAGTTGCGCAAGTCGAGACATAAAAGGGAAAGGCTCCAGAAAATTATCAGAAAAAGTCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                    <  2:75508/223‑1 (MQ=255)
                                                                                                                           aaaGGGTAAATGGTTGGGCGTTATATCGTAAGCCTTTTTTCAAGGAAACCAACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGAAAAATATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATaa                                                                                             <  1:301614/229‑1 (MQ=255)
                                                                                                                               ggTTAACTGGTTGGGTGTATTATGTTCATCCCTTTTTTCAGTAAACCTACGACGTTTAAGTTGCCGAAGTCGAGAATTTAAGGGAAATAGGTCCTGAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATTGCAAGTCGTTAGGCTCTGGGAAAAATAAAGa                                                                                          <  1:35232/229‑1 (MQ=255)
                                                                                                                                    aaTGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGTATTGAGGAGCATAATACGGTCAAACCATCATCACCGGCCTTCCCGCCTGGGCAGATGTTTGGGCTCTGGGCAAATAAAAACTggg                                                                                     >  1:524614/1‑228 (MQ=255)
                                                                                                                                     atgggTTGGGCGCTATATCCTCAACCGTTTTTCAAGGAAACCAACGCCGCTTTGGTGGCCGAATCCGAGACATTAAAGGAAAAAGATCCGAAAAATTACCAGACAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGt                                                                                     <  2:624598/226‑1 (MQ=255)
                                                                                                                                      gtGGTTGGGCACTATCAGCTACTCCCTGGTTGAAAGAACCATACGCAGCTTGAATTGCCGAAGGCGAGACAGTAAAGGGAAAAGGTCCTAAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTc                                                                                   <  1:34831/228‑1 (MQ=255)
                                                                                                                                        ggTTGGGCTCTATATGCTCATCCCTGTTTTCAGAAAACCTACGAGGTTTTAGTTGCAGAAGTCGGGACATTAAAGGAAAAAGACCTTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTAATGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCAc                                                                                 <  1:728126/229‑1 (MQ=255)
                                                                                                                                             aaggacccgatGCTTTTCCCTGTTTAGAGGAAACCAACCACGATTGATGTGGCAAAGTCGAGAAATTAAAGGAAAAAGCCCCAGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGt                                                                             <  2:253052/220‑1 (MQ=255)
                                                                                                                                              gccgttataTGTACACCCCTGTTTTAAGGAAACTTACGACATATTAGCTGCAGAAAGGGAGAACTAAAAGGGAAAAAGTTCAAAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTa                                                                            <  2:752300/224‑1 (MQ=255)
                                                                                                                                                         gccaaaccttGTTTTAAGGAAACTTACGACGCTTAAGTTCCAGAATCCGAGACTTTAAAGGGAAAAGATCCTGAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                  <  1:219269/220‑1 (MQ=255)
                                                                                                                                                         gCTCACCCGTGTTTTAAGGAACCCTTCGCTGATTAAGTTGCCGAAGGAGAGCCATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGTAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                  <  1:279643/229‑1 (MQ=255)
                                                                                                                                                          gacATCACTGTTTTGAGGAAACCTACGACCTTTAGGTTGCCAAGGTCGAGACATTAAGGGGAAAATACCCTGAAAATTTTCAGAAAAAGTCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                  <  1:174348/226‑1 (MQ=255)
                                                                                                                                                              tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAAGAATTCAGAGGCTTAAACGGTGAATCCATCATCACCCGCAATACGTCAAGGCAAAGCGCTTGGCTCCGGGTAAAAAAAAAAATCGTCACGGGGAAAAATTTGGGATGGt                                                             >  2:219269/1‑229 (MQ=255)
                                                                                                                                                              tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGGTGGCCATAAAATGATTTAGGAGCATATCACGGTCAATCCATAATCACCGGCAATCCGTACCGGGAAGTCGGTAAGGCATGGGGAAAATTAAAAACGGTTACGGGGAAAAATTgggtgtggt                                                             >  1:108607/1‑222 (MQ=255)
                                                                                                                                                                ccTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAAGATTGAGGAGCAAATCACGGTCAAACCATCATCACCGGCATTCCCTCATGGCAAGGCGTTAGGCTCTGGGGAAAATAAAAACTGGTCACGGGGCAAATTTTGTGCTTGt                                                             >  1:529410/1‑227 (MQ=255)
                                                                                                                                                                  tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAATGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTTCGTCCTGGCAAGTCGTTAGGCCCTGGGGAAAATAAAAAATGGGCACGGGGAAAAATTTGGGCTgtgtct                                                          >  1:253052/1‑223 (MQ=255)
                                                                                                                                                                  tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAAGGTTGAGGAGCATATAACAGGCAATCACTCATCACCCGCCTTCCCTCATTGGAAGACGTTTGGGGCTGGGGAAAAAAAAAACTAGTACCCGGGGAAAATTTGTGATGGGCGtt                                                         >  2:569612/1‑229 (MQ=255)
                                                                                                                                                                  tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCCTAAAGAGATTTAGGGGCATATCACGCGCAATCCAGCATCCCCCCCAATCCGTCATTGCAAGACGTTAGGGCCTGGGTAAAAAAAAAACAGGGCCCGGGGaaaaatagtggttggtgggt                                                         >  2:532081/1‑213 (MQ=255)
                                                                                                                                                                    ttttCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGGCATAAAATGATTGAGGAGCATAACACGCTCAATACATCATAACCGGCATTCCGTCATAGGAAACCGTTAGGCTCTGGGGAAAAAAAAAACTGGGCACCGGGAACAATTTGTGCTgggtgctat                                                       >  2:711865/1‑222 (MQ=255)
                                                                                                                                                                    ttttAAGGAAACCTACGACGTTTTAGTTGCGGAAGTCGGAACATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                         <  2:246773/135‑1 (MQ=255)
                                                                                                                                                                       tCAGGAAACCTACGACGCTTTAGTCGCCGAAGTCGAGACATTAAAGGGAAAATATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGt                                                    <  1:507552/229‑1 (MQ=255)
                                                                                                                                                                         taaGGAACCCTAGTACGTTTAAGTTGCCAAATGCAGGAAATAAAGGGGAAAAGACCCTAAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGt                                                    <  1:227779/227‑1 (MQ=255)
                                                                                                                                                                         taaGGAAAACTACCGCGTTTTAGGTGCCGAAGGCGAAAAATTAAGGGGAAAAGGTTCTGAAAAATAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGt                                                    <  2:290456/227‑1 (MQ=255)
                                                                                                                                                                                         cgttttAGTTGCCGAAGTCGAGAAATTAAAGGGAAAAGACCCGGAAAATTATCAGAGAAAAGCCGGCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTta                                      <  1:77360/225‑1 (MQ=255)
                                                                                                                                                                                           tgtttatgtggcGAAGTCGAAAAGTAAAAGGGAAAAGACCCTGAAAATTATCAGAGAAAGGCCGCCACAAAGTTATTGGCTGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTta                                      <  1:110922/216‑1 (MQ=255)
                                                                                                                                                                                                        aagaggagaAATAAAAGGAAAAAGATCCGGAAAATAATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGaa                      <  2:740378/224‑1 (MQ=255)
                                                                                                                                                                                                           tCGAGAAATTAAAGGAAAAAGATCCTGAAATTTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGt                   <  1:468815/229‑1 (MQ=255)
                                                                                                                                                                                                             agaaCATTAAGGGAAAAAGATCGTGAAAATTAACAGGAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                 <  1:711866/226‑1 (MQ=255)
                                                                                                                                                                                                              agaCATTAAAGTGAAAAGACGCTGAAAATTAAAAGAGAAAACCCGCCCCAAAGTTATTGGCGGTAGTCCATAAAGTGAATGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                 <  1:723252/228‑1 (MQ=255)
                                                                                                                                                                                                                   taaaaGGGAAAATATCCAGAAAATGAAAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGCCCATAAAGGGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCAtt            <  2:529410/226‑1 (MQ=255)
                                                                                                                                                                                                                   taaaaGGAAAAAGACCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCAtt            <  1:636275/226‑1 (MQ=255)
                                                                                                                                                                                                                    taaagaggaaaaagccgaaaaaaTTAGCAGAGAAAAACCGCCACAAAGTTATTGGGGGGAAGCCATAAAGGGATTGAGGAGCATATCACGGGCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCGCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTc           <  2:576610/210‑1 (MQ=255)
                                                                                                                                                                                                                    tAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTTGTGCATCAAAGGATTTAGGGGGCAATCACGGACAAATCATCCATCCCGGCGTTCCGTCTTTGCAAGTCCGTAGGCTTGGGGAAAAAACAAAAATGGGCACGGGGAAAAATTTGTGGACGTCGTTGACGTCCATTATTTTTTTTTAGTTGAAAAAGGAAAAgaaccaagta           >  2:495101/1‑220 (MQ=255)
                                                                                                                                                                                                                            aaaGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGACCCTAAAATAGTTGGGGAGCAAATCACGGTCAATCCATCCTCACCGGCATTCCTTCATGGCGAGTAGTTAGGCTCTGGGCAACATCAAGACTGGGGACGGGGAAACATTTGTGCTGGGCGTTATCGTTTCTTCTTTCGTTTTAATGAAAAAAAGAAAAGTATCATTCTGggatggg  >  1:421645/1‑228 (MQ=255)
                                                                                                                                                                                                                                    cTGAAAATTAGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTg                                                                  <  1:65512/159‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                         
ATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGA  >  NC_000913/3277226‑3277665

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