Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,899,745 A→G 7.3% L278P (CTT→CCT)  yieL ← putative xylanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,899,7450AG 7.3% 175.7 / 3.5 55L278P (CTT→CCT) yieLputative xylanase
Reads supporting (aligned to +/- strand):  ref base A (30/20);  new base G (2/2);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAACGTCCGGTATCACC  >  NC_000913/3899471‑3900028
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                           
cAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAgcgcg                                                                                                                                                                                                                                                                               <  2:104768/289‑1 (MQ=255)
        gATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATc                                                                                                                                                                                                                                                                       >  1:286015/1‑289 (MQ=255)
                    aaTCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACg                                                                                                                                                                                                                                                           <  2:286015/289‑1 (MQ=255)
                                            cAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACGAGAGGAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                                                                                                                                                                                                                                   <  1:269613/289‑1 (MQ=255)
                                            cAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                                                                                                                                                                                                                                   <  1:190977/289‑1 (MQ=255)
                                                               ggCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTtcatc                                                                                                                                                                                                                >  1:118390/1‑289 (MQ=255)
                                                                                ttACGACATCTTTATCCCCCACAACCACCGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCg                                                                                                                                                                                                <  2:118390/288‑1 (MQ=255)
                                                                                         cTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCAt                                                                                                                                                                                       >  1:119910/1‑288 (MQ=255)
                                                                                                   cacaACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCgg                                                                                                                                                                                                                                                                     >  1:414207/1‑200 (MQ=255)
                                                                                                   cacaACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCgg                                                                                                                                                                                                                                                                     <  2:414207/200‑1 (MQ=255)
                                                                                                   cacaACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCCGTTCGCGATCCGCCGCTTTAGCATTCAGCGGATaa                                                                                                                                                                            >  2:430430/1‑289 (MQ=255)
                                                                                                   cacaACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGAGGCCATCGGCATCTTTACGGACCTTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGTACTGCCGCTTTAGCATTCACCGGTTaa                                                                                                                                                                            >  2:170061/1‑289 (MQ=255)
                                                                                                   cacaACCACAGCAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGCGGTAACACCACTGAATGTGGCCAGCAAGCCAAAGCCTTCCAGATGATTCATTCCGGAAACCCGCGCCAGGTACCCGCCCTGCGCAAGCCCTGCCAGCGCGCTGCCTTCGGCATCTTTACGCACATAAAAAACCATGCGCAACCGAGGGTTAATACCGTCCATCCGTTCACGATCGGCCCCTGTAGcctctgacggtaaa                                                                                                                                                                            >  2:309061/1‑276 (MQ=255)
                                                                                                    acaACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATaaa                                                                                                                                                                           >  1:29425/1‑289 (MQ=255)
                                                                                                                      cGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCttt                                                                                                                                                            >  1:137695/1‑286 (MQ=255)
                                                                                                                        tAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCCTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGCTACAAGACTTTACGTTTTTCCTGa                                                                                                                                                       >  2:78078/1‑289 (MQ=255)
                                                                                                                             gCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAGGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCa                                                                                                                                                                                    >  2:132797/1‑255 (MQ=255)
                                                                                                                             gCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAGGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCa                                                                                                                                                                                    <  1:132797/255‑1 (MQ=255)
                                                                                                                             gCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                                                                                                                                                                                                                                   >  1:483757/1‑208 (MQ=255)
                                                                                                                             gCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                                                                                                                                                                                                                                   <  2:483757/208‑1 (MQ=255)
                                                                                                                               tGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACTGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACg                                                                                                                                                 <  1:78078/288‑1 (MQ=255)
                                                                                                                                     aTAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCCTCAGTTCGTGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCt                                                                                                                                           >  1:283011/1‑288 (MQ=255)
                                                                                                                                     aTAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCtt                                                                                                                                          >  1:293801/1‑289 (MQ=255)
                                                                                                                                                            ggCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAACGCTTTCCAGATGATTCCTTCCGGAAACCATCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCACCGCGCGGCCATCGGCCTCTTTCCGGACATTAAAACGCTTGCTAATCAGCCGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCAGATAAAAGACTTTACGTCTTTCCTGAGTCTCGAACTCTTCGGGACTAATGCCCTTCGCATCg                                                                                                                   >  2:34940/1‑289 (MQ=255)
                                                                                                                                                                gagACACCTACATCCAGTACGGTCGTCGTGCTAACACCACTGAATGTGGCCATCCAGACAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAATCGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGttt                                                                                                               <  1:309061/287‑1 (MQ=255)
                                                                                                                                                                 cgACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAGGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCTCTTCGCATCGGtt                                                                                                                >  2:35298/1‑287 (MQ=255)
                                                                                                                                                                                 tACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCg                                                                                                      <  2:271739/281‑1 (MQ=255)
                                                                                                                                                                                 tACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCg                                                                                                      >  1:271739/1‑281 (MQ=255)
                                                                                                                                                                                 tACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCGGTACCCGGGATccc                                                                                               >  2:240106/1‑286 (MQ=255)
                                                                                                                                                                                     gTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCaccacc                                                                                           >  2:66080/1‑288 (MQ=255)
                                                                                                                                                                                           gTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                                                                                                                                                                                                                                   <  1:352580/146‑1 (MQ=255)
                                                                                                                                                                                           gTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAg                                                                                                                                                                                                                                   >  2:352580/1‑146 (MQ=255)
                                                                                                                                                                                                   accacTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCCTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGTACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAAtt                                                                            >  2:30505/1‑289 (MQ=255)
                                                                                                                                                                                                   accacTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCCTCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAACGACTTTACGTTTTTCCTCAGGCACGAAATCTTCGGGAATAATGCCCTTAGCAACGGTTTCGTTATCCGGGATCACCACCAGCATCTGTTTAAtt                                                                            >  2:459838/1‑289 (MQ=255)
                                                                                                                                                                                                        tGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTc                                                                         >  1:274249/1‑287 (MQ=255)
                                                                                                                                                                                                          aaTGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCt                                                                      >  1:169482/1‑288 (MQ=255)
                                                                                                                                                                                                           aTGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAGGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCTCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCtt                                                                     <  1:35298/288‑1 (MQ=255)
                                                                                                                                                                                                            tgtgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCAAGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGAGATCACCACCAGCATCGGTTTAATTTTCCCttaa                                                                   >  1:483952/1‑287 (MQ=255)
                                                                                                                                                                                                             gtgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGc                                                                                                                                                                                  >  2:126197/1‑177 (MQ=255)
                                                                                                                                                                                                             gtgGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGc                                                                                                                                                                                  <  1:126197/177‑1 (MQ=255)
                                                                                                                                                                                                                      cagccaAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAgg                                                          <  1:444581/288‑1 (MQ=255)
                                                                                                                                                                                                                           aaaGCTTGCCAGAGGATTCAGTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATcc                                                    <  2:137695/289‑1 (MQ=255)
                                                                                                                                                                                                                                        tggttcattcCGGAAACCATGGCCTGGTACCCGCCTTGTGAACGCCCTGCCAGCGCGCTGCCATCGGCATCTTTACTGACATTAAAACGCTTGCTAGTCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGTATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGa                                       <  2:283011/286‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaaCCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAgg                                                          >  2:322561/1‑257 (MQ=255)
                                                                                                                                                                                                                                                     aaaCCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAgg                                                          <  1:322561/257‑1 (MQ=255)
                                                                                                                                                                                                                                                            cgcCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCg                   >  1:87279/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               cTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGa                 <  2:293801/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGa                                                                                                                                                                                                                                                         >  2:289050/1‑53 (MQ=255)
                                                                                                                                                                                                                                                                  gTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGa                                                                                                                                                                                                                                                         <  1:289050/53‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAACGTCCGGt        <  2:29425/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAACGTCCGGt        <  2:119910/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCgg                                                                                                                                                                                                                                >  2:19503/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                        gCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCgg                                                                                                                                                                                                                                <  1:19503/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCGTCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTCCCTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCCGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCGGATCGCTAGCGGAACGTCCGGTATCAc   >  2:143977/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             gtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAACGTCCGGTATCAcc  >  2:8186/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                           
CAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGGTTGATAGCTGCCGGATCGTTCAGTCGGGCCGCGACACCTTCATCCGGTACGGTTGTCGTGGTAACACCACTGAATGTGGCCAGCCAGCCAAAGCTTTCCAGATGATTCATTCCGGAAACCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAACGTCCGGTATCACC  >  NC_000913/3899471‑3900028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: