Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,867,576	G→T	intergenic (‑154/‑152)	ibpA ← / → yidQ	small heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,867,576	0	G	T	100.0%	235.1 / NA	66	intergenic (‑154/‑152)	ibpA/yidQ	small heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (31/35);  total (31/35)

ATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGG  >  NC_000913/3867330‑3867816                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 aTTACTCTGGCTCTGGTTGCTATCTAAGTGGTTAAACACACGGACAAATCCAATAGCAGAACGGTAAAGCGGGAATAAATCAACGTTACGCATAATCAGTAGCTCCTGAACTCAGCGAGAACGTAAGACCTTCCACACTGGACAGGTCAGGAAGCCAGAACACCAATAAGGCGCGACCTCATCGGCTACGTTGTAACAATGGGTCTGGAAATGCATTTTNCAAGCCTTCTCAGCGGACTTTTTTTTTCAGt                                                                                                                                                                                                                                              >  2:362152/1‑251 (MQ=255)   tACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACCCCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGttt                                                                                                                                                                                                                                            >  2:185918/1‑251 (MQ=255)   tACTCTGGCTATGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAAGGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGtt                                                                                                                                                                                                                                             <  2:229923/250‑1 (MQ=255)    tctctggctctggTTGATTTCTAAGTGGTTAAACAAACGTTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGACAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGttt                                                                                                                                                                                                                                            <  2:304283/249‑1 (MQ=255)         ggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACTCCCATCAGGCGCGCCCTCATCGGCTACGATGTAAAAATTGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTc                                                                                                                                                                                                                                      <  2:98422/251‑1 (MQ=255)         ggctctggTTGTTTTCTAAGGGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCGCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGTCTACGATGTAAAAATGGGTCGGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATgg                                                                                                                                                                                                                                        <  1:419222/249‑1 (MQ=255)          gctcttgTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGTTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCt                                                                                                                                                                                                                                     <  1:720856/251‑1 (MQ=255)            tctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTAt                                                                                                                                                                                                                                   >  1:161548/1‑251 (MQ=255)                  tgttTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAGCAAAGTTACGCATAATCAATAGCGCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGTACACCCATCAGGCGGGTCCTCATCGGCTCCGATGTAAAAATGGGTCTGGAATTGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCAt                                                                                                                                                                                                                              <  1:353454/250‑1 (MQ=255)                      ttCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCCAAGTTACGCATAATCAATAGCTACTGAAATCAACGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGAAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGt                                                                                                                                                                                                                           >  1:413785‑M1/1‑248 (MQ=255)                          aaGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtaga                                                                                                                                                                                                                        <  1:663843‑M1/248‑5 (MQ=255)                           aGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCAGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttg                                                                                                                                                                                                                    >  2:381870‑M1/1‑243 (MQ=255)                               gTTAAACCCACGGACAAATCCAATAGCAGAACGGACAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGAGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCAAGAACACCCAACAGGAGCGTCCGCATCGGCTACGATGTAAAAATGGGCATGGAAATGAACTTTTAAAACCTTATCACCCGAATTTTTTTTTCAGTTTATGGTCTATTGCATAgttaatatgcggg                                                                                                                                                                                                                >  2:102640/1‑239 (MQ=255)                               gTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                                >  1:345610‑M1/1‑239 (MQ=255)                                 tAAACAAACGGTCAAAACCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggga                                                                                                                                                                                                               <  2:108625‑M1/250‑14 (MQ=255)                                  aaacaaacGGGCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAATTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTTTGGAAATGAACTTTTCAATCCTTATCTGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagattttcgggat                                                                                                                                                                                                              <  2:640883‑M1/250‑15 (MQ=255)                                  aaaaaaacggtcAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAATACCTGCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTAACAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagagttgcgggatg                                                                                                                                                                                                             <  2:139249‑M1/247‑16 (MQ=255)                                  aaaaaaacGGTCAAATCCAATAGCAGAACGGTATAGCAGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTAAAATCAGCGAGAATGTAAGACCTTCCCAAATGTGCAGGTCAGGTGGTCAGAACACCCACCGGGCGGGTCCTCATCGTTTACGATGTAAAAATGTTTCTGGAAATGATCTTTTCAAGACTTATCAGCTGACTGTTTTTGTCAGTTTTTGTTCTATTGCATAGaaggcttgcgggat                                                                                                                                                                                                              <  2:646151/246‑14 (MQ=255)                                             cAAATCCAATAGCAGATCGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTTGAAATGAACTTTTCAAGCCTTATCAGCGGACGTTTTTGTTCAGTTTATGGTCTATTGCATAGtatatttgcgggatgataattgcgg                                                                                                                                                                                                   <  1:605867‑M1/250‑26 (MQ=255)                                             cAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAACGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGAACTTCCACAATGGACAGGTCAGGAAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGCAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcgg                                                                                                                                                                                                   >  2:605867‑M1/1‑225 (MQ=255)                                                aTCCAATAGCAGAACGGTAAAGCGGGGATAAATCATAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattg                                                                                                                                                                                               >  1:33884‑M1/1‑222 (MQ=255)                                                  ccAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgc                                                                                                                                                                                              <  2:274266‑M1/250‑31 (MQ=255)                                                       aGCAGAACGGAAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaat                                                                                                                                                                                        >  1:546499‑M1/1‑215 (MQ=255)                                                         cAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGCCCTTCCACAATGGACAGGTCATGTAGCCAGAACCCCCATCAGGCGCGTCCTCATCGGCTGCGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCATCGTACTTTTTTTTTCCGGTTATGGTCTATTGCATAGtagattgtcgggatgataattgcggattgcagaaataa                                                                                                                                                                                      <  2:473652‑M1/251‑39 (MQ=255)                                                                gTAAAGCGGGGATAAACCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggc                                                                                                                                                                               <  2:227097‑M1/251‑46 (MQ=255)                                                                        gggATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatg                                                                                                                                                                        >  2:276756‑M1/1‑198 (MQ=255)                                                                          gATAAATCAAAGTTACGCATAATCAATAGCTNCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatc                                                                                                                                                                     <  2:348696‑M1/251‑56 (MQ=255)                                                                               aTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgc                                                                                                                                                                >  1:349301‑M1/1‑191 (MQ=255)                                                                                 cAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcc                                                                                                                                                               <  1:512871‑M1/250‑62 (MQ=255)                                                                                  aaaGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccct                                                                                                                                                             <  1:680952‑M1/251‑64 (MQ=255)                                                                                  aaaGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccct                                                                                                                                                             >  1:483661‑M1/1‑188 (MQ=255)                                                                                        aCGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttct                                                                                                                                                        <  1:188696‑M1/250‑69 (MQ=255)                                                                                        aCGCAGAATCAATAGCCCCTGAAATCAGCGAGAATGTCAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCTTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTGTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttctt                                                                                                                                                       <  2:47900‑M1/251‑70 (MQ=255)                                                                                                      gCTCCTGAAATCAGCGAGAATTTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccg                                                                                                                                         <  2:656557‑M1/251‑84 (MQ=255)                                                                                                          cTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatg                                                                                                                                     >  1:428362‑M1/1‑164 (MQ=255)                                                                                                               aTCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcg                                                                                                                                <  2:669941‑M1/251‑93 (MQ=255)                                                                                                                   gCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCCGAACACCCATAAGGCGCGTCATCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAATTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaacgggcggagatgatctccgcccttttcttatagcatcttgccggatgcggggtga                                                                                                                             >  2:281365‑M1/1‑155 (MQ=255)                                                                                                                       gaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgcc                                                                                                                        <  1:156610‑M1/251‑101 (MQ=255)                                                                                                                          tGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagctttgcgggctgataattgcggcttgcagaaataaaaagggcggagatgctctccgcccttttcttatagcttcttgccggatgcggcgtgaacccctt                                                                                                                      >  2:108288‑M1/1‑148 (MQ=255)                                                                                                                                       acaATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaa                                                                                                         >  1:602245‑M1/1‑135 (MQ=255)                                                                                                                                        caATGGACAGGTCAGGTAGCCGAAACACCCAACAGGCGCGTCCTCAACGTCTTCGATGGAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGGCTTGTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaa                                                                                                       <  2:565747‑M1/251‑118 (MQ=255)                                                                                                                                        caATGGACAGGTCAGGTAGCCAGAACACCCANCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaa                                                                                                       <  2:345610‑M1/251‑118 (MQ=255)                                                                                                                                         aaTGGAAAGGTCAGGAAGCCAGAACACACAGCAGGCGCGTCATCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCGTATCAGCGGACTTTTTTTTTCAGATTATGGTCTATTGCATAGtagatgtgcgggatgataattgcagattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcgtcttggcggatgcggcgtgaacgccttatccggcctacaaaat                                                                                                      <  2:199891‑M1/251‑119 (MQ=255)                                                                                                                                          aTGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatc                                                                                                     <  1:677690‑M1/251‑120 (MQ=255)                                                                                                                                             gACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatg                                                                                                  >  2:406018‑M1/1‑129 (MQ=255)                                                                                                                                               aaggtcaggtaGCCAAAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaa                                                                                                <  2:460662‑M1/250‑125 (MQ=255)                                                                                                                                                             aGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattccatacattg                                                                                  >  1:623527‑M1/1‑113 (MQ=255)                                                                                                                                                                               gTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcct                                                                <  1:616757‑M1/251‑157 (MQ=255)                                                                                                                                                                                        ggCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgt                                                       >  1:361349‑M1/1‑86 (MQ=255)                                                                                                                                                                                               aTGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcacgattttcgtaggcctgataagcgtagcgcat                                                >  1:243518‑M1/1‑79 (MQ=255)                                                                                                                                                                                               aTGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcat                                                <  1:338344‑M1/251‑173 (MQ=255)                                                                                                                                                                                                      aaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcac                                         >  1:144685‑M1/1‑72 (MQ=255)                                                                                                                                                                                                         gggTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCAGAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgct                                      <  2:16079‑M1/251‑183 (MQ=255)                                                                                                                                                                                                          ggTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctc                                     >  1:183876‑M1/1‑68 (MQ=255)                                                                                                                                                                                                           gTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcc                                    >  1:524712‑M1/1‑67 (MQ=255)                                                                                                                                                                                                             cTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggcttgcagaaataaaaagggcggaactgatctccgcccttttcttatagcttcttgccggatgcggcgggaaccccttatccggcctacaaaatcataaaaattccatacattgcaagattttctgaggcctgaaaagcttaacgcatccggcccgctcgc                                   >  2:487202‑M1/1‑65 (MQ=255)                                                                                                                                                                                                                   aTGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggtt                             >  1:339714‑M1/1‑59 (MQ=255)                                                                                                                                                                                                                      aaCTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcg                         <  1:361655‑M1/251‑196 (MQ=255)                                                                                                                                                                                                                           ttCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgctaattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccatt                    >  1:467815‑M1/1‑51 (MQ=255)                                                                                                                                                                                                                           ttCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcctgaaaattcaataccttgcaagattttggtaggcctgataagcgaagcgcataaggcacgcgcgcctggttagcgccatt                    >  2:36773‑M1/1‑51 (MQ=255)                                                                                                                                                                                                                           ttCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccat                     <  1:253127‑M1/250‑200 (MQ=255)                                                                                                                                                                                                                            tCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgaccttttcttatagattcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcccgctcgcatggttagcgccatt                    >  1:415779‑M1/1‑50 (MQ=255)                                                                                                                                                                                                                                  cTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatccccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatc             >  1:341952‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                                   ttATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggctgcggcgtgaacgccttatccggcctaccacatcatgaacattcactacattgcaagattttcgtaggcctgatacgcgtcacgcatcaggcacgataccatgtttagcgcccttcaatatcg            >  2:82118‑M1/1‑43 (MQ=255)                                                                                                                                                                                                                                     aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcga           <  1:647673‑M1/250‑210 (MQ=255)                                                                                                                                                                                                                                             aCTTTTTTTTTCAGTTTATGGTCTATTGCATNGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgccg  <  2:349301‑M1/251‑219 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGG  >  NC_000913/3867330‑3867816

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