Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	922,366	C→T	*75* (TAG→TAA)	cspD ←	inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	922,366	0	C	T	100.0%	130.2 / NA	41	*75* (TAG→TAA)	cspD	inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/22);  total (19/22)

GTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAA  >  NC_000913/922236‑922490                                                                                                                                   |                                                                                                                             gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGAc                                                                                                                      <  2:386662/139‑1 (MQ=255) gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGAc                                                                                                                      >  2:1434597/1‑139 (MQ=255)         cACGAAATGCGGCACGACTTGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCtt                                                                                                              >  2:1072593/1‑139 (MQ=255)         cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCtt                                                                                                              <  1:36942/139‑1 (MQ=255)                      aCGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                 <  1:97024/139‑1 (MQ=255)                           tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACaa                                                                                            <  1:1072593/139‑1 (MQ=255)                                         gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTATGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGa                                                                              <  1:457934/139‑1 (MQ=255)                                         gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACATTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGa                                                                              >  1:1154825/1‑139 (MQ=255)                                         gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGa                                                                              >  2:335479/1‑139 (MQ=255)                                              ggCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTa                                                                                                                            >  1:1585921/1‑88 (MQ=255)                                              ggCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTa                                                                                                                            <  2:1585921/88‑1 (MQ=255)                                                               tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGa                                                        >  1:1509923/1‑139 (MQ=255)                                                                  gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                   <  2:1395612/77‑1 (MQ=255)                                                                  gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                   >  1:1395612/1‑77 (MQ=255)                                                                  gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                           <  2:1379957/133‑1 (MQ=255)                                                                  gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                           >  1:1379957/1‑133 (MQ=255)                                                                     aaaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                  <  2:828029/123‑1 (MQ=255)                                                                     aaaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                  >  1:828029/1‑123 (MQ=255)                                                                             cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATg                                                                                                                          <  1:1425816/59‑1 (MQ=255)                                                                             cAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATg                                                                                                                          >  2:1425816/1‑59 (MQ=255)                                                                                  gATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                   <  1:1486234/61‑1 (MQ=255)                                                                                  gATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                   >  2:1486234/1‑61 (MQ=255)                                                                                      gggCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACGTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCa                                 <  1:804488/139‑1 (MQ=255)                                                                                             cATTTTGCCTTTAGGATGTACACAATGAGACAGAAGTGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCTCCTGGTGGACATCAAACTGAAAGGATTGTCAAGCTTTTa                          >  1:567984/1‑139 (MQ=255)                                                                                             cATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTa                          <  2:697816/139‑1 (MQ=255)                                                                                                tttGCCTTTAGGATGTACACAATGAGAGAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTCGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCg                       <  1:928469/139‑1 (MQ=255)                                                                                                  tGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                                              <  1:204873/66‑1 (MQ=255)                                                                                                  tGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                                              >  2:204873/1‑66 (MQ=255)                                                                                                  tGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGtt                     >  1:1334840/1‑139 (MQ=255)                                                                                                                acacAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAAccatccat       <  2:90161/139‑1 (MQ=255)                                                                                                                acacAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAAccatccat       <  2:551983/139‑1 (MQ=255)                                                                                                                acacAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAAccatccat       <  2:1509923/139‑1 (MQ=255)                                                                                                                 cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACcatccatc      >  1:14217/1‑139 (MQ=255)                                                                                                                   caATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTg    <  1:1434597/139‑1 (MQ=255)                                                                                                                    aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGa   >  2:206327/1‑139 (MQ=255)                                                                                                                         gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTg                                                                           <  2:467181/62‑1 (MQ=255)                                                                                                                         gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTg                                                                           >  1:467181/1‑62 (MQ=255)                                                                                                                         gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa  <  1:1251679/135‑1 (MQ=255)                                                                                                                         gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa  >  2:1251679/1‑135 (MQ=255)                                                                                                                              aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa  >  2:92908/1‑130 (MQ=255)                                                                                                                              aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa  <  1:92908/130‑1 (MQ=255)                                                                                                                                   |                                                                                                                             GTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAA  >  NC_000913/922236‑922490

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