Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,406,735	T→C	G273G (GGT→GGC)	accC →	biotin carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,406,735	0	T	C	100.0%	116.6 / NA	34	G273G (GGT→GGC)	accC	biotin carboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/19);  total (15/19)

GCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGTACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAA  >  NC_000913/3406605‑3406849                                                                                                                                   |                                                                                                                   gCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCAc                                                                                                                  <  2:292825/133‑1 (MQ=255)   tCCATGCATCGCCGCGACTAGAGAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACtt                                                                                                                <  1:215095/133‑1 (MQ=255)       tGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAACGCGTGTGTTGATATCGCCTATCGCGGGGCATGCACTTTCGa                                                                                                            >  2:215095/1‑133 (MQ=255)               cgccACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGt                                                                                                    >  2:329166/1‑133 (MQ=255)               cgccACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGt                                                                                                    >  1:353743/1‑133 (MQ=255)                              gTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCCCTTTCGAATTCCTGATCGAAACCGGCGAGt                                                                                     >  1:222608/1‑133 (MQ=255)                              gTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGt                                                                                     >  2:51736/1‑133 (MQ=255)                              gTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGt                                                                                     >  1:311718/1‑133 (MQ=255)                              gTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGt                                                                                     >  1:360861/1‑133 (MQ=255)                               tCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGtt                                                                                    >  1:302994/1‑133 (MQ=255)                                  aagaagCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTa                                                                                 >  1:85931/1‑133 (MQ=255)                                          ccAGCACCCGGCATTCCCCCGGCACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCg                                                                         <  1:162546/133‑1 (MQ=255)                                                     cATTACCCCGTAACTCCGTCGCTACATCAGCGCACGTTGCGGTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACAcc                                                              <  2:160035/133‑1 (MQ=255)                                                     cATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACAcc                                                              <  1:516093/133‑1 (MQ=255)                                                          ccccGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAACGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTAt                                                         <  1:374622/133‑1 (MQ=255)                                                               gAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTTCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAgg                                                    >  2:318861/1‑133 (MQ=255)                                                                     cGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAac                                              >  2:318815/1‑133 (MQ=255)                                                                          cTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCCATTTCATCCAAATGAAAACCCGTATTCAGGTAGAACACCCg                                         >  2:227661/1‑133 (MQ=255)                                                                            aCATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGt                                       <  1:51736/133‑1 (MQ=255)                                                                              aTCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTa                                     <  2:353743/133‑1 (MQ=255)                                                                                  gCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATAGAAATGAAAACCGGTAGTCAGGTAGAACACCCGGTTACAGa                                 >  2:140659/1‑133 (MQ=255)                                                                                     aaCGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGaa                                <  1:277165/131‑1 (MQ=255)                                                                                     aaCGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGaa                                >  2:277165/1‑131 (MQ=255)                                                                                      aCGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATg                             <  2:384268/133‑1 (MQ=255)                                                                                        gTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGAt                           <  2:302994/133‑1 (MQ=255)                                                                                          tGCGCTAACGCGTGTATTCATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCa                         <  1:140659/133‑1 (MQ=255)                                                                                                    cGTGTGTTGATATCGGCTATCGCGGTTCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGa               <  1:318861/133‑1 (MQ=255)                                                                                                     gtgtgtTGATATCGGTTATCGCGGAGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGAc              <  2:222608/133‑1 (MQ=255)                                                                                                     gtgtgtTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGAc              <  1:318815/133‑1 (MQ=255)                                                                                                      tgtgtTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGAcc             >  2:68440/1‑133 (MQ=255)                                                                                                       gtgtTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCt            <  1:307445/133‑1 (MQ=255)                                                                                                         gtTGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGa          >  2:463041/1‑133 (MQ=255)                                                                                                          ttGATATCGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGAt         <  2:311718/133‑1 (MQ=255)                                                                                                                 cgacTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGaa  <  2:360861/130‑1 (MQ=255)                                                                                                                 cGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGaa  <  1:329166/133‑1 (MQ=255)                                                                                                                 cGGCTATCGCGGTGCAGGCACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGaa  <  2:85931/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                   GCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGCTATCGCGGTGCAGGTACTTTCGAGTTCCTGTTCGAAAACGGCGAGTTCTATTTCATCGAAATGAACACCCGTATTCAGGTAGAACACCCGGTTACAGAAATGATCACCGGCGTTGACCTGATCAAAGAA  >  NC_000913/3406605‑3406849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑