Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,640,455	C→T	N18N (AAC→AAT)	yjjY →	protein YjjY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,640,455	0	C	T	100.0%	112.5 / NA	33	N18N (AAC→AAT)	yjjY	protein YjjY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/19);  total (14/19)

AACTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAACATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGACAT  >  NC_000913/4640326‑4640590                                                                                                                                  |                                                                                                                                        aaCTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCt                                                                                                                               <  1:146487/140‑1 (MQ=255)       aTCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCg                                                                                                                         >  1:514795/1‑140 (MQ=255)           aaaCAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTggg                                                                                                                     <  1:343479/140‑1 (MQ=255)                 gAAGTACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACttt                                                                                                               >  2:184103/1‑140 (MQ=255)                     tACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCt                                                                                                           <  1:480545/140‑1 (MQ=255)                     tACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCt                                                                                                           <  1:507302/140‑1 (MQ=255)                               ccTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGa                                                                                                 >  2:326291/1‑140 (MQ=255)                                       gcctgATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGAc                                                                                                <  2:543618/133‑1 (MQ=255)                                       gcctgATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGAc                                                                                                >  1:543618/1‑133 (MQ=255)                                       gcctgATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCa                                                                                         >  1:99051/1‑140 (MQ=255)                                           gATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCaaaaa                                                                                     <  2:514795/140‑1 (MQ=255)                                                  gctgcAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGt                                                                              >  2:164701/1‑140 (MQ=255)                                                   ctgcAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTc                                                                             >  2:192938/1‑140 (MQ=255)                                                      cAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATc                                                                          >  2:393048/1‑140 (MQ=255)                                                         aTTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTa                                                                       <  1:356863/140‑1 (MQ=255)                                                              cATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCAt                                                                  >  2:293837/1‑140 (MQ=255)                                                                  aTCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTa                                                              <  2:110677/140‑1 (MQ=255)                                                                    cGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAAc                                                            <  2:126689/140‑1 (MQ=255)                                                                         tAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCaa                                                       >  1:95359/1‑140 (MQ=255)                                                                            cATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACAt                                                    >  2:293487/1‑140 (MQ=255)                                                                              tGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATaa                                                  <  2:60041/140‑1 (MQ=255)                                                                              tGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATaa                                                  <  1:293487/140‑1 (MQ=255)                                                                              tGACTAAAGTACGTAATTGCGTTCTTGATGCACGTTCGATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATaa                                                  <  1:500276/140‑1 (MQ=255)                                                                               gACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAAc                                                 <  1:326291/140‑1 (MQ=255)                                                                                           tAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAggg                                     <  2:99051/140‑1 (MQ=255)                                                                                                 cGTTCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGa                               >  1:242867/1‑140 (MQ=255)                                                                                                   ttCTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTC‑AAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGacac                            >  2:283096/1‑140 (MQ=255)                                                                                                     cTTGATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACAcc                           <  2:242867/140‑1 (MQ=255)                                                                                                        gATGCACTTTCCATCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCaa                        <  1:393048/140‑1 (MQ=255)                                                                                                                    aTCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCt            <  2:425759/140‑1 (MQ=255)                                                                                                                     tCAACGTCAACAATATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCtt           >  1:149984/1‑140 (MQ=255)                                                                                                                            caacaaTATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGAc    <  1:433142/140‑1 (MQ=255)                                                                                                                            caacaaTATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGAc    <  1:164701/140‑1 (MQ=255)                                                                                                                             aacaatATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTC‑AAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGACat  <  1:283096/140‑1 (MQ=255)                                                                                                                              acaataTCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGACat  <  1:192938/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        AACTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCATGCTGCAAATTAACATGATCGGCGTAACATGACTAAAGTACGTAATTGCGTTCTTGATGCACTTTCCATCAACGTCAACAACATCATTAGCTTGGTCGTGGGTACTTTCCCTCAGGACCCGACAGTGTCAAAAACGGCTGTCATCCTAACCATTTTAACAGCAACATAACAGGCTAAGAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGACAT  >  NC_000913/4640326‑4640590

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