Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,066,175	T→C	intergenic (‑411/+130)	cobU ← / ← insH‑6	cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,066,175	0	T	C	100.0%	158.7 / NA	44	intergenic (‑411/+130)	cobU/insH‑6	cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/21);  total (23/21)

CGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGC  >  NC_000913/2066040‑2066290                                                                                                                                        |                                                                                                                    cGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCgggg                                                                                                                 <  1:246302/140‑1 (MQ=255) cGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCgg                                                                                                                   >  1:164470/1‑138 (MQ=255)    gCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGaaa                                                                                                              <  1:54730/140‑1 (MQ=255)        tgtgTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAAttct                                                                                                          >  2:949800/1‑140 (MQ=255)            ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                                                                                                      >  2:935679/1‑140 (MQ=255)             tCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACCATGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                                                                                                     <  2:682086/140‑1 (MQ=255)              cTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                                                                                    >  2:813970/1‑140 (MQ=255)                   aTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                                                               <  2:228600/140‑1 (MQ=255)                   aTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                                                               <  1:722790/140‑1 (MQ=255)                                            cGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                      >  2:47924/1‑140 (MQ=255)                                               cAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                   <  1:813970/140‑1 (MQ=255)                                                    cTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAt                                                              >  2:660082/1‑140 (MQ=255)                                                            tAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                      >  2:703287/1‑140 (MQ=255)                                                            tAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                      >  1:112033/1‑140 (MQ=255)                                                            tAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                      >  1:649246/1‑140 (MQ=255)                                                            tAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                      >  2:1039051/1‑140 (MQ=255)                                                               aaGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCg                                                   <  1:949800/140‑1 (MQ=255)                                                                      tGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                            >  1:900462/1‑140 (MQ=255)                                                                      tGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                            >  2:724452/1‑140 (MQ=255)                                                                       gTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                           <  2:166549/140‑1 (MQ=255)                                                                            acaacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                      <  1:1039051/140‑1 (MQ=255)                                                                            acaacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                      <  1:1090355/140‑1 (MQ=255)                                                                            acaacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                      <  1:172332/140‑1 (MQ=255)                                                                                cGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                  <  1:724452/140‑1 (MQ=255)                                                                                cGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                  >  1:546991/1‑140 (MQ=255)                                                                                 gTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                 >  1:853040/1‑140 (MQ=255)                                                                                 gTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                 <  2:649246/140‑1 (MQ=255)                                                                                 gTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                 <  2:104516/140‑1 (MQ=255)                                                                                          ttATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGa                        >  2:880446/1‑140 (MQ=255)                                                                                           tATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                       <  2:900462/140‑1 (MQ=255)                                                                                            aTCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                      >  2:302482/1‑124 (MQ=255)                                                                                            aTCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                      <  1:302482/124‑1 (MQ=255)                                                                                            aTCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                      <  1:621252/140‑1 (MQ=255)                                                                                            aTCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                      >  1:313596/1‑140 (MQ=255)                                                                                             tCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                  <  1:875199/127‑1 (MQ=255)                                                                                             tCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                  >  2:875199/1‑127 (MQ=255)                                                                                             tCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCg                     >  2:1008246/1‑140 (MQ=255)                                                                                                  aaTGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttct                >  2:1069051/1‑140 (MQ=255)                                                                                                    tggtggCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc              <  1:660082/140‑1 (MQ=255)                                                                                                    tggtggCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc              <  1:880446/140‑1 (MQ=255)                                                                                                            aCTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa      >  1:258169/1‑140 (MQ=21)                                                                                                            aCTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa      >  2:144681/1‑140 (MQ=18)                                                                                                              tATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg    <  1:1035605/140‑1 (MQ=18)                                                                                                                tGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc  >  2:83727/1‑140 (MQ=11)                                                                                                                                        |                                                                                                                    CGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGC  >  NC_000913/2066040‑2066290

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: ‑