Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,277,500 T→G 17.7% V55G (GTG→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,5000TG17.7% 146.4 / 6.7 76V55G (GTG→GGG) yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (11/46);  new base G (9/5);  total (25/51)
Fisher's exact test for biased strand distribution p-value = 1.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCA  >  NC_000913/3277269‑3277730
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
caacaaGGAAAAAAACTTGTCGCTGGCATGGACCTCAAAATTGATGATGAGGTTGGGGACGACAAATAATGGTTTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgag                                                                                                                                                                                                                                      <  2:837743/239‑1 (MQ=255)
caacaaGAAAAAAATATTGTCGTGGACTGGGACTTCAACATGGATGAGGAGAGTGGCGACGAAGAATAATGGTTTTTCCACAAAGGGTTAATGGTGGGGTTTTTTATGCTCATCCCTGTTTTCAGGAAACCAACTACGCTTTAGTTGCCGAAGTGGAAACATAAAAGGGAAAAGATCCTAAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgag                                                                                                                                                                                                                                      <  2:469887/239‑1 (MQ=255)
caaaaaGGAAAAAAACTTGTCTCTGGAATGGACGTCAACAGTGAGGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTGGGGCGCTATAGGCTCACCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGAAATTAAAGGGAAAAGACCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgag                                                                                                                                                                                                                                      <  2:866713/239‑1 (MQ=255)
aaacaaGGAAAGAAATTTGTCGCGGGAGAGGACGACAACATTGAGGAGTAGTTTGGCGACCAGAAATAAGTGATTTTCCAAAAAGGTTTAATGGTTGGGCGTTATATGCCCATCCCTGTTTCAAGGAAACCTCCGACGCTTTAGTTGCCGAAGTGGATACATAAAAGGAAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTgag                                                                                                                                                                                                                                      <  2:612668/238‑1 (MQ=255)
    aaGGAAAGAACCTTTCCGCTGGCAGGGACCTAAAAATTGGTGGAGAGATTGGCGACGAAGAATAATGGATTTTCCACAAAGGGTTAAGGGTTGGGCGATATATGCTCACCCCTGTTTCCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGc                                                                                                                                                                                                                                  <  2:1015247/239‑1 (MQ=255)
      ggAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGAtttgaggggat                                                                                                                                                                                                                                >  1:660743/1‑230 (MQ=255)
           gaaaCTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATAGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATAGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAACGGGAAAAGATCCTGCAAATTATCCGAGACAAGCCGCCCCAAACTTATTGGCGGCAGGCCATCAAATGATTgagtagacatttcc                                                                                                                                                                                                                           >  1:971723/1‑231 (MQ=255)
               cttgtgctgggCATGGACGTAAAAATGGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGTTATATGCTCATCCCTGTTTTCAGGAACCCTACGACCCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTc                                                                                                                                                                                                                       <  1:575341/230‑1 (MQ=255)
                   gggCTGGGATGGGACTTCAACTGTGATGGTGGGGTGGGCGCCAAAGAATAAGGGATTTTCCACAAAGGGTTAATGGTTGGGCGTTATATGCTCATCCCTGTTTGAAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGAAATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAAt                                                                                                                                                                                                                    <  2:221934/237‑1 (MQ=255)
                            tGGGCGTAAACATGGATGATGTGTTTGGCGACGACAAAAAATGGATTTTCCACAAGGGGTTAATGGTTGGGCGCTAAATGCTCAGCCCTGTTTACAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCAc                                                                                                                                                                                                          <  2:732899/239‑1 (MQ=255)
                                    caaattttatgatgcAATTGGCGACGACGAATAATGGATTTTCCAAAAAGGGTTAATGGTGGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAATTCGAGACATTAAAGGGAAAAGTCCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCAtt                                                                                                                                                                                                   <  1:113648/231‑1 (MQ=255)
                                    cAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAACCCGCCACAAAGTTATTGTCGGTAGTCCATAACAAGATTTAGGGGCAATACACGGTCAATCCACCATCACCTGcct                                                                                                                                                                                                    >  1:952034/1‑236 (MQ=255)
                                                tgaGATTGGAGACGGCGAATAAGGGATTTTCCACAAGGGGTTAATGGTTGGGCGCTATAGGCTCATCCCTGTTTTCAGGAAACCTACGACGTTTTAGTTGCCGAAGGCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCa                                                                                                                                                                                        <  1:738924/238‑1 (MQ=255)
                                                 gatagtGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAg                                                                                                                                                                                      <  1:589002/234‑1 (MQ=255)
                                                 gagaTTGGCGACAACAAAAAAGGGATTTACCACAAGGGGTTAATGTTTGGGCGTTATATGCCCATCCCTGTTTACAGGAAACCTACGACGCTTTAGTTGCCGAAGGCGAGACATTAAAGGGAAAAGATCCTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAg                                                                                                                                                                                      <  1:680850/239‑1 (MQ=255)
                                                      tGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCAAAAAAGGCTTGAGGAGCCAATAAAGGTCAATACATCATCAACGGCATTTCGCCATGGCAAGTCGtt                                                                                                                                                                                 >  1:331828/1‑239 (MQ=255)
                                                           tctacgaATTATGTTTTTTCCAAAAAGGGTTAAGGGTTGGGCGTTATTTGCTCATCCCGTTTTTAAGGAAACCTACCACGCTTTAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGATCCTGAAAATTAGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGct                                                                                                                                                                           <  2:716838/236‑1 (MQ=255)
                                                            cgacgaATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGGAGTCCCTAAAAGGGTTTGGGGGGCAATACCGGGCAATCCATCATCACCGGGCTTTCCCAATGGCAAATGGTTAggg                                                                                                                                                                            >  1:215978/1‑236 (MQ=255)
                                                                  atatggttttttCCAAAAGGGGTTAATGGTGGGGGCTTATAGGCCCATCCCTGTTTCAAGGAAACCTACCACGTTTTAGTTGCCGAAGGCGGGACATTAAAGGGAAAAGCGCCGGAAAATGATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                                    <  2:476774/230‑1 (MQ=255)
                                                                  atagtttaTTTTCCACAAGGGGTATATGTTGGGGGGTTTTGGTCTCATCCGTTTTTAGAAAAACACTACGGCGTTTTTGTTGCCAAAGCAGAAAAATGAAAGGAAAAAGACCGTAAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCAGCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                                    <  2:384585/231‑1 (MQ=255)
                                                                  ataatGGATTTTCCCAAAAGGGTTAATGGTTGGGCCTTATATGCTCACCCCTTTTTTCAGGAAACCTACGATGCTTTAGTTGCCAAAGTCGAGACATTAAAGGGAAAAGTTCCTAAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGa                                                                                                                                                                    <  1:1076625/238‑1 (MQ=255)
                                                                    aatGGATTTTCCAAAAAGGGTTAAGGGTTGGGCGCTATATGCTCATCCTTGTTTTCAGGAACCCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCcatcat                                                                                                                                                                                                             <  1:543344/197‑1 (MQ=255)
                                                                          ttttCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGGAATTAGGAGCCTAATAATGGCAATCCATCATCACCGGCATTCCGTCATGGCCAGGCGGTAGGCCCTGGGGAAAAATaaa                                                                                                                                                            >  1:118317/1‑239 (MQ=255)
                                                                           ataggagaaaGGGTTAATGGTGGGGCGTTTTATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCGAAAGTCGAGAAATTAAAGGAAAAAGATCCTGAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGa                                                                                                                                                          <  1:661832/232‑1 (MQ=255)
                                                                             tCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGACCCTAAAGGATTTGAAGAGAATATAATGGGCAATCCATCAACACCTGCATTCCTGCATTGCCAGTCGGTAGGCTCTGGGGAAAATaaaaaac                                                                                                                                                        >  1:44050/1‑235 (MQ=255)
                                                                                  aaaGGGGTAAGTGGTTGGGGACTATTGCCCCATCCCTTTTTGAAGAAAACTACGGGCGTTTTAGGTGCCGAAGTCCAGCCATTAAGGGGAAAAGACCCTGAAAATTATCGGAAAAACGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTc                                                                                                                                                   <  2:659671/239‑1 (MQ=255)
                                                                                      ggTTAATGGTTGGGCGTTATATGCTCATCCCTTTTTACAGGAAACCTCCGACGTTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGt                                                                                                                                                                                                                        <  2:193355/168‑1 (MQ=255)
                                                                                                 tGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGGGATTGAGGAGCATATCACGGTCAATCcatca                                                                                                                                                                                                              >  1:193355/1‑168 (MQ=255)
                                                                                                 tGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGGAGTCCATAAAGAAATTTAGGAGCAAAACACAGGCAAACCATCATCACCCGCATTCCCTTCAGGCAAGGCGTTTGGGTCCGGGGAAAATAAAAACAGGTTCACAGGAAAAATTTg                                                                                                                                     >  1:23267/1‑239 (MQ=255)
                                                                                                  gggCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAGAGATTGAGGAGCATATTACGGTTAATCCATCATCACCGGGATTCCGTCATTGCAAGACTTTAGGGTCTGGGAAAAATAAAGACTGGGCACCGGGAAAATTTTGt                                                                                                                                   >  1:105750/1‑239 (MQ=255)
                                                                                                  gggCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAATGATTGAGGAGCATATAACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATTAAAACTGGGCACGGGGAAAATTTggg                                                                                                                                   >  1:217870/1‑238 (MQ=255)
                                                                                                        tataTGCTCATCCCTGTTTTCAGGAACCCTACGACGTTTTAGTTGCCGAAGTCGAAACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGa                                                                                                                                                          <  1:525664/211‑1 (MQ=255)
                                                                                                        tataTGCTCACCCCTGTTTTGAGGAAACCTCCGCCGTTTTAGTGGCGAAAGTCGAGAAATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTgg                                                                                                                              <  1:930374/238‑1 (MQ=255)
                                                                                                               tcacctcttttttCAGGAAACTACCGCAGTTTTAGTTCCAGAAGACAAGACATTAAAGGAAAAAGATCCAGAAAATTACAATACAAACGCCGCCCCAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATc                                                                                                                     <  2:707822/230‑1 (MQ=255)
                                                                                                               tCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAATGATTGAGGGGCACATCACGGTCAATCCCTCCTCACCCGCATTCCGCCAAGGCAAGTCGTTAAGACCCTGGGAAAATAAAGACTGGGCACGGGGAAAAATTGGGGTCGGTCGGTATc                                                                                                                     >  1:274357/1‑239 (MQ=255)
                                                                                                                  tCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGAGATTGAGGAGCATATCCCGGTCCATCCATCATCACCGGCATTCCGTCATGGCAAGGAGGTAGGCTCTGGGTAAAATAAAGACTGGTCCCGGGGAAAATTTGTTGCTGGTCGttttat                                                                                                                    >  1:785717/1‑233 (MQ=255)
                                                                                                                      tGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAAGGATTTAGGGGAAAATCACGGTCAAACCATCATCACCGGCCTTCCGTCATGGCAAAACGGTAGGCTCTGGGAAAAAAAAAAACAGGGCACGGGGAAAATTTGGGGCTGGGGGTTATCGTCcctt                                                                                                              >  1:947182/1‑239 (MQ=255)
                                                                                                                       gTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAAGGATTTAGGAGACTATCACGGGCAATCCCTAATCACCGGCACTTCGTCATTGCAAATGGGTAGGCGCTGGGGAAAAAAAAAACTGGTGACGGGGAAAAATTGGTGTTGGTCGtttttcgctcttt                                                                                                             >  2:971771/1‑228 (MQ=255)
                                                                                                                             agaaaaCCTCCGACGCTTTAGTTGCCGAAGTCGAGAAATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGt                                                                                                       <  1:1012267/236‑1 (MQ=255)
                                                                                                                                 aaCCTACGACGCTTTAGTTGCCGAAGGCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTta                                                                                                     <  1:389903/237‑1 (MQ=255)
                                                                                                                                              tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAGTGATTTTGGGGCAAATCACGGGCAATCCATCATTACCCGCCTTCCGTCATGGCACGTCGTTGGGGTCTGGGGAAAAAAACAAATAGGCACGGGGAAAATTTTGGGCTGGTGGGTTATGTCTCTTCTTTTTTTATAGAGaaaaaagaa                                                                                      >  2:1005752/1‑235 (MQ=255)
                                                                                                                                              tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCCTAAAAGGGTTGGGGGGCATATACCGGTCCATCCATCATCACCGGCAATCCCTCATGGCAAGACGTTAGACTCTGGGGAAAATAAAAAATGGATCCGGGGAAAATTTTGTGGTGGGCGTTATTTTCTCTTTTTTTGCTATGAGTaaaaaaaaa                                                                                       >  2:370956/1‑236 (MQ=255)
                                                                                                                                              tttAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGAGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCCAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGGCACGGGGAAAATTTGGGGCTGGTCGTTATCGTCCCTTTTTTCGTTATAGTTaaaaaaaaa                                                                                      >  1:228469/1‑235 (MQ=255)
                                                                                                                                                tAGTTGCAGAAGTCGAGACATAAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGaaa                                                                                    <  1:309135/239‑1 (MQ=255)
                                                                                                                                                   gtGCCAAAGCCGAAACTGTAAAGGGAAAAGAACCTAAAAATTACCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGtca                                                                                <  1:289748/238‑1 (MQ=255)
                                                                                                                                                     gCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAATGATTGAGGAGCATATCACGGTCAAACCATCATCACCGGCATTCCGTCATGGCAAAGCGGTAGGCTCTGGGAAAAATAAAGACTTGTCCCGGGGAAAATTTTGGGCTGGTGGTGTTCGTCTCTTCTTTTTTTTTAGGGAAAAAAAGaaaagcat                                                                               >  1:1074119/1‑234 (MQ=255)
                                                                                                                                                         aagtggaaaatttAAAGGGAAAAGTCTCTGAAATTTATCAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCAtt                                                                           <  2:215978/228‑1 (MQ=255)
                                                                                                                                                             cgggaaaTTAAAGGAAAAAGGTCCAAAAAATTAAAGAAAAAAAGCCGCCACAAAGTTATTGGGGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTgg                                                                       <  2:947182/233‑1 (MQ=255)
                                                                                                                                                                aataTTAAAGGGAAAATATCCTAAAAATTTTAAGAGAAAGGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGgga                                                                     <  2:32470/235‑1 (MQ=255)
                                                                                                                                                                      aaaGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCAAGGTCAATCCATCATCACCAGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGGCACGGGGAAAATTTGGTGGTGGTGGTTATTGGCTCTTTTTTTGTTATAGGGAAAAAAAGAAAATAATAATTATggggaggggga                                                              >  2:298402/1‑231 (MQ=255)
                                                                                                                                                                        aGGTAAAAGTACCTGAAAATTCTCAAAAAAATGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAAc                                                            <  1:232394/239‑1 (MQ=255)
                                                                                                                                                                        aGGAAAAATATCCAAAAAATTAACAGAAAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAAc                                                            <  2:1074119/239‑1 (MQ=255)
                                                                                                                                                                        aGGAAAAAGATCCTGAAAATTATAAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAAc                                                            <  1:57618/239‑1 (MQ=255)
                                                                                                                                                                        aGGAAAAAGATCCAAAAAATTAACAGAAAAAAGCCGCCACAAAGTTTATGGCGGTTGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAAc                                                            <  2:118317/239‑1 (MQ=255)
                                                                                                                                                                                 atctaaaaatTTTTCAGAGAAACGCCCCCACAAAGTTATTGGCGGTAGTCAAAAAAGTTATGGGGGAGCATCTTACGGTCAAGCCATCAACGCCGGCATTCCGTCATGGCAAGTCGTTAGGCGCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAac                                                   <  2:815187/233‑1 (MQ=255)
                                                                                                                                                                                 aTCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCAAATCACGGGCAATCCATCATCACCGGGATTCCGTCATGGGAAGGCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGGGCTGGGCGTTATCGGCTCTTCTTTCGTTATATTGAAAAAGAGAAAATCATCATTCTTGGATGGATCAACGATGaaaaa                                                   >  1:114907/1‑238 (MQ=255)
                                                                                                                                                                                    cTGAAAATTAACAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGGCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACAct                                                <  1:720264/239‑1 (MQ=255)
                                                                                                                                                                                       aaaaTCAGCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATACCACGGTCAAGCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTg                                             <  1:53758/239‑1 (MQ=255)
                                                                                                                                                                                             attggaaaaaaGCCGCCTCAAAGTTATGGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCCCGGTACATCCAGCATCACCGGCACTCCGTCATGGCAAGGCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTGCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCAcc                                       <  2:752446/232‑1 (MQ=255)
                                                                                                                                                                                                agaaAAAAGCCGCCAAAAAGTTATTTGCGGTATTCCATAAATGGATTGAAGCGCATATCCCGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGGCTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTAc                                    <  1:425076/239‑1 (MQ=255)
                                                                                                                                                                                                   aaaaaCCCCCCAAAAAGTTATTGGCGGGAGTCCATAAAGTGATTTAGGAGCATAACCCGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTAc                                    <  2:897907/235‑1 (MQ=255)
                                                                                                                                                                                                    aaaaGCCCCCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTAc                                    <  2:1044102/236‑1 (MQ=255)
                                                                                                                                                                                                      aaGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAATGATTGAGGGGCATATCACGGTCAATCCATCATCACCGGCCTTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAAAAAAAACTGGTCACGGGGAAAATTTGGTGCTGGGCGGTATCGTCTCTTCTTTCGTTATAGGGAAAAAGAGAAAATCATCATTCTGGGGTGGATGAAACATGAAAAAACCTTGCGCCCCTCCGGTaaa                              >  1:1003045/1‑239 (MQ=255)
                                                                                                                                                                                                           cggccaaaaGGTTTTGCCGGGTGTCCATAAAGGTAATGAGGGGCTAATCACGGGCAAGCCCTCATTCACGGCATTCCGTCATGGGAAGTCCTTAGGCTCTGGGAAAAATAAAGACTGGTCACCGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaa                            <  1:425868/231‑1 (MQ=255)
                                                                                                                                                                                                             gcacaAAGTTTTTGGGGGTAGTCCAAAAAGTTATTTGGGAGCATTGCACGGTCAAGCCATCATCACCGGCATTCCGTCATGGCAAGGCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                        <  2:592831/238‑1 (MQ=255)
                                                                                                                                                                                                             cctcaAATTTTTTGGCGGTATTCCAAAAAGTACTTGAGGAACATTACAGGGTACATCCAACATCAACGGCATACCGGCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                        <  2:705727/236‑1 (MQ=255)
                                                                                                                                                                                                             ccccaAAGTTATTGGCGGTAGTCCATAAAGTGAATGAGGGGCATATCACGGTCAATCCATCAGCACCGGCATTCCGGCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                        <  2:873103/236‑1 (MQ=255)
                                                                                                                                                                                                             ccACAAAGTTATTGGCGGTAGTGCATAAAAGAGTTTAGGGGGATATCACGGTCAATCCCTCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGGTCTGGGGAAAAAAAAGAATGGTCACGGGGAAAATTTTGTGCTGGTCGTTATCGTCTCTTCTTTTGTTATAGTGAAAAAAAGAAAAGAATCATTTTGTGGTGGATGAACGATAGAAAAACTCTTGGCCCCTAAGGGgaaaaaaaa                        >  2:1081596/1‑237 (MQ=255)
                                                                                                                                                                                                             ccACAAAGTTATTGGCGGTAGTCCATAAAATGATTGAGGGGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCCAATCGTTAGGCTCTGGGGAAAATAAAAACTGGGCACCGGGAAAAATTGGTGCTGGGAGTTTACGTTCCTTTTTTTGTTATAGTGAAAAAGAGAAAATCATCATTTTGGGATGGGTGAACGTTGAAAACAATCTGCGGACCCACGGTaaaaaaaac                        >  1:222981/1‑237 (MQ=255)
                                                                                                                                                                                                              cacaAAGTTAATGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCAACATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACa                        <  1:38279/238‑1 (MQ=255)
                                                                                                                                                                                                                       tttGGCGGTTGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATGTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCtata              <  2:429885/238‑1 (MQ=255)
                                                                                                                                                                                                                       aTTGGCGTTAGGCCATAAAGGGATTGAGGGGCATATACCGGTCAATCCATCATCACCGGACTTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGCCGGGTCACGGGTAAAAGTGGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCtata              <  2:249127/239‑1 (MQ=255)
                                                                                                                                                                                                                       aTTGGCGGTAGTCCATAAAGGGATTGAGGAGCATTGCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCtata              <  1:319912/239‑1 (MQ=255)
                                                                                                                                                                                                                        ttGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAAATTGGGGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAAAGAAAGTCATCATTTTGGGATGGATGAACGATGAAAAAACTCTTCGCACCTACGGTAAAAAAACAAATGCCtcta              >  1:109962/1‑238 (MQ=255)
                                                                                                                                                                                                                             ggTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTa         <  1:370956/238‑1 (MQ=255)
                                                                                                                                                                                                                                   ccATAAAGGGATGGAGGAGGATTTCACGGGTAAGCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACGGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCa  <  2:952034/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
CAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCA  >  NC_000913/3277269‑3277730

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