Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,892,006	A→C	intergenic (‑108/‑208)	cysJ ← / → queD	sulfite reductase, flavoprotein subunit/6‑carboxy‑5,6,7,8‑tetrahydropterin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,892,006	0	A	C	100.0%	118.3 / NA	42	intergenic (‑108/‑208)	cysJ/queD	sulfite reductase, flavoprotein subunit/6‑carboxy‑5,6,7,8‑tetrahydropterin synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/23);  total (19/23)

CAACGCGGAAGGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAA  >  NC_000913/2891869‑2892141                                                                                                                                          |                                                                                                                                        cAACGCGGAAGGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTACAGACGGGCTAACg                                                                                                                                    >  1:359690/1‑143 (MQ=255)  aaCGCGGAAGGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGCAGCGAAGGTAAGGATAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGa                                                                                                                                   >  2:115322/1‑143 (MQ=255)            gTGGGACTTGTGTCGTCATGCGTCGTTAGGTTCCAGTAAGCAAAGCTGTTTCTGTGTCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACtaataa                                                                                                                         <  1:780770/143‑1 (MQ=255)             tGGCACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCTATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAAc                                                                                                                        <  1:621954/143‑1 (MQ=255)               ggACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTACTAACCa                                                                                                                      >  2:643851/1‑143 (MQ=255)                 aCCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGACCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGTGCTAGCGATAACTAATAACCaaa                                                                                                                    >  1:743411/1‑143 (MQ=255)                     gtgtCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCTCTGTTAGTGCTCATAGAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGa                                                                                                                <  1:460034/143‑1 (MQ=255)                          gtcATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAAc                                                                                                           <  2:45343/143‑1 (MQ=255)                            cATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCCTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCATACGGTGTAGATTAAAGAAAGGCTAGCGATAACAAATAAACAAATCGACTAACCt                                                                                                         >  1:274186/1‑143 (MQ=255)                                          ttCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGAGCAAACAGTAAAAATAAAAAACAGGCAAGCAATAACTAATAACCAAATCAACTAACCTGTTTTAGCAATAGt                                                                                           >  1:624816/1‑143 (MQ=255)                                           tCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGTCCATAAATCAGAAGAGAAGTTAAGGTTAACGGGGAAAAAGGTGTGGAGTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTc                                                                                          <  1:630042/143‑1 (MQ=255)                                                taattaaaGCTGTTTCTGCGCCCTGTCAGGGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTaa                                                                                     <  2:786054/138‑1 (MQ=255)                                                   gCAGAGCTGTTGCTGCGCCCGGGCAGCGGCCATAAAACAGAAGAGAAGGTAAGGGTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTaacaa                                                                                  <  2:519557/143‑1 (MQ=255)                                                            tttCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTATTAGCAATAGTCTTTAACAACAAAATAGa                                                                         >  1:750279/1‑143 (MQ=255)                                                             ttCTGTGTCCTGTCTTCGCCCATAATATAGAAGAGAAGGTAAGATAAACGGGGGAAACGGTGGGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGAt                                                                        <  1:333181/143‑1 (MQ=255)                                                             ttCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGAt                                                                        >  1:461202/1‑143 (MQ=255)                                                                 gcgcACTATCAGCGCGCATAAAACAGAAGAGAAGGTAAGGGTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTaac                                                                    <  2:118493/143‑1 (MQ=255)                                                                    ctcTGCCAGCGTCCATAAAGCAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTCAATTAGAGAGGTGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATTGTCTTTAATAACAAAATAGATTaaccaa                                                                 <  2:276991/141‑1 (MQ=255)                                                                         tttggtcccATAAAACAGAAGGGAAGGTAAGTTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTaa                                                            <  2:267342/137‑1 (MQ=255)                                                                              gcacATCAAACACAAGAGAAGGAAAGGTTTACGGGGCAAACGGTGTGGATTTAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGaaa                                                       <  2:760205/140‑1 (MQ=255)                                                                                ccATAAAACAGAAGAGAAGGTAAGGTTTACGGGGCAAACGGTGTTGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTCGCAATCGTCTTTAACAACAAAATAGATTAACCAACCTAATGaaaaa                                                     <  1:607234/143‑1 (MQ=255)                                                                                  aTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTCACCTGTTTTAACAATAGTATATAACAAAAAAATAGATAAACAAACCTAATGAAAAACa                                                   >  1:70875/1‑143 (MQ=255)                                                                                  aTAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACa                                                   >  2:452738/1‑143 (MQ=255)                                                                                     aaaCAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCGAGCGATAACTAATAACCAAATCGACTAAACTGATTTAGCAATAGTCTTTGACAACCAAATAACTAAACCAACCTAACGAAAAACaaaa                                                >  2:432695/1‑142 (MQ=255)                                                                                     aaaCAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACAGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTATAGAAAACGACTGAAACAACAAAATAGAGTAACCACACTACGGAAAAGCaaaa                                                >  2:434079/1‑142 (MQ=255)                                                                                            agagaagGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTa                                         <  1:764152/143‑1 (MQ=255)                                                                                             gagaagGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTGTAGCAATAGTCATTAACAAAAAAATAGATTAACCAACCTAATGAAAAACAGATGAATTTAg                                        >  2:831935/1‑143 (MQ=255)                                                                                              agaagGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTACCCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATCAAACAACCAAAAGAAAAACAAAAGAAATTaac                                       >  1:616369/1‑141 (MQ=255)                                                                                                  ggtaaggtTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAAc                                                                                                           >  2:396615/1‑71 (MQ=255)                                                                                                  ggtaaggtTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAAc                                                                                                           <  1:396615/71‑1 (MQ=255)                                                                                                   gtaaggtTAACGGGGCAAACTGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATc                                  <  1:461054/143‑1 (MQ=255)                                                                                                   gtaaggtTAACGGGGCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATc                                  <  1:148102/143‑1 (MQ=255)                                                                                                        gtTAACGGTGCAAACGGTGTGGATTAAAGACTGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTaa                             <  1:373560/143‑1 (MQ=255)                                                                                                          tAACGGGGCTAACCGTGTGGATTAAAGGCGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGa                           <  2:675781/143‑1 (MQ=255)                                                                                                                ggCAAACGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAACGAATTGAGCAAATCATAAAGATcaaag                     >  2:749106/1‑141 (MQ=255)                                                                                                                 gCAAACGGTGTGGATTAAAGACGGGCTAGCTATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCa                    <  1:810027/143‑1 (MQ=255)                                                                                                                   caaCGGTGTGGATTAAAGACGGGCTAGCGATAACTAATTACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGc                  <  1:548387/142‑1 (MQ=255)                                                                                                                    aaCGGTGTGGATTAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCGTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCg                 >  2:386619/1‑143 (MQ=255)                                                                                                                         tgtgGATTAAAGACGTGCTAGCGATAACTAATAATCAAATCTACTATCCTGTTTTAGTAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGAtttt            <  1:245871/143‑1 (MQ=255)                                                                                                                             gATTAAAGACGGGCTAGCGATGACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTTACTTCCTCATCGCCCAACCAAATTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTgcgc        >  1:410751/1‑143 (MQ=255)                                                                                                                               ttAAAGACGGGCTAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATAAACCAACCTAATGAAAAACGAAAGAATTTAGCCAATCATGAAGATCAAGCAGCGATTGTGCGcag      >  2:278355/1‑142 (MQ=255)                                                                                                                                   aGACGAGCTAGTGATAACTAATAACTTAATCGACTAACCTGTTTTAGCAATAGTCTAGAACAACAAAAAAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACaaa  <  1:650211/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        CAACGCGGAAGGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAGACGGGATAGCGATAACTAATAACCAAATCGACTAACCTGTTTTAGCAATAGTCTTTAACAACAAAATAGATTAACCAACCTAATGAAAAACAAATGAATTTAGCCAATCATTAAGATAAATCAGCGATTTTGCGCAACAAA  >  NC_000913/2891869‑2892141

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