Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,261,926	G→T	P147Q (CCG→CAG)	fruK ←	1‑phosphofructokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,261,926	0	G	T	100.0%	146.0 / NA	42	P147Q (CCG→CAG)	fruK	1‑phosphofructokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (27/15);  total (27/15)

TTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCGAGAAGTTGAAGTCGGT  >  NC_000913/2261789‑2262056                                                                                                                                          |                                                                                                                                   tttCCAGCTCGCGGCGGTTAGGGTTCACCAGCCACCGGGCCGCTTTCAAACCTGCTACTAAAGCTTCACCGCTACTATCAAAAATAATGCAAGGACACTGACTACGCAGGCGAGGCATCCAGTCGGGGAACGCTTCCtgg                                                                                                                                  >  2:214246/1‑140 (MQ=255)          cgcgGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGa                                                                                                                              <  1:170704/135‑1 (MQ=255)          cgcgGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGa                                                                                                                              >  2:170704/1‑135 (MQ=255)                    aGGTTTCACCAGCCCCGGGGCCGCTTTTAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGGCATCCAGTCGGGGAACGCTTCCTGGCTGACGCCTGACGGGAAGc                                                                                                               >  2:827536/1‑140 (MQ=255)                     ggTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCt                                                                                                              <  2:176368/140‑1 (MQ=255)                     ggTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCt                                                                                                              <  2:127937/140‑1 (MQ=255)                         tCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCg                                                                                                          <  1:144042/140‑1 (MQ=255)                                ccACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGacac                                                                                                   <  1:827536/140‑1 (MQ=255)                                ccACGGCGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGacac                                                                                                   >  1:104059/1‑140 (MQ=255)                                   cGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGa                                                                                                <  1:763777/140‑1 (MQ=255)                                    ggTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAc                                                                                               >  2:965512/1‑140 (MQ=255)                                           cTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCg                                                                                        >  2:137651/1‑140 (MQ=255)                                             ttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGaa                                                                                      >  1:510876/1‑140 (MQ=255)                                              tCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                                                                                     >  1:890105/1‑140 (MQ=255)                                                      gCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAc                                                                                 <  2:201643/136‑1 (MQ=255)                                                      gCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAc                                                                                 >  1:201643/1‑136 (MQ=255)                                                       ctactaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGaa                                                                                      <  1:448275/130‑1 (MQ=255)                                                       ctactaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGaa                                                                                      >  2:448275/1‑130 (MQ=255)                                                          ctaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCAGATCGAACTGACCGAGCCAg                                                                         <  2:510876/140‑1 (MQ=255)                                                                 ttCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCagag                                                                  >  2:930545/1‑140 (MQ=255)                                                                 ttCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCagag                                                                  >  1:713597/1‑140 (MQ=255)                                                                    aCGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAAt                                                               <  1:930545/140‑1 (MQ=255)                                                                         tactaTCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTc                                                          >  2:638912/1‑140 (MQ=255)                                                                           ctaTCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAg                                                                         >  2:90991/1‑123 (MQ=255)                                                                           ctaTCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTcac                                                        >  2:959441/1‑140 (MQ=255)                                                                           ctaTCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTGCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAg                                                                         <  1:90991/123‑1 (MQ=255)                                                                                  aGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGt                                                 >  1:21928/1‑140 (MQ=255)                                                                                       aTGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCa                                            <  1:203960/140‑1 (MQ=255)                                                                                            aGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTcgg                                       <  2:618460/140‑1 (MQ=255)                                                                                              gACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTcggcg                                     >  2:774544/1‑140 (MQ=255)                                                                                                  cTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCggggg                                 >  1:168610/1‑140 (MQ=255)                                                                                                   tGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGt                                >  2:953024/1‑140 (MQ=255)                                                                                                        aCGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACtt                           >  2:396879/1‑140 (MQ=255)                                                                                                        aCGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACtt                           >  1:1042905/1‑140 (MQ=255)                                                                                                             ggCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCaaaa                      >  2:804225/1‑140 (MQ=255)                                                                                                             ggCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCaaaa                      >  1:1009046/1‑140 (MQ=255)                                                                                                              gCGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAAc                     <  1:396879/140‑1 (MQ=255)                                                                                                               cGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGGCACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGGGACTTCACAAcc                    >  2:530531/1‑140 (MQ=255)                                                                                                               cGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAAcc                    >  2:64881/1‑140 (MQ=255)                                                                                                               cGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAAcc                    >  1:574287/1‑140 (MQ=255)                                                                                                               cGAGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAAcc                    >  2:961887/1‑140 (MQ=255)                                                                                                                 aGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCg                  >  2:1054815/1‑140 (MQ=255)                                                                                                                 aGTCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCg                  >  2:159772/1‑140 (MQ=255)                                                                                                                   tCATCCAGTCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCgag                >  1:617003/1‑140 (MQ=255)                                                                                                                           tCGGTGAACGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCGAGAAGTTGaa        <  1:953024/140‑1 (MQ=255)                                                                                                                                 aaCGCTTCCTGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCGAGAAGTTGAAGTCGGt  <  2:1042905/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   TTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAGAATCAGTCACAAAGCGTTCCCAGTCGGCGGGGGTGACTTCAAAACCCGAGAAGTTGAAGTCGGT  >  NC_000913/2261789‑2262056

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑