Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,175,151	(A)6→5	intergenic (‑125/‑237)	coaA ← / → thrU	pantothenate kinase/tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,146	0	A	.	100.0%	184.4 / NA	41	intergenic (‑120/‑242)	coaA/thrU	pantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (20/21);  total (20/21)

ATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGT  >  NC_000913/4175001‑4175288                                                                                                                                                  |                                                                                                                                               aTTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGaaa                                                                                                                                              >  2:13936/1‑148 (MQ=255)       gTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATctct                                                                                                                                      <  1:13936/149‑1 (MQ=255)                tttATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAg                                                                                                                             >  3:144781/1‑149 (MQ=255)                       tCATAAGCATGTTTCTTGCGGTTATAGGTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATATTAAG‑AAAAATCTCTGc                                                                                                                                    <  4:54216/135‑1 (MQ=255)                       tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGc                                                                                                                                    >  3:54216/1‑135 (MQ=255)                             gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGaaaa                                                                                                                <  4:144781/149‑1 (MQ=255)                                 gTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCg                                                                                                                       >  1:158476/1‑138 (MQ=255)                                 gTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCg                                                                                                                       <  2:158476/138‑1 (MQ=255)                                              aTAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAt                                                                                                           <  2:185940/137‑1 (MQ=255)                                              aTAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAt                                                                                                           >  1:185940/1‑137 (MQ=255)                                                       tCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                         <  3:134025/146‑1 (MQ=255)                                                              aGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAg                                                                                >  2:106560/1‑148 (MQ=255)                                                                                    cacaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGa                                                                                                                   <  6:47856/91‑1 (MQ=255)                                                                                    cacaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGa                                                                                                                   >  5:47856/1‑91 (MQ=255)                                                                                      caGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGa                                                              >  5:11148/1‑142 (MQ=255)                                                                                      caGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGa                                                              <  6:11148/142‑1 (MQ=255)                                                                                         ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaa                                                                                    <  3:149156/117‑1 (MQ=255)                                                                                         ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaa                                                                                    >  4:149156/1‑117 (MQ=255)                                                                                            tGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                       >  6:84523/1‑95 (MQ=255)                                                                                             gTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCt                                                                                                                           <  4:124565/74‑1 (MQ=255)                                                                                             gTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCt                                                                                                                           >  3:124565/1‑74 (MQ=255)                                                                                                 aaTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaa                                                                                   >  5:46772/1‑110 (MQ=255)                                                                                                 aaTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaa                                                                                   <  6:46772/110‑1 (MQ=255)                                                                                                            gggCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGa                                                              <  4:15486/120‑1 (MQ=255)                                                                                                                aggaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaa                                                                                    <  3:119822/94‑1 (MQ=255)                                                                                                                aggaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaa                                                                                    >  4:119822/1‑94 (MQ=255)                                                                                                                      gTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGt                                                                       >  8:121024/1‑101 (MQ=255)                                                                                                                      gTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGt                                                                       <  7:121024/101‑1 (MQ=255)                                                                                                                                    ggaggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGtt         >  3:99180/1‑149 (MQ=255)                                                                                                                                     gaggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGt                                                                       <  6:153022/86‑1 (MQ=255)                                                                                                                                     gaggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGt                                                                       >  5:153022/1‑86 (MQ=255)                                                                                                                                     gaggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAg                                                 <  7:116020/108‑1 (MQ=255)                                                                                                                                     gaggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAg                                                 >  8:116020/1‑108 (MQ=255)                                                                                                                                      aggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGcc                >  2:34079/1‑140 (MQ=255)                                                                                                                                      aggaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGcc                <  1:34079/140‑1 (MQ=255)                                                                                                                                        gaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTa                                                                    >  3:93628/1‑86 (MQ=255)                                                                                                                                        gaATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTa                                                                    <  4:93628/86‑1 (MQ=255)                                                                                                                                          ataataAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAg                                                                   >  6:119155/1‑85 (MQ=255)                                                                                                                                          ataataAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAg                                                                   <  5:119155/85‑1 (MQ=255)                                                                                                                                           taataaG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGt  <  4:99180/149‑1 (MQ=255)                                                                                                                                            aataaG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGtt                                     <  5:86914/113‑1 (MQ=255)                                                                                                                                            aataaG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGtt                                     >  6:86914/1‑113 (MQ=255)                                                                                                                                                  |                                                                                                                                               ATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTGT  >  NC_000913/4175001‑4175288

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑