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%	166.1 / NA	51	*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 (26/25);  total (26/25)

ACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTG  >  NC_000913/922257‑922495                                                                                                              |                                                                                                                                  aCGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                      <  2:1235328/139‑1 (MQ=255)             aGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACt                                                                                          >  2:1097063/1‑139 (MQ=255)                  aTGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                        >  2:55982/1‑104 (MQ=255)                  aTGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGc                                                                                                                        <  1:55982/104‑1 (MQ=255)                    gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGa                                                                                   >  1:1292863/1‑139 (MQ=255)                       ctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTg                                                                                >  1:947251/1‑139 (MQ=255)                            aCGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCttt                                                                           <  2:30073/139‑1 (MQ=255)                             cGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCtttt                                                                          <  1:1097063/139‑1 (MQ=255)                              gagaGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                         <  2:947251/139‑1 (MQ=255)                                       ttttGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt                                                                <  2:1292863/139‑1 (MQ=255)                                          tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGAc                                                                                                         <  2:671188/95‑1 (MQ=255)                                          tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGAc                                                                                                         >  1:671188/1‑95 (MQ=255)                                           gagaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                              >  2:363415/1‑137 (MQ=255)                                           gagaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                              <  1:363415/137‑1 (MQ=255)                                                 aaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                                                    >  2:695707/1‑125 (MQ=255)                                                 aaaaTGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                                                    <  1:695707/125‑1 (MQ=255)                                                              aTCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                                                           >  2:273648/1‑121 (MQ=255)                                                              aTCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                                                           <  1:273648/121‑1 (MQ=255)                                                                     tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                       >  1:1067892/1‑102 (MQ=255)                                                                     tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTGCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                       <  2:1067892/102‑1 (MQ=255)                                                                                         gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                              <  1:1335616/123‑1 (MQ=255)                                                                                         gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                              >  2:357341/1‑123 (MQ=255)                                                                                         gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                              >  2:1335616/1‑123 (MQ=255)                                                                                         gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                              <  1:357341/123‑1 (MQ=255)                                                                                          tACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                                                    <  2:60635/84‑1 (MQ=255)                                                                                          tACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                                                    >  1:60635/1‑84 (MQ=255)                                                                                            cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                                                   >  1:933630/1‑51 (MQ=255)                                                                                            cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                                                   <  2:933630/51‑1 (MQ=255)                                                                                            cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACcatccatc           >  1:429124/1‑139 (MQ=255)                                                                                            cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACcatccatc           >  1:579832/1‑139 (MQ=255)                                                                                               aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGa        >  1:176026/1‑139 (MQ=255)                                                                                                aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa       <  2:176026/139‑1 (MQ=255)                                                                                                 tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTa                               >  2:65703/1‑114 (MQ=255)                                                                                                 tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTa                               <  1:65703/114‑1 (MQ=255)                                                                                                 tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa       <  1:412570/138‑1 (MQ=255)                                                                                                 tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa       >  2:412570/1‑138 (MQ=255)                                                                                                 tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAt      >  1:527745/1‑139 (MQ=255)                                                                                                 tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAt      >  1:1435136/1‑139 (MQ=255)                                                                                                   agaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtgg    <  2:527745/139‑1 (MQ=255)                                                                                                    gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCATAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                             >  1:984488/1‑97 (MQ=255)                                                                                                    gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                       <  1:964204/71‑1 (MQ=255)                                                                                                    gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                       >  2:964204/1‑71 (MQ=255)                                                                                                    gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                             >  1:977713/1‑97 (MQ=255)                                                                                                    gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                             <  2:977713/97‑1 (MQ=255)                                                                                                    gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                             <  2:984488/97‑1 (MQ=255)                                                                                                       agaagaGTTATGCGACTGCCGCTTCTACCTCGACGGGCACaa                                                                                                 <  1:806354/42‑1 (MQ=255)                                                                                                       agaagaGTTATGCGACTGCCGCTTCTACCTCGACGGGCACaa                                                                                                 >  2:806354/1‑42 (MQ=38)                                                                                                         aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                       >  2:750201/1‑66 (MQ=255)                                                                                                         aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                       <  1:750201/66‑1 (MQ=255)                                                                                                         aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggtg  >  2:789756/1‑135 (MQ=255)                                                                                                         aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggtg  <  1:789756/135‑1 (MQ=255)                                                                                                              |                                                                                                                                  ACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTG  >  NC_000913/922257‑922495

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