Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	3,904,745	A→C	intergenic (‑242/+503)	dppA ← / ← ECW_RS25300	dipeptide ABC transporter substrate‑binding protein DppA/RtT sRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	3,904,745	0	A	C	100.0%	206.3 / NA	51	intergenic (‑242/+503)	dppA/ECW_RS25300	dipeptide ABC transporter substrate‑binding protein DppA/RtT sRNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (32/19);  total (32/19)

ACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATT  >  NC_017635/3904698‑3904791                                                |                                               aCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCacacac                                              <  2:1830923/50‑1 (MQ=255)  cGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             <  2:1379147/50‑1 (MQ=255)  cGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCacacac                                              >  2:594635/1‑49 (MQ=255)   gTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  2:1655500/1‑49 (MQ=255)    tCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  2:1481683/1‑48 (MQ=255)    tCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCAcacacaaa                                           <  2:46687/50‑1 (MQ=255)     cAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  1:482311/1‑47 (MQ=255)      aaTACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  1:122651/1‑46 (MQ=255)      aaTACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  2:1652536/1‑46 (MQ=255)      aaTACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  1:1848364/1‑46 (MQ=255)      aaTACAGCCAACCGGGATTTACACCAACGGTGAGAATCcacacaca                                             >  1:900744/1‑46 (MQ=255)      aaTACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAg                                          <  2:919884/49‑1 (MQ=255)       aTACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGAt                                        <  2:1287584/50‑1 (MQ=255)       aTACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGAt                                        <  2:1426550/50‑1 (MQ=255)        tACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGAt                                        >  2:1725248/1‑49 (MQ=255)           aGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTaa                                     >  1:575248/1‑49 (MQ=255)               aaCCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTaaaaaaa                                <  1:591597/50‑1 (MQ=255)               aaCCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTaaaaaaa                                <  2:1297042/50‑1 (MQ=255)               aaCCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTaaaaaaa                                >  2:1812072/1‑50 (MQ=255)                aCCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAAc                               <  1:778216/50‑1 (MQ=255)                 ccGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACt                              >  2:1613679/1‑50 (MQ=255)                 ccGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACt                              <  2:145090/50‑1 (MQ=255)                  cGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACtt                             >  2:1860590/1‑50 (MQ=255)                  cGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACtt                             >  2:1897324/1‑50 (MQ=255)                   gggATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACtt                             >  1:17959/1‑49 (MQ=255)                   gggATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTc                            <  2:1609786/50‑1 (MQ=255)                    ggATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCa                           <  2:1744876/50‑1 (MQ=255)                    ggATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCa                           <  1:1016459/50‑1 (MQ=255)                    ggATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCa                           <  1:2054802/50‑1 (MQ=255)                      aTTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCaaa                         >  2:1222878/1‑50 (MQ=255)                      aTTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCaaa                         <  2:1055124/50‑1 (MQ=255)                       tttACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCaaa                         >  1:1219550/1‑49 (MQ=255)                        ttACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCaaa                         >  1:851910/1‑48 (MQ=255)                          acacCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAg                      <  2:2072429/49‑1 (MQ=255)                           cacCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGc                     >  2:1789268/1‑49 (MQ=255)                             ccAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAt                  >  1:1192179/1‑50 (MQ=255)                             ccAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAt                  <  1:2023164/50‑1 (MQ=255)                              cAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAtt                 >  2:1807450/1‑50 (MQ=255)                              cAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAtt                 <  1:1570762/50‑1 (MQ=255)                              cAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAtt                 <  2:2069029/50‑1 (MQ=255)                               aaCGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAttt                >  1:399151/1‑50 (MQ=255)                               aaCGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAttt                >  2:1857299/1‑50 (MQ=255)                               aaCGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTAttt                >  2:264088/1‑50 (MQ=255)                                aCGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTg               >  1:1993256/1‑50 (MQ=255)                                 cGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTg               >  2:1966434/1‑49 (MQ=255)                                 cGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTg               >  2:2074494/1‑49 (MQ=255)                                 cGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTg               >  1:131880/1‑49 (MQ=255)                                 cGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTg               >  1:875746/1‑49 (MQ=255)                                  ggTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTg               >  1:1517245/1‑48 (MQ=255)                                   gTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            >  2:916327/1‑50 (MQ=255)                                   gTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            <  2:863903/50‑1 (MQ=255)                                   gTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            <  1:80330/50‑1 (MQ=255)                                   gTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            >  2:2011020/1‑50 (MQ=255)                                   gTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            >  2:244605/1‑50 (MQ=255)                                    tGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            >  2:566748/1‑49 (MQ=255)                                     gagaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTCTTTgcag            >  1:398861/1‑48 (MQ=39)                                     gagaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcag            >  2:1858162/1‑48 (MQ=255)                                       gaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcagca          >  2:1538002/1‑48 (MQ=255)                                       gaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcagca          >  1:1973152/1‑48 (MQ=255)                                        aaTCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTgcagca          >  1:889846/1‑47 (MQ=255)                                          tCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAgaga      >  2:650709/1‑49 (MQ=255)                                          tCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGaa     <  1:410056/50‑1 (MQ=255)                                            cacacacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAt   >  1:975303/1‑50 (MQ=255)                                            cacacacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAt   <  2:1234709/50‑1 (MQ=255)                                            cacacacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAt   >  1:787025/1‑50 (MQ=255)                                            cacacacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAt   >  1:401240/1‑50 (MQ=255)                                             acacacaaAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAtt  >  2:906744/1‑50 (MQ=255)                                             acacacaaAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAtt  >  1:496000/1‑50 (MQ=255)                                                |                                               ACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATT  >  NC_017635/3904698‑3904791

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