Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,277,509 A→G 16.0% 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,5090AG16.0% 166.8 / 5.5 76E58G (GAG→GGG) yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (17/45);  new base G (9/4);  total (26/50)
Fisher's exact test for biased strand distribution p-value = 8.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAA  >  NC_000913/3277273‑3277740
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
aaGGAAAGAACCTTTCCGCTGGCAGGGACCTAAAAATTGGTGGAGAGATTGGCGACGAAGAATAATGGATTTTCCACAAAGGGTTAAGGGTTGGGCGATATATGCTCACCCCTGTTTCCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATAAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGc                                                                                                                                                                                                                                            <  2:1015247/239‑1 (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)
                                                                                                                                                                                                                                         taTTGGGGGGCAATTCACGGGTAATCCAATATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGtt     <  2:228469/236‑1 (MQ=255)
                                                                                                                                                                                                                                         gATTGAGGAGCATAGCACGGTCAACCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCGGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCTCCTACGGTAAAAAAACa                                  <  1:107221/208‑1 (MQ=255)
                                                                                                                                                                                                                                          aTTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGGGCACCTACGGTAAAAAAACAGATGCCTATAACGAATTCAGCAAAATGTTaaa  >  1:956924/1‑239 (MQ=255)
                                                                                                                                                                                                                                          aTTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa                                >  2:6653/1‑209 (MQ=255)
                                                                                                                                                                                                                                          aTTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGa                                <  1:6653/209‑1 (MQ=255)
                                                                                                                                                                                                                                          aTTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaaa  >  1:332292/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
AAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATGAGATTGGCGACGACGAATAATGGATTTTCCACAAAGGGTTAATGGTTGGGCGCTATATGCTCATCCCTGTTTTCAGGAAACCTACGACGCTTTAGTTGCCGAAGTCGAGACATTAAAGGGAAAAGATCCTGAAAATTATCAGAGAAAAGCCGCCACAAAGTTATTGGCGGTAGTCCATAAAGTGATTGAGGAGCATATCACGGTCAATCCATCATCACCGGCATTCCGTCATGGCAAGTCGTTAGGCTCTGGGAAAAATAAAGACTGGTCACGGGTAAAATTTGGTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAA  >  NC_000913/3277273‑3277740

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: