Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,681,163	A→C	D63A (GAC→GCC)	folM →	dihydromonapterin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,681,163	0	A	C	100.0%	188.8 / NA	52	D63A (GAC→GCC)	folM	dihydromonapterin reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/24);  total (28/24)

GCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGACGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATC  >  NC_000913/1681031‑1681299                                                                                                                                     |                                                                                                                                         gCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGtgtg                                                                                                                                   >  2:678716/1‑140 (MQ=255)           cgcATGGCATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTg                                                                                                                         >  1:404885/1‑140 (MQ=255)                 gCATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATg                                                                                                                   <  2:104250/140‑1 (MQ=255)                 gCATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATg                                                                                                                   <  2:404885/140‑1 (MQ=255)                  cATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGa                                                                                                                  >  1:316825/1‑140 (MQ=255)                  cATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGa                                                                                                                  >  1:372418/1‑140 (MQ=255)                    tttCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAg                                                                                                                <  2:502418/140‑1 (MQ=255)                      tCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTa                                                                                                              >  1:915723/1‑140 (MQ=255)                      tCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTa                                                                                                              >  1:857910/1‑140 (MQ=255)                            aTCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTaaaa                                                                                                        >  2:910446/1‑140 (MQ=255)                              cAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      <  1:386522/140‑1 (MQ=255)                              cAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      <  2:857910/140‑1 (MQ=255)                              cAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      <  1:569974/140‑1 (MQ=255)                              cAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      <  2:316825/140‑1 (MQ=255)                              cAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      <  2:228086/140‑1 (MQ=255)                                 aaGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCAc                                                                                                   >  1:211879/1‑140 (MQ=255)                                  aGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATg                                                                                                                   >  1:1089146/1‑123 (MQ=255)                                  aGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATg                                                                                                                   <  2:1089146/123‑1 (MQ=255)                                          gTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCt                                                                                          <  2:1028331/140‑1 (MQ=255)                                           tGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      <  1:181913/127‑1 (MQ=255)                                           tGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAg                                                                                                      >  2:181913/1‑127 (MQ=255)                                                     cTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTAttt                                                                               >  1:1065403/1‑140 (MQ=255)                                                       aTCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTg                                                                             <  2:915723/140‑1 (MQ=255)                                                         cGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCa                                                                           <  2:1065403/140‑1 (MQ=255)                                                          ggACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCAt                                                                          <  1:999548/140‑1 (MQ=255)                                                               cacTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGc                                                                     <  1:910446/140‑1 (MQ=255)                                                                acTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGcc                                                                    <  1:86371/140‑1 (MQ=255)                                                                   aTCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGt                                                                 >  2:564447/1‑140 (MQ=255)                                                                          cATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGa                                                          >  2:616773/1‑140 (MQ=255)                                                                           aTTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGAt                                                         >  2:550493/1‑140 (MQ=255)                                                                           aTTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGAt                                                         <  1:649193/140‑1 (MQ=255)                                                                            ttGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATg                                                        >  2:827541/1‑140 (MQ=255)                                                                              gATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGc                                                      >  2:902987/1‑140 (MQ=255)                                                                               aTGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCg                                                     <  1:564447/140‑1 (MQ=255)                                                                                   aCTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGaaa                                                 >  2:618376/1‑140 (MQ=255)                                                                                      gATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAc                                              >  1:369952/1‑140 (MQ=255)                                                                                      gATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAc                                              >  2:123283/1‑140 (MQ=255)                                                                                       aTTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             >  2:353984/1‑140 (MQ=255)                                                                                       aTTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             >  2:676433/1‑140 (MQ=255)                                                                                       aTTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             >  2:842936/1‑140 (MQ=255)                                                                                       aTTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAAcc                                             <  1:902987/140‑1 (MQ=255)                                                                                            tGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTg                                        <  2:821288/140‑1 (MQ=255)                                                                                                 gTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGAGGATGGCGGAAAAACCGGGTGCGCCa                                   >  1:765968/1‑140 (MQ=255)                                                                                                   gcgcAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACt                                 <  2:617161/140‑1 (MQ=255)                                                                                                   gcgcAGTGTATTCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACt                                 <  1:676433/140‑1 (MQ=255)                                                                                                             ttCAGGCTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTa                       <  1:550493/140‑1 (MQ=255)                                                                                                                   cTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGc                              >  1:1001804/1‑127 (MQ=255)                                                                                                                   cTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGc                              <  2:1001804/127‑1 (MQ=255)                                                                                                                   cTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCt                 >  2:712690/1‑140 (MQ=255)                                                                                                                   cTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCt                 >  2:630701/1‑140 (MQ=255)                                                                                                                   cTGATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCt                 >  2:403914/1‑140 (MQ=255)                                                                                                                     gATTTTTCGACCAACGCCGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTg               >  2:326822/1‑140 (MQ=255)                                                                                                                                  acgccgGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc  <  1:75848/140‑1 (MQ=255)                                                                                                                                  acgccgGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATc  <  2:369952/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                         GCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTGATTGTCAGCTATCGGACACACTATCCAGCCATTGATGGACTGATTAATGCAGGTGCGCAGTGTATTCAGGCTGATTTTTCGACCAACGACGGTGTGATGGCGTTTGCCGATGAAGTACTAAAAAGCACCCATGGTCTGCGTGCTATTTTGCATAACGCCAGTGCGTGGATGGCGGAAAAACCGGGTGCGCCACTGGCCGACGTACTGGCTTGCATGATGCAGATC  >  NC_000913/1681031‑1681299

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