Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867,112	C→G	A182G (GCC→GGC)	iaaA →	Isoaspartyl peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867,112	0	C	G	100.0%	169.7 / NA	47	A182G (GCC→GGC)	iaaA	Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (33/14);  total (33/14)

CTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTAC  >  NC_000913/866965‑867256                                                                                                                                                    |                                                                                                                                                 cTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCg                                                                                                                                                >  1:347435/1‑150 (MQ=255)      ggCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGggcgtggcgt                                                                                                                                          <  1:414087/151‑1 (MQ=255)               cgcgTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAg                                                                                                                                 <  2:139834/151‑1 (MQ=255)                      cGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCa                                                                                                                           >  2:342463/1‑150 (MQ=255)                      cGCCGGAGATTTTCTCCACATCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCa                                                                                                                           <  1:428805/150‑1 (MQ=255)                        ccGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAAtt                                                                                                                        >  2:398931/1‑151 (MQ=255)                        ccGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAAtt                                                                                                                        <  1:372951/151‑1 (MQ=255)                          ggAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTg                                                                                                                      <  2:51279/151‑1 (MQ=255)                           gagaTTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTgg                                                                                                                     <  1:91071/151‑1 (MQ=255)                                tttCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTggcggc                                                                                                                 >  2:129427/1‑150 (MQ=255)                                 ttCTCCACGTCTTTGCGCTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTggcggc                                                                                                                 >  2:410628/1‑149 (MQ=255)                                         gTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTcc                                                                                                       >  2:60608/1‑151 (MQ=255)                                               gCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGc                                                                                                 >  1:348681/1‑151 (MQ=255)                                                  ttATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCgg                                                                                               >  1:369164/1‑150 (MQ=255)                                                   tATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGa                                                                                              >  2:556655/1‑150 (MQ=255)                                                    aTGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAAt                                                                                            >  2:344372/1‑151 (MQ=255)                                                          aacTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCa                                                                                       >  2:375843/1‑150 (MQ=255)                                                          aacTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCa                                                                                       >  2:288657/1‑150 (MQ=255)                                                             tactGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATa                                                                                    >  1:432902/1‑150 (MQ=255)                                                               ctGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAAt                                                                                 >  2:275268/1‑151 (MQ=255)                                                                 ggCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAAtt                                                                                >  1:536664/1‑150 (MQ=255)                                                                    aGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTAcgc                                                                            <  2:142542/151‑3 (MQ=255)                                                                        cgcAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGAc                                                                        <  1:177071/151‑1 (MQ=255)                                                                         gcAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACg                                                                       >  1:380791/1‑151 (MQ=255)                                                                           aaggaaggGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACg                                                                       >  2:183380/1‑149 (MQ=255)                                                                           aaggaaggGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACg                                                                       >  1:554836/1‑149 (MQ=255)                                                                             ggaaggGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGtt                                                                   <  1:313522/151‑1 (MQ=255)                                                                                    gCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATa                                                            <  1:570828/151‑1 (MQ=255)                                                                                     cAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGAt                                                             >  1:63446/1‑149 (MQ=255)                                                                                          gTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCt                                                      >  1:358359/1‑151 (MQ=255)                                                                                             cTCGACCATAGCGGTGCGCCATTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCtt                                                     >  2:589764/1‑149 (MQ=255)                                                                                                 aCCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTggg                                               >  1:239810/1‑151 (MQ=255)                                                                                                        cGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGa                                        <  2:259975/151‑1 (MQ=255)                                                                                                         ggTGCGCCACTGGGTGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGAt                                       <  1:505625/151‑1 (MQ=255)                                                                                                           tGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGc                                     >  2:530910/1‑151 (MQ=255)                                                                                                              gcCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTAc                                  >  1:127512/1‑151 (MQ=255)                                                                                                                 aCTGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGcc                               >  1:612796/1‑151 (MQ=255)                                                                                                                   tGGATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCa                              >  2:62135/1‑150 (MQ=255)                                                                                                                    ggATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAAt                            <  1:239456/151‑1 (MQ=255)                                                                                                                     gATGAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAAt                            >  1:573213/1‑150 (MQ=255)                                                                                                                        gAAAAACAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACg                        >  1:12104/1‑151 (MQ=255)                                                                                                                            aaCAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGcc                      >  1:104789/1‑149 (MQ=255)                                                                                                                              cAGAAAATGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAgtg                  >  1:251378/1‑151 (MQ=255)                                                                                                                                 aaaaTGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTgg               >  1:577450/1‑151 (MQ=255)                                                                                                                                 aaaaTGGGCACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGAGCCGGATGCTACGCCAATAACGCCAGTGTgg               >  1:267836/1‑151 (MQ=255)                                                                                                                                         cACCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCt       <  2:463599/151‑1 (MQ=255)                                                                                                                                          aCCGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCtt      >  1:603485/1‑151 (MQ=255)                                                                                                                                           ccGTGGGGGGCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTg     <  2:502027/151‑1 (MQ=255)                                                                                                                                               ggggggCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTAc  <  2:127509/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CTCGTGGCATGGAGCGCGTCTCGCCGGAGATTTTCTCCACGTCTTTGCGTTATGAACAACTACTGGCAGCGCGCAAGGAAGGGGCAACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTAC  >  NC_000913/866965‑867256

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