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%	103.8 / NA	34	*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 (17/17);  total (17/17)

GGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGT  >  NC_000913/922231‑922494                                                                                                                                        |                                                                                                                                 ggCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTAt                                                                                                                               >  2:214703/1‑139 (MQ=255)      gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGAc                                                                                                                          >  2:405859/1‑139 (MQ=255)            cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGc                                                                                                                    <  1:735333/139‑1 (MQ=255)                         gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                       <  2:69338/139‑1 (MQ=255)                           aCGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                     <  1:781817/139‑1 (MQ=255)                             gACTGGATCCGGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGca                                                                                                   >  2:227069/1‑139 (MQ=255)                                   aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa                                                                                             >  1:723406/1‑139 (MQ=255)                                         tAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACgg                                                                                                      >  2:679631/1‑124 (MQ=255)                                         tAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACgg                                                                                                      <  1:679631/124‑1 (MQ=255)                                                  tGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                     <  1:189306/116‑1 (MQ=255)                                                  tGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                     >  2:189306/1‑116 (MQ=255)                                                                  tttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtg                                                              <  1:214703/139‑1 (MQ=255)                                                                   ttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                             <  1:38294/139‑1 (MQ=255)                                                                     gagaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCAGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAc                                                           >  1:484848/1‑139 (MQ=255)                                                                               tGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                                                 <  1:634762/139‑1 (MQ=255)                                                                                                  cATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTTGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTa                              <  1:4315/139‑1 (MQ=255)                                                                                                      ttGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGt                          >  2:73879/1‑139 (MQ=255)                                                                                                                gATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAAcca                <  2:678619/139‑1 (MQ=255)                                                                                                                          aTGAGTCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACGCTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa      <  2:720162/139‑1 (MQ=255)                                                                                                                          aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac                                                                                            >  2:733788/1‑53 (MQ=255)                                                                                                                          aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac                                                                                            <  1:733788/53‑1 (MQ=255)                                                                                                                          aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa      <  1:73879/139‑1 (MQ=255)                                                                                                                           tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa      <  2:35350/138‑1 (MQ=255)                                                                                                                           tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa      >  1:35350/1‑138 (MQ=255)                                                                                                                           tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAt     >  2:159814/1‑139 (MQ=255)                                                                                                                           tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCCTTTAGCGTTCTGTAACCATCCATCTGAAt     >  2:609996/1‑139 (MQ=255)                                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggt  <  1:405859/139‑1 (MQ=255)                                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggt  >  2:530949/1‑139 (MQ=255)                                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggt  >  2:366905/1‑139 (MQ=255)                                                                                                                              gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggt  >  2:212368/1‑139 (MQ=255)                                                                                                                                 agtagagTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                     >  2:674549/4‑37 (MQ=255)                                                                                                                                 agtagagTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                     <  1:674549/34‑1 (MQ=255)                                                                                                                                   aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGt                          >  2:407647/1‑110 (MQ=255)                                                                                                                                   aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGt                          <  1:407647/110‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 GGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGT  >  NC_000913/922231‑922494

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