Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,221,860 T→G K44Q (AAA→CAA)  bglX ← beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,221,8600TG97.3% 248.3 / NA 74K44Q (AAA→CAA) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (39/33);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAG  >  NC_000913/2221631‑2222087
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
aaCACCGTGCGCTGACCGTTCAGCAAGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATGCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAg                                                                                                                                                                                                                               <  2:428828/236‑1 (MQ=255)
aaCACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAg                                                                                                                                                                                                                               <  1:158948/236‑1 (MQ=255)
     cGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGtt                                                                                                                                                                                                                          >  2:100679/1‑236 (MQ=255)
       tGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCg                                                                                                                                                                                                                        >  2:246229/1‑236 (MQ=255)
       tGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCg                                                                                                                                                                                                                        >  1:252661/1‑236 (MQ=255)
         cgcTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGt                                                                                                                                                                                                                      >  1:612457/1‑236 (MQ=255)
                          gtcgtAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCcgc                                                                                                                                                                                                      >  1:179597/1‑235 (MQ=255)
                            cgtAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCt                                                                                                                                                                                                   <  1:417789/236‑1 (MQ=255)
                             gtAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCtt                                                                                                                                                                                                  <  2:381343/236‑1 (MQ=255)
                                  cAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGc                                                                                                                                                                                             <  1:246229/236‑1 (MQ=255)
                                     agaaaagaGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGtt                                                                                                                                                                                          >  2:590775/1‑236 (MQ=255)
                                     agaaaagaGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGtt                                                                                                                                                                                          >  1:291475/1‑236 (MQ=255)
                                      gaaaagaGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTa                                                                                                                                                                                         >  1:425825/1‑236 (MQ=255)
                                                aTTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatggat                                                                                                                                                                                   <  1:368019/232‑1 (MQ=255)
                                                 ttttCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTg                                                                                                                                                                              >  1:40520/1‑236 (MQ=255)
                                                 ttttCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATCTGCTTATGCAGTTCGGTGACGAACGCATTCCGCGCTTCGGGCGTTAATGGATGGTTg                                                                                                                                                                              >  2:558966/1‑236 (MQ=255)
                                                     cAGGCGGCTGAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGa                                                                                                                                                                          <  2:432697/236‑1 (MQ=255)
                                                                         ccTGATCCTGCATGGCGCGGATATCCTGGCGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgct                                                                                                                                                      <  1:763282/236‑1 (MQ=255)
                                                                              tccgtcATGGCCGGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggctg                                                                                                                                                 <  1:105625/231‑1 (MQ=255)
                                                                              tCCTGCATGGCGCGGATTTCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTTACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggctg                                                                                                                                                 <  2:189336/236‑1 (MQ=255)
                                                                              tCCTGCATGGCCGGGATTGCCTGACGGGTTCCGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggctg                                                                                                                                                 <  2:40520/236‑1 (MQ=255)
                                                                               ccTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                                                                                                                                                              >  2:760329/1‑190 (MQ=255)
                                                                               ccTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                                                                                                                                                              <  1:760329/190‑1 (MQ=255)
                                                                                cTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCa                                                                                                                                               <  1:64829/236‑1 (MQ=255)
                                                                                 tGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCa                                                                                                                                               <  1:100679/235‑1 (MQ=255)
                                                                                        cgcgGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAg                                                                                                                                        >  2:206761/1‑235 (MQ=255)
                                                                                        cgcgGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAg                                                                                                                                        >  2:617734/1‑235 (MQ=255)
                                                                                        cgcgGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGa                                                                                                                                       >  1:652278/1‑236 (MQ=255)
                                                                                                 ccTGACGGGTTACGGTGTTGAAAATAGCCCCAACCTGCCCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGCCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGGGACGAACGACTCCCGCGCTTCGGGCGTCAATGGATGGGTGCCGAATAAATCATCCTCCAGGACTGGCTGAAGGGGCCGACTctgcgc                                                                                                                               >  2:505169/1‑230 (MQ=255)
                                                                                                     aCGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGAtt                                                                                                                           >  1:420522/1‑235 (MQ=255)
                                                                                                             cGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTa                                                                                                                       <  1:631807/231‑1 (MQ=255)
                                                                                                                  ttGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACCTAg                                                                                                             >  2:217717/1‑236 (MQ=255)
                                                                                                                   tGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCCGGGCCAGACTTACCTCGATTCCTACTGAACATAGc                                                                                                            >  2:602202/1‑236 (MQ=255)
                                                                                                                   tGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGc                                                                                                            >  1:479324/1‑236 (MQ=255)
                                                                                                                     acaaTCGCCCCAACCTGACCGTCTTTGCTCATCTCGCGGTTCGCCTCTTTCGGGTTATCCGGGCCTACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGcc                                                                                                           <  2:417811/233‑1 (MQ=255)
                                                                                                                      aattCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCa                                                                                                          <  1:409569/232‑1 (MQ=255)
                                                                                                                            ccccAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCg                                                                                                                                                                                                                        <  1:143150/119‑1 (MQ=255)
                                                                                                                            ccccAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCg                                                                                                                                                                                                                        >  2:143150/1‑119 (MQ=255)
                                                                                                                            ccccAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGAt                                                                                                                            <  1:110792/211‑1 (MQ=255)
                                                                                                                            ccccAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGAt                                                                                                                            >  2:110792/1‑211 (MQ=255)
                                                                                                                            ccccAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCata                                                                                                   >  1:673545/1‑236 (MQ=255)
                                                                                                                                      aCCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGAt                                                                                                                            >  2:41297/1‑201 (MQ=255)
                                                                                                                                      aCCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGAt                                                                                                                            <  1:41297/201‑1 (MQ=255)
                                                                                                                                                    tctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGccc                                                                           >  1:339286/1‑236 (MQ=255)
                                                                                                                                                     ctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCg                                                                          >  1:570021/1‑236 (MQ=255)
                                                                                                                                                                 ctctTTCGGGTTATCCGGGCCTACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGa                                                              >  1:626965/1‑236 (MQ=255)
                                                                                                                                                                  tctTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGAGTTTTCTTCCGTTGTTGCCCGACGTGGCAGCgta                                                             >  2:89580/1‑234 (MQ=255)
                                                                                                                                                                              tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGCGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCaaaa                                                 >  1:16377/1‑236 (MQ=255)
                                                                                                                                                                              gccGGGCCGACGCTGATTAAGCGCAGCTGACCAATTGTCTCATAAGCTCTCATTTTCTTAAGCAGTTCGGTAACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATTGATTTTTCTTCCTTGGTTGCCCGACGTGGCAGCGAAAATGGTGCaaa                                                  <  2:41007/234‑1 (MQ=255)
                                                                                                                                                                                   gCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGt                                            >  1:518187/1‑236 (MQ=255)
                                                                                                                                                                                   gCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGt                                            >  2:461067/1‑236 (MQ=255)
                                                                                                                                                                                       aCGCTGATTAAGCGCAGCTGCCCAATTTCCTCCTCAACTGTCATTTGCTTAGGCAGTTCGGTGACGAACGCATCCCGCGTTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGt                                         <  2:252661/235‑1 (MQ=255)
                                                                                                                                                                                       aCGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGt                                         <  2:632703/235‑1 (MQ=255)
                                                                                                                                                                                       aCGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGt                                         <  2:179597/235‑1 (MQ=255)
                                                                                                                                                                                           ttaTTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGt                                         <  2:612457/229‑1 (MQ=255)
                                                                                                                                                                                                     cAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCatga                          >  1:643861/1‑236 (MQ=255)
                                                                                                                                                                                                       gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAAt                                                          <  2:184305/202‑1 (MQ=255)
                                                                                                                                                                                                       gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAAt                                                          >  1:184305/1‑202 (MQ=255)
                                                                                                                                                                                                       gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCatgatg                        >  1:239304/1‑236 (MQ=255)
                                                                                                                                                                                                       gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCatgatg                        >  2:322836/1‑236 (MQ=255)
                                                                                                                                                                                                       gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCatgatg                        >  2:73921/1‑236 (MQ=255)
                                                                                                                                                                                                           aCCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTg                                                                                                                    >  1:648049/1‑140 (MQ=255)
                                                                                                                                                                                                           aCCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTg                                                                                                                    <  2:648049/140‑1 (MQ=255)
                                                                                                                                                                                                            ccATTTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGga                      <  1:313895/233‑1 (MQ=255)
                                                                                                                                                                                                            ccAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGgaga                    <  1:89603/235‑1 (MQ=255)
                                                                                                                                                                                                                ttttCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGaa               <  2:410525/236‑1 (MQ=255)
                                                                                                                                                                                                                ttttCCAATCAACTGTAATTGTGTTAAGCGTTTACGGAACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGaa               <  2:749182/236‑1 (MQ=255)
                                                                                                                                                                                                                     tcatcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGGGacaagg           >  1:130856/1‑233 (MQ=255)
                                                                                                                                                                                                                       atcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGaa         >  1:179767/1‑235 (MQ=255)
                                                                                                                                                                                                                         caACTGTCATTTGCATAAGCAGTTCGGTGACGAACGCATCCCTCGCTTCGGGAGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAGAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGc      <  2:705757/236‑1 (MQ=255)
                                                                                                                                                                                                                           acgttaaTTTGCATAAGCAGGTGGGTAACGAACGCATCCCGCGTTTCGGGCGTTAAGGGATGGTTGCCGAATAATTCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTc    <  2:130856/230‑1 (MQ=255)
                                                                                                                                                                                                                           aCTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTc    <  1:461067/236‑1 (MQ=255)
                                                                                                                                                                                                                           aCTGGCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTc    <  1:590775/236‑1 (MQ=255)
                                                                                                                                                                                                                                 aTTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAg  >  1:318813/1‑232 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
AACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACTCACCGCGATTCCTACTGAACATAGCCATTTCATATGGATTTTTCTTCCTTTGTTGCCCGACGTGGCAGCGAAAATGGTGCAAAAACCGTAGTTTGCCATAAGCATGATGGAGAGAGAAAAAGAATGCTCAG  >  NC_000913/2221631‑2222087

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