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%	199.2 / NA	55	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 (32/23);  total (32/23)

ACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAG  >  NC_000913/3867333‑3867815                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                aCTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCCCAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGttt                                                                                                                                                                                                                                           <  1:96171/250‑1 (MQ=255)      ggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGTTAGCCAGAACACCCATCAGGCTCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGt                                                                                                                                                                                                                                      <  1:530879/250‑1 (MQ=255)      ggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGTCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTc                                                                                                                                                                                                                                     <  1:441714/251‑1 (MQ=255)        ctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCtt                                                                                                                                                                                                                                   >  1:311927/1‑250 (MQ=255)           tggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCTTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTAtt                                                                                                                                                                                                                                 <  2:214801/250‑1 (MQ=255)            ggTTGTTTTCTAAGTGGTTAAACAAACGGTAAAAACCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTg                                                                                                                                                                                                                                <  2:94917/250‑1 (MQ=255)                         gTGGTTAAAAAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTAAGCATAATAAATAGCTCCTGAAATCAGCGAGAATGAAAGACCTTCCACAATGGACAGGTCCGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTAAAGCCTTATCAGCGGACTTTTTTTTTCAGGTTATGGTCTATTGCATAGtagatttgc                                                                                                                                                                                                                  >  2:455371‑M1/1‑242 (MQ=255)                            gTTAAACAAACGGTCAAATCCAATAGCAGAACGGAAAAGAGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGaagatttgcggg                                                                                                                                                                                                               >  1:428906‑M1/1‑239 (MQ=255)                                 acaaccGGCCAAATCCAAAAACATAACGGTCAAGCGGGGCTAAATCACAGTTACGCATAATCAATAGCACATGAAATCAGCGAAAAGGTAAGACCTTCCACAAAGGACAGGTCAGGAAGCCAGAACACCCATCAGGCGCGTCCTCATCGACTACGGTGTAAAAATGGGTCTGGAAATGAACTTTACAAGCCTTAACAGCGGACTTTATTTTTCAGTTTATGGTATATTGCATAgtagctttgcggggtgat                                                                                                                                                                                                          >  2:493844/1‑234 (MQ=255)                                 acaaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgat                                                                                                                                                                                                          >  1:164559‑M1/1‑234 (MQ=255)                                   aaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTAt                                                                                                                                                                                                                                         >  1:419573/1‑218 (MQ=255)                                   aaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgata                                                                                                                                                                                                         >  1:210195‑M1/1‑232 (MQ=255)                                   aaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTGCCGCATAATCAATAGCTACTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGTTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCAGCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataa                                                                                                                                                                                                        <  2:417904‑M1/251‑20 (MQ=255)                                      cGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattg                                                                                                                                                                                                     >  1:177218‑M1/1‑229 (MQ=255)                                          cAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcgg                                                                                                                                                                                                  <  1:90792‑M1/250‑26 (MQ=255)                                               ccAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCATCGAGAATGTAAGACCTTCCCCAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGTTACGATGTAAAATTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgc                                                                                                                                                                                             <  2:56227‑M1/250‑31 (MQ=255)                                                  aTAGCAGAACGGGAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATAAGCGAGAATGGAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaa                                                                                                                                                                                         <  2:474651‑M1/251‑35 (MQ=255)                                                    agtagaaCGTTAAAGCGGGGATAAATCAAATGTACGCCTAATCAATAGCTCCTGAAAACAGCGAGAAGGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACCCCCATCAGGCGCGTCCTCATCGGCTAAGATGTAAAAAGGGGTCTGGAAATGAACTTGTCAAGCCTTATCAACGGACTTTTTTTTTCAGTTTATGGGCTAATGCATAGtagatttgcgtgatgataagtgcggattgcagaaat                                                                                                                                                                                       <  2:23043‑M1/248‑37 (MQ=255)                                                      cAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaata                                                                                                                                                                                      >  1:230898‑M1/1‑213 (MQ=255)                                                             gTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCAAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggc                                                                                                                                                                              <  1:136650‑M1/251‑46 (MQ=255)                                                                     gggATAAATCAAAGTTACGCATAATCAATAGCTCCTGACATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGCACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagctttgcgggatgataattgcggattgcagaaatacaaatggcggaacttt                                                                                                                                                                      >  1:349805‑M1/1‑198 (MQ=255)                                                                           aaTCACAGTTACGCATAATCAATAGCTCCTGAAATCAGCGGGAATGTAAGACCTTACACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCACGTCATCATCGGCTCAGATGTAAAAATGGGCCCGAAAATGAACTTTTCAAGCCTTATCAGCGGAATTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggcttgcagaaataacaagggcggggatgatacccg                                                                                                                                                                >  2:207085‑M1/1‑192 (MQ=255)                                                                            aTCAACGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGCCAGGTCAGGTAGCCAGACCACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgaaaattgcggattgcagaaataaaaagggcggagatgatctccgc                                                                                                                                                               >  2:329028‑M1/1‑191 (MQ=255)                                                                               aaaGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccct                                                                                                                                                            >  1:308286‑M1/1‑188 (MQ=255)                                                                                 aGTTACGCATAAAAAATAGCTCCTGAAAGCAGCGATAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttt                                                                                                                                                          <  2:164559‑M1/251‑66 (MQ=255)                                                                                    tACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                               >  1:323973‑M1/1‑183 (MQ=255)                                                                                    tACGCAAAAACAATAGCTCCGGAAATCAGCGAGAATGTAAGACCTGCCACAAGGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAAGGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                               <  2:323973‑M1/195‑13 (MQ=255)                                                                                     aCGCATAATCAATAGCTCCTGAAATCAGCGAGAATGGAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagctgatctccgcccttttct                                                                                                                                                       <  2:491121‑M1/250‑69 (MQ=255)                                                                                        cATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttat                                                                                                                                                    >  2:419813‑M1/1‑179 (MQ=255)                                                                                                  aGCTCCTGAAAGCAGCGAGAATGTAAGCCCTTCCCCAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgcc                                                                                                                                         <  2:330634‑M1/251‑83 (MQ=255)                                                                                                    cTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGAAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccgg                                                                                                                                       >  1:326407‑M1/1‑167 (MQ=255)                                                                                                                 cGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCATCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgattattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaac                                                                                                                          >  2:563287‑M1/1‑154 (MQ=255)                                                                                                                     aaTGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGCCTGGAACTGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggctgataattgcggattgcagaaataaaaagggcggagatgctctccgcccttttcttatagcttctttccggattcggcgtgaccgccc                                                                                                                      >  2:447746‑M1/1‑150 (MQ=255)                                                                                                                        gTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGCTGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACCTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttattatagcttcttgccggatgcggcgtgaaccccttat                                                                                                                   >  2:195529‑M1/1‑147 (MQ=255)                                                                                                                        gTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgcacgccttat                                                                                                                   >  2:47886‑M1/1‑147 (MQ=255)                                                                                                                          aaGACCTTCCACAAGGGACGGGTAAGGTAGCCAGAACACCCTTCACGCCCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCATGTTATGGTCTATTTCATAGtagatttgctggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggattcggcgttaacgccttatc                                                                                                                  <  2:487227‑M1/250‑106 (MQ=255)                                                                                                                                      aaTGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaat                                                                                                     <  1:207483‑M1/251‑119 (MQ=255)                                                                                                                                       aTGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatc                                                                                                    >  1:568834‑M1/1‑132 (MQ=255)                                                                                                                                                   ggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattca                                                                                         >  1:134700‑M1/1‑120 (MQ=255)                                                                                                                                                          aGAACACCCATCAGGCGAGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTATATGGTCTATTGCATAGtagatttgcgggatgataattgcggggtgcagaaatacaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggagtgaaccccttttccggccaacaaaaacatgacacttcaatacattt                                                                                 >  2:558123‑M1/1‑113 (MQ=255)                                                                                                                                                                      aGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCCGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggggatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgt                                                                     >  1:474968‑M1/1‑101 (MQ=255)                                                                                                                                                                            gTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgaggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtacgcct                                                               >  1:251180‑M1/1‑95 (MQ=255)                                                                                                                                                                              ccTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctg                                                              <  1:330203‑M1/250‑158 (MQ=255)                                                                                                                                                                                   tCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCGGTTTATGGTCTATTGCATAGtagctttgcgggatgataattgcggattgcagaaataaaaagggcggcgatgatctacgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgaaagattttcgtcagcctgata                                                           >  1:325801‑M1/1‑88 (MQ=255)                                                                                                                                                                                     ggCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtggacgccttatccggcctacaaaatcatgaaaattcaatacattgccagattttcgtgtgcctgataagcgt                                                      >  2:206745‑M1/1‑86 (MQ=255)                                                                                                                                                                                         aCGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTAGCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagc                                                   <  1:147206‑M1/250‑169 (MQ=255)                                                                                                                                                                                          cGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgc                                                 >  1:85158‑M1/1‑81 (MQ=255)                                                                                                                                                                                                   aaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcac                                        <  1:439839‑M1/251‑180 (MQ=255)                                                                                                                                                                                                       ggTCTGGAAATTAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctc                                    <  1:223304‑M1/251‑184 (MQ=255)                                                                                                                                                                                                                  gAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCTGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgcgggatgcggcgtgaacgccttatccgtcctacaaaatcatgaaaattcaatacagtgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttag                          <  1:28573‑M1/250‑194 (MQ=255)                                                                                                                                                                                                                       tttCAAGCCTTATCATCGGTCTGTTTTTTTCAGTGTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatcgcagaaagaaaaagagcggagatgatctccgcccttttcttatagcttcgtgccggatgcggcgtgaacgccttatccggcttacaaaatcatgaaaattcaattcatttcaagattttcgtaggcctgctaagcgtagcgcatcaggcacgccggcatggttcgcgcc                      <  2:428906‑M1/249‑198 (MQ=255)                                                                                                                                                                                                                           aaGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaa                 >  1:264318‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                  aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagacataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatcaggcctacaaaatcatgaaaattcaatacattgcaagattttcgtatgcctgataagcgtagcgcatcaggcacgctcgcctggttaggcccattaaatatcgat         >  1:432740‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                                  aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgat         >  1:257594‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                                         gACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggctgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgaaggcctgataagcgtagcgcatcaggcccgctcgcatggttagcgccattaaatatcgatattcgcc  >  1:268910‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                ACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAG  >  NC_000913/3867333‑3867815

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