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%	191.1 / NA	53	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 (26/27);  total (26/27)

ACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGGAAA  >  NC_000913/3867333‑3867820                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     aCTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGCCTTTTTTTTTCAGttt                                                                                                                                                                                                                                                >  2:78137/1‑250 (MQ=255) aCTCGGGCTCTGGTTGTTTTCTAGTTGGTTAATCAAACGGTCAAATCCAATAGCAGAACGTTAACGCGGGGATAAATCAAAGTTACGCTTAATCAATAGCCCCTGAAATCAGCGATAATGTAATACCTTCCACAATGGACAGGTCAGGTAGCCAGAACTCCCATCAGGCGCGTCCTCATCGGCTACGATGCAAAAATGGGTCGGGAATTGAACTTTTCAAGCCTTATCAGCGGACTTGTTTTTTCAGttt                                                                                                                                                                                                                                                <  2:361412/250‑1 (MQ=255)    ctgtccgtggttgtgttataacgtgTTAANCAAACGGGCAAAACCAACACTGGAACCGTAAAGCGGGGATAAAACAAAGTTAGGCATAATCAATAGCGCCTGAAATCAGCGAGAATGTATGGCCATCCACAATGGACAGGTCAGGTAGCGAGGACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGTGTCTGGAAACGAACTTTTCAAGCCTTATCAGCGGACGTTTTTTTTCAGTTTAt                                                                                                                                                                                                                                              <  2:226821/225‑1 (MQ=255)                      tAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTTTGGTCTATTGCCTAGtaggtt                                                                                                                                                                                                                          >  2:249441‑M1/1‑245 (MQ=255)                       gagTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGATCCGATGTCAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGTTATATTGCATAGtagattt                                                                                                                                                                                                                         >  2:431492‑M1/2‑244 (MQ=255)                       aaGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagattt                                                                                                                                                                                                                         <  1:442968‑M1/251‑8 (MQ=255)                            ggtaaACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGTCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                                    <  2:459045‑M1/249‑13 (MQ=255)                               aaacaaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatg                                                                                                                                                                                                                 >  1:224766‑M1/1‑236 (MQ=255)                                 acaaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgat                                                                                                                                                                                                               <  1:4954‑M1/251‑18 (MQ=255)                                     acGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataatt                                                                                                                                                                                                           >  1:249620‑M1/1‑230 (MQ=255)                                            aaTCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGGGAATGTAAGACCTTCCACAATGGCAAGGTCAGGTAGCCAGAATACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAATTGGGTCTGGAAATGAAATTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCAGAGtagatttgcgggatgataattgcggatt                                                                                                                                                                                                    <  2:242078‑M1/251‑29 (MQ=255)                                               ccAATAGCAGAACGGTAAAGCGGGGATAAATGAAAGTTACGCAGAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATTAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgc                                                                                                                                                                                                  <  2:431108‑M1/250‑31 (MQ=255)                                                cAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCGCCTGAAATCAGGGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGGCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcag                                                                                                                                                                                                <  1:21538‑M1/251‑33 (MQ=255)                                                  agagccgaacggtAAAGCGGGGATAAATCAAAGGTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGAAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACAACCATCAGGCGCGTCCTCATCGGCTACGATGGAAAAATGGGACTGGAAATGAACTTTTCAAGCCTTATCATCGGACTTTTTTTTTAAGTTTATGGTCTATTGCTTAGtagatttgcgggatgataattgcggattgcaaa                                                                                                                                                                                               >  2:38789‑M1/7‑217 (MQ=255)                                                       aGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataa                                                                                                                                                                                          >  1:271733‑M1/1‑212 (MQ=255)                                                                  gCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGGGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcgg                                                                                                                                                                                 >  1:254279‑M1/1‑201 (MQ=255)                                                                       gATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatc                                                                                                                                                                         <  1:355938‑M1/251‑56 (MQ=255)                                                                            aTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggggtgataattgcggattgcagaaataaaaagggcggagaggatctccgc                                                                                                                                                                    >  2:270372‑M1/1‑191 (MQ=255)                                                                                   ttACGCATAATTAATAGCTACTGAAATCAGCGAGAATGTAAGAACTTCAACAATTGACAGGTCAGGACGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGCACTTTTCAAGCCTTATCCGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatgtgcggggtgataattgcggattgcagaaataaaaagggcggagctgataccagcccttttc                                                                                                                                                             >  2:195899‑M1/1‑184 (MQ=255)                                                                                          tAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagc                                                                                                                                                      >  1:446860‑M1/1‑177 (MQ=255)                                                                                           aaTCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagct                                                                                                                                                     >  1:99013‑M1/1‑176 (MQ=255)                                                                                              cAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATATAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttct                                                                                                                                                  <  2:346140‑M1/251‑79 (MQ=255)                                                                                                 tAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGCCAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGGTGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagctttgcgggatgataattgcggattgcagaaataaaaaggtcggcgatgatctccgcccttttcttatagcttcttgc                                                                                                                                               >  2:303562‑M1/1‑170 (MQ=255)                                                                                                         gAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcg                                                                                                                                       <  2:454400‑M1/251‑90 (MQ=255)                                                                                                          aaaTCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcgg                                                                                                                                      >  1:250476‑M1/1‑161 (MQ=255)                                                                                                               aGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtg                                                                                                                                  >  1:415739‑M1/1‑156 (MQ=255)                                                                                                                  gagaATGTAAGCCCCTCCACAATGGACAGGTCAGGTAGCCAGAACCCCCAGCAGGCGCGTCCTCATCGGCTACGATGAAAAAATGGGTCTGGAAATGAACTTTACAAGCCTTATCAGCAGTCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacg                                                                                                                              <  1:222992‑M1/251‑99 (MQ=255)                                                                                                                  gagaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                      >  2:57468‑M1/1‑153 (MQ=255)                                                                                                                  gagaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                      <  1:57468‑M1/211‑59 (MQ=255)                                                                                                                   agaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgc                                                                                                                             <  2:101276‑M1/251‑100 (MQ=255)                                                                                                                   agaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgc                                                                                                                             <  1:283090‑M1/251‑100 (MQ=255)                                                                                                                        gTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgagctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttat                                                                                                                        <  2:406576‑M1/251‑105 (MQ=255)                                                                                                                               cTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaaccccttatccggcct                                                                                                                 >  1:393532‑M1/1‑140 (MQ=255)                                                                                                                                   cacaATGGACAGGGCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaa                                                                                                             <  1:204565‑M1/251‑116 (MQ=255)                                                                                                                                              ggtcaggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                      <  2:417564‑M1/183‑59 (MQ=255)                                                                                                                                              ggtcaggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                                      >  1:417564‑M1/1‑125 (MQ=255)                                                                                                                                               gtcaggtAGCCAGACCACCCAGCAGTCGCTTCCTCATCGTCTCCGATGTAAAAATGTGTCTGGAAATGAACTTTTCAAGCCTTATCCGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatcgccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaat                                                                                                 <  2:45877‑M1/251‑128 (MQ=255)                                                                                                                                                          aGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattg                                                                                      >  2:1871‑M1/1‑113 (MQ=255)                                                                                                                                                              cacCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcgggttgcagaaataaaaagggcggagctgatctccgcccttttcttatagcttcttgccggatgcggagtgaacgccttatccggcctacaaaatcaagaaaattcaatacattgca                                                                                    >  2:144573‑M1/1‑109 (MQ=255)                                                                                                                                                                cccATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagat                                                                                <  2:110655‑M1/251‑145 (MQ=255)                                                                                                                                                                 ccATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagatt                                                                               <  1:242331‑M1/251‑146 (MQ=255)                                                                                                                                                                         cgcgTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtagg                                                                       >  1:351590‑M1/1‑98 (MQ=255)                                                                                                                                                                                                   aaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggattcggcgtgaacgccttatccggcctacaaaatcatgaaaattccatacattgcaagattttcgtaggcctgataagcgtcgcgcatcaggcac                                             <  1:257332‑M1/251‑180 (MQ=255)                                                                                                                                                                                                   aaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTAAGTTTATGGTCTATTGCATAGtagattttcgggatgttaattgcggtttgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccgaatgcggcgtgaaccccttatccggcatccaaaataatgaaatttcaataacttgcaaaattttcctgggcctgataagcgtagcgcaacaggcac                                             >  2:296953‑M1/1‑72 (MQ=255)                                                                                                                                                                                                    aTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaaaaatgaaaattcaatacattgcaagaatttcgtaggcctgataagcgccgcgcatcaggcacc                                            >  1:277985‑M1/1‑71 (MQ=255)                                                                                                                                                                                                               aaTGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgctaagcgtagcgcatcaggcacgctagcatggtt                                 >  2:76649‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                           aaGCCTTATCCGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggctgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttctagcttcttgccggaggcggcgtgaacgccttatccggtctacaaaatcatgaaaattcaatacaatgcaagattttcgtaggcctgataagcgtaacgcatccggcacgcacgcatggttagcgccattaa                      >  2:411771‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                            aGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATGGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatgtttagcgccattaaat                    <  1:131456‑M1/251‑205 (MQ=255)                                                                                                                                                                                                                                  aTCAGCGTACTTTTTTTTTCAGTTTATGGTCNATTGCATAGgagattggcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggaggcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgat              <  2:207972‑M1/251‑211 (MQ=255)                                                                                                                                                                                                                                      gCGGCCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattc          <  2:254279‑M1/251‑215 (MQ=255)                                                                                                                                                                                                                                       cGGACTTTTTTTTTCAGTTTATGGTCTATTGNATAGtagatttgcgggatggtaattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggccttataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcg         <  2:218796‑M1/251‑216 (MQ=255)                                                                                                                                                                                                                                        ggACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgc        <  2:415739‑M1/251‑217 (MQ=255)                                                                                                                                                                                                                                              ttttttTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggcttgcagaaataaaaagggcggagatgatctccgcacttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgttgacctgataagcgtagcgcatcaggcacgctcgcatggttagcgcccttaaatatcgatattagccgcttt  >  1:54330‑M1/1‑29 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     ACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGGAAA  >  NC_000913/3867333‑3867820

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: ‑