Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867,196	A→C	D210A (GAT→GCT)	iaaA →	Isoaspartyl peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867,196	0	A	C	100.0%	251.4 / NA	67	D210A (GAT→GCT)	iaaA	Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (42/25);  total (42/25)

AACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGA  >  NC_000913/867050‑867340                                                                                                                                                   |                                                                                                                                                 aaCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGt                                                                                                                                              >  1:178571/1‑151 (MQ=255)  aCCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTc                                                                                                                                             >  2:434993/1‑151 (MQ=255)     gTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCt                                                                                                                                          >  2:769476/1‑151 (MQ=255)         tCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGt                                                                                                                                      >  1:11473/1‑151 (MQ=255)          cGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAg                                                                                                                                       <  2:245844/149‑1 (MQ=255)          cGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTg                                                                                                                                     >  2:535342/1‑151 (MQ=255)           gACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTgg                                                                                                                                    >  1:607681/1‑151 (MQ=255)            aCCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTggg                                                                                                                                   >  1:648028/1‑151 (MQ=255)                 aGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCg                                                                                                                              >  2:908303/1‑151 (MQ=255)                    ggTGCGCCACCGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGAt                                                                                                                           <  1:577361/151‑1 (MQ=255)                            aCTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGcc                                                                                                                   <  2:405605/151‑1 (MQ=255)                            aCTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGcc                                                                                                                   >  1:911881/1‑151 (MQ=255)                              tGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCa                                                                                                                  >  2:759993/1‑150 (MQ=255)                              tGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCa                                                                                                                  >  1:758262/1‑150 (MQ=255)                               ggATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAAt                                                                                                                <  1:495371/151‑1 (MQ=255)                                   gAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAAc                                                                                                             >  2:81778/1‑150 (MQ=255)                                   gAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAAc                                                                                                             <  1:115742/150‑1 (MQ=255)                                    aaaaaCAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGc                                                                                                           >  1:198251/1‑151 (MQ=255)                                    aaaaaCAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGc                                                                                                           >  2:785887/1‑151 (MQ=255)                                      aaaCAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGcc                                                                                                          >  1:463339/1‑150 (MQ=255)                                       aaCAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGcc                                                                                                          >  2:420795/1‑149 (MQ=255)                                          aGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAgtgt                                                                                                     >  2:805182/1‑151 (MQ=255)                                              aaTGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCg                                                                                                 <  2:216204/151‑1 (MQ=255)                                                 gggCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGtt                                                                                              >  2:227194/1‑151 (MQ=255)                                                       cGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGt                                                                                        >  1:293571/1‑151 (MQ=255)                                                        gTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGt                                                                                        >  1:213674/1‑150 (MQ=255)                                                        gTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGt                                                                                        >  1:221466/1‑150 (MQ=255)                                                        gTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTa                                                                                       <  1:211923/151‑1 (MQ=255)                                                             ggCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCgg                                                                                   >  2:684933/1‑150 (MQ=255)                                                              gCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCgg                                                                                   >  2:586257/1‑149 (MQ=255)                                                                  tGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACggg                                                                             >  1:573068/1‑151 (MQ=255)                                                                      gTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCg                                                                           >  2:125190/1‑149 (MQ=255)                                                                       ttGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGc                                                                            >  2:456919/1‑147 (MQ=255)                                                                              tAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCa                                                                  <  2:905065/150‑1 (MQ=255)                                                                              tAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCa                                                                  >  1:705031/1‑150 (MQ=255)                                                                              tAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAAGCTTCa                                                                  >  2:465296/1‑150 (MQ=255)                                                                                gACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATcc                                                               <  1:467751/151‑1 (MQ=255)                                                                                  cGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCcgc                                                             <  2:43956/151‑1 (MQ=255)                                                                                   ggCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCcgc                                                             >  1:900797/1‑150 (MQ=255)                                                                                      aaTTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTTTACCGGCACGGGCGAAGACTTTATCCccgcgc                                                         >  2:503537/1‑151 (MQ=255)                                                                                       aTTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCgcgcgc                                                         >  2:670723/1‑150 (MQ=255)                                                                                       aTTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCt                                                        <  1:716896/151‑1 (MQ=255)                                                                                                 aGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCatat                                              >  1:889813/1‑151 (MQ=255)                                                                                                  gCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATg                                             <  2:646153/151‑1 (MQ=255)                                                                                                   ccACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGa                                            <  1:465308/151‑1 (MQ=255)                                                                                                    cACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTTTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCatat                                              >  2:681649/1‑148 (MQ=255)                                                                                                     aCGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGAc                                           >  1:731708/1‑150 (MQ=255)                                                                                                       gTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATc                                        <  1:504625/151‑1 (MQ=255)                                                                                                           acaGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATcgcc                                     >  2:842300/1‑150 (MQ=255)                                                                                                             aGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCcgcg                                  <  1:129752/151‑1 (MQ=255)                                                                                                                  gAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAAt                             <  1:477895/151‑1 (MQ=255)                                                                                                                  gAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAAt                             <  2:258481/151‑1 (MQ=255)                                                                                                                  gAATGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAAt                             <  1:59317/151‑1 (MQ=255)                                                                                                                   aaTGACCAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGACATATGACATCGCCGCGTTAATg                            <  1:374494/151‑1 (MQ=255)                                                                                                                         cAATACATTACCCGGACGCGTTGGCGCTCGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCCCGCGCTGGCGGCCTATGACATCGCCGCGTTACTGGAtt                        >  2:403309/1‑149 (MQ=255)                                                                                                                         cAATAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTAc                      <  2:606041/151‑1 (MQ=255)                                                                                                                            tAAATTACCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTAcgg                    >  1:202916/1‑150 (MQ=255)                                                                                                                               aTTGCCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTAcggc                   <  1:894685/148‑1 (MQ=255)                                                                                                                                  aCCCGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGAtt              <  1:525977/150‑1 (MQ=255)                                                                                                                                    ccGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTAcggcggggtta            >  2:777921/1‑145 (MQ=255)                                                                                                                                    ccGGACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGGtt              >  1:317510/1‑148 (MQ=255)                                                                                                                                       gACGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGGCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGGTTTAggc         >  2:356859/1‑148 (MQ=255)                                                                                                                                        aCGAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTTCGGCGGGttta            >  2:712369/1‑144 (MQ=255)                                                                                                                                          gAGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTcgc     <  1:718999/151‑1 (MQ=255)                                                                                                                                           aGTTGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGGTTAAGTCTcgc     >  1:49891/1‑150 (MQ=255)                                                                                                                                              tGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGa  <  2:821106/150‑1 (MQ=255)                                                                                                                                              tGGCGCTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGa  <  2:49891/150‑1 (MQ=255)                                                                                                                                                gcgcTAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGGTTAAATCTCGCGGa  >  2:781420/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AACCGTCCTCGACCATAGCGGTGCGCCACTGGATGAAAAACAGAAAATGGGCACCGTGGGGGCCGTGGCGTTGGATTTAGACGGCAATTTGGCGGCAGCCACGTCCACAGGCGGAATGACCAATAAATTACCCGGACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGA  >  NC_000913/867050‑867340

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