Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,515,726	T→C	Y63C (TAT→TGT)	dam ←	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,515,726	0	T	C	100.0%	202.5 / NA	55	Y63C (TAT→TGT)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/24);  total (31/24)

CCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTA  >  NC_000913/3515592‑3515851                                                                                                                                       |                                                                                                                              cccGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACagag                                                                                                                           >  2:930166/1‑139 (MQ=255)            ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                <  1:266060/139‑1 (MQ=255)                        ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGatat                                                                                                    >  2:1304421/1‑139 (MQ=255)                           ttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCgg                                                                                                 <  1:930166/139‑1 (MQ=255)                                ctctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                      >  2:1663365/1‑129 (MQ=255)                                ctctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTGCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                      <  1:1663365/129‑1 (MQ=255)                                   ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGata                                                                                         >  2:1709487/1‑139 (MQ=255)                                           aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGatat                                                                                                    >  1:553692/1‑120 (MQ=255)                                           aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGatat                                                                                                    <  2:553692/120‑1 (MQ=255)                                           aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACgaga                                                                                 >  1:1647394/1‑139 (MQ=255)                                             tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                               >  2:1450882/1‑123 (MQ=255)                                             tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                               <  1:1450882/123‑1 (MQ=255)                                                  gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                          <  2:241651/139‑1 (MQ=255)                                                    aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgt                                                                        >  1:174686/1‑139 (MQ=255)                                                         tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGAt                                                                                                                     <  1:1416343/89‑1 (MQ=255)                                                         tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGAt                                                                                                                     >  2:1416343/1‑89 (MQ=255)                                                          cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                 <  1:1235847/92‑1 (MQ=255)                                                          cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                 >  2:1235847/1‑92 (MQ=255)                                                             cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                <  2:315321/90‑1 (MQ=255)                                                             cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                >  1:315321/1‑90 (MQ=255)                                                             cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                               >  2:824183/1‑107 (MQ=255)                                                             cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                               <  1:824183/107‑1 (MQ=255)                                                                                 aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaa                                           <  1:881228/139‑1 (MQ=255)                                                                                      cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                     >  2:780383/1‑76 (MQ=255)                                                                                      cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                     <  1:780383/76‑1 (MQ=255)                                                                                      cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                      >  1:115612/1‑139 (MQ=255)                                                                                      cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                      >  2:1028673/1‑139 (MQ=255)                                                                                       tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                 >  1:268070/1‑63 (MQ=255)                                                                                       tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                 <  2:268070/63‑1 (MQ=255)                                                                                            tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                >  1:531553/1‑139 (MQ=255)                                                                                              cGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAg                              >  1:483442/1‑139 (MQ=255)                                                                                                   tGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                >  2:392364/1‑52 (MQ=255)                                                                                                   tGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                <  1:392364/52‑1 (MQ=255)                                                                                                     tacgtacTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                       <  2:174686/139‑1 (MQ=255)                                                                                                          acTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                >  1:127845/1‑61 (MQ=255)                                                                                                          acTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                <  2:127845/61‑1 (MQ=255)                                                                                                          acTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                          <  1:1524906/83‑1 (MQ=255)                                                                                                          acTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                          >  2:1524906/1‑83 (MQ=255)                                                                                                           cTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                <  2:349861/44‑1 (MQ=255)                                                                                                           cTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                >  1:349861/1‑44 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                                                    >  1:1703558/1‑100 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                                                    <  2:1703558/100‑1 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                       <  2:1027854/129‑1 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                       >  1:1027854/1‑129 (MQ=255)                                                                                                                caGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaa            >  2:914015/1‑139 (MQ=255)                                                                                                                        aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt    >  1:1207340/1‑139 (MQ=255)                                                                                                                        aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt    <  2:1207340/139‑1 (MQ=255)                                                                                                                        aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt    >  2:57966/1‑139 (MQ=255)                                                                                                                          cTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTa  >  1:1411378/1‑139 (MQ=255)                                                                                                                             cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCa                               <  1:720803/107‑1 (MQ=255)                                                                                                                             cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCa                               >  2:720803/1‑107 (MQ=255)                                                                                                                                  tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                            <  1:292040/57‑1 (MQ=255)                                                                                                                                  tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                            >  2:292040/1‑57 (MQ=255)                                                                                                                                  tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg         >  2:1642465/1‑124 (MQ=255)                                                                                                                                  tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg         <  1:1642465/124‑1 (MQ=255)                                                                                                                                       |                                                                                                                              CCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTA  >  NC_000913/3515592‑3515851

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