Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,576,388	A→G	S132P (TCC→CCC)	ydeN ←	sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,576,388	0	A	G	100.0%	142.5 / NA	44	S132P (TCC→CCC)	ydeN	sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/21);  total (23/21)

AATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAG  >  NZ_CP009273/1576302‑1576471                                                                                       |                                                                                    aaTACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAac                                                                                  <  1:1211637/90‑1 (MQ=255) aaTACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAac                                                                                  <  1:1248250/90‑1 (MQ=255)   tACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                                <  2:834209/90‑1 (MQ=255)   tACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                                <  2:207767/90‑1 (MQ=255)        tCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                                                           <  1:894458/90‑1 (MQ=255)          cTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                         <  1:49849/90‑1 (MQ=255)            gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                       >  1:1192883/1‑90 (MQ=255)            gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                       >  2:117957/1‑90 (MQ=255)            gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                       >  1:878106/1‑90 (MQ=255)                aTCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATa                                                                   <  2:917655/90‑1 (MQ=255)                 tCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAg                                                                  >  1:414699/1‑90 (MQ=255)                 tCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAg                                                                  >  1:1548568/1‑90 (MQ=255)                   ggTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGcc                                                                >  1:1056354/1‑90 (MQ=255)                      aTTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCg                                                                              >  2:1449004/1‑73 (MQ=255)                      aTTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCg                                                                              <  1:1449004/73‑1 (MQ=255)                            aTAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTaaa                                                       >  1:850572/1‑90 (MQ=255)                              aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGtt                                                             >  2:932888/1‑82 (MQ=255)                              aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGtt                                                             <  1:932888/82‑1 (MQ=255)                              aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                     >  2:265012/1‑90 (MQ=255)                                 cacCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                  >  2:273123/1‑90 (MQ=255)                                   ccAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACg                                                >  2:513055/1‑90 (MQ=255)                                     aaaGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGcc                                              >  2:1205713/1‑90 (MQ=255)                                           ggcgggAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATc                                        >  1:809702/1‑90 (MQ=255)                                                      gACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCg                                                               <  1:952210/56‑1 (MQ=255)                                                      gACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCg                                                               >  2:952210/1‑56 (MQ=255)                                                      gACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaaa                             <  1:779906/90‑1 (MQ=255)                                                          ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCg                                                                              <  2:63179/37‑1 (MQ=38)                                                          ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCg                                                                              >  1:63179/1‑37 (MQ=255)                                                          ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACAt                                                                    <  2:1572251/47‑1 (MQ=255)                                                          ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACAt                                                                    >  1:1572251/1‑47 (MQ=255)                                                          ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAg                         >  2:590297/1‑90 (MQ=255)                                                              attattGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                     <  1:217349/90‑1 (MQ=255)                                                              attattGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                     <  2:1398114/90‑1 (MQ=255)                                                              attattGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                     <  2:177243/90‑1 (MQ=255)                                                                      ggcgCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa             >  1:439184/1‑90 (MQ=255)                                                                         gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGtt                                                             <  2:1438883/39‑1 (MQ=39)                                                                         gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGtt                                                             >  1:1438883/1‑39 (MQ=255)                                                                         gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt          >  2:206581/1‑90 (MQ=255)                                                                         gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt          >  2:381604/1‑90 (MQ=255)                                                                         gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt          >  1:150936/1‑90 (MQ=255)                                                                              aGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtgtg     <  1:69696/90‑1 (MQ=255)                                                                               ggggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc    <  1:908538/90‑1 (MQ=255)                                                                                gggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCa   <  1:513055/90‑1 (MQ=255)                                                                                 ggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAg  <  1:1179252/90‑1 (MQ=255)                                                                                       |                                                                                    AATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAG  >  NZ_CP009273/1576302‑1576471

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