Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,039	0	T	G	39.0%	113.2 / 10.9	61	N108T (AAC→ACC)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base G (23/0);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 2.03e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677890‑2678187                                                                                                                                                         |                                                                                                                                                     ccaGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGtt                                                                                                                                                     <  1:437768/151‑1 (MQ=255)       gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGc                                                                                                                                               <  2:983077/151‑1 (MQ=255)         ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGGTTTCCGccc                                                                                                                                             >  2:288799/1‑150 (MQ=255)             gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTTACCGAGGCAGA‑GGGTTTTCGCCACTGt                                                                                                                                         >  2:886034/1‑151 (MQ=255)             gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTg                                                                                                                                          <  1:805068/150‑1 (MQ=255)               gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCa                                                                                                                                       <  1:313376/151‑1 (MQ=255)                 cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑GCCCCGACCGGGGGCGA‑GGGGTTTCGCCCCTGTCAg                                                                                                                                      >  1:100962/1‑150 (MQ=255)                 cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTGCCCCAGGCAGC‑GGGTTTTCGcccc                                                                                                                                            >  1:671843/1‑142 (MQ=255)                 cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTCACCGCGGCAGG‑GGGTTTTCGCCCCTCTCCGc                                                                                                                                     >  1:742483/1‑151 (MQ=255)                 cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCCGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAg                                                                                                                                      >  1:911049/1‑150 (MQ=255)                   ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGCGG‑CCCCT‑ACAGAGGGGGA‑GGGTTTTCGCCCCTGTCCGccc                                                                                                                                   >  1:406193/1‑151 (MQ=255)                      ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑GCCCTGGCCGAGGGAGA‑GGGTTTTCGCCACTGTCAGCCCCTc                                                                                                                                >  2:192735/1‑151 (MQ=255)                       tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGGTTTCCGCCCCTGTCAGCCCTTcc                                                                                                                               >  2:1025044/1‑151 (MQ=255)                          cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTccacc                                                                                                                            <  2:68118/151‑1 (MQ=255)                           gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTccacca                                                                                                                           <  2:725130/151‑1 (MQ=255)                             tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTTCCCCCCCGTCAGCCCCTccccc                                                                                                                            >  2:628052/1‑148 (MQ=255)                                 aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTa                                                                                                                     <  1:838588/151‑1 (MQ=255)                                         ttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTg                                                                                                              <  1:900196/150‑1 (MQ=255)                                         ttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACAGACCGAGGGAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCCGACTTTACTcc                                                                                                                >  1:900727/1‑148 (MQ=255)                                          tCAGGACTACCGTCAACGCCCAGCGCCTGGAAGCACCGTTGGGCGTCCTGCGCCCTTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACGGTCAGCCCTTCCACCAGACTTAACTCCTGa                                                                                                             <  1:862111/150‑1 (MQ=255)                                          tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGa                                                                                                             <  2:325985/150‑1 (MQ=255)                                          tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGa                                                                                                             <  1:720302/150‑1 (MQ=255)                                           cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACCGACCCAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCCTCCCCCAAACTTTACTCCTggt                                                                                                            >  2:902284/1‑148 (MQ=255)                                             ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACCCACCGACGCAGA‑GGTTTTCCGCCCCCGTCAGCCCTTCCCCCCAACCTAACCCCCGATaa                                                                                                          >  2:568306/1‑150 (MQ=255)                                                 tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTGACCGAGGGAGA‑GGGTTTTCGCCACCGTCAGCCCCTCCCCCCAACTTAACTCCCTATAAAcc                                                                                                       >  2:508609/1‑149 (MQ=255)                                                  aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgc                                                                                                     <  1:533130/150‑1 (MQ=255)                                                     gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTGACCCAGGGCGA‑GGGTTTTCGCCCCCTGCAGCCCTTCCACCCGACTTAAATCCCGATAAAAcccgcg                                                                                                  >  2:291226/1‑150 (MQ=255)                                                      tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑GCACTGACCGAGGGAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCCCCACACTTTACTCCTGATaaaacccccgc                                                                                                 >  2:350665/1‑142 (MQ=255)                                                       cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑GCCCTGACCGAGGGAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCCCACTTAACTCCTTATAAAACCCGCGCAc                                                                                               >  2:928657/1‑151 (MQ=255)                                                         aCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGGTTTCCGCCCCTGTCAGCCCTTCCCCCCGACTTAAATCCTTATAAAACGCGCGCAccc                                                                                             >  2:909466/1‑151 (MQ=255)                                                          cGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GGCCCCCGGACCG‑GGCAGG‑GGGTTTTCGCCCCCGTCCGCCCCTCCCCCCAACCTAAATCCCTATAAAACGCGCGCAcccc                                                                                            >  1:509890/1‑151 (MQ=255)                                                          cGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAcccc                                                                                            <  2:367693/151‑1 (MQ=255)                                                          cGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAcccc                                                                                            <  1:714637/151‑1 (MQ=255)                                                             ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                         <  2:523147/151‑1 (MQ=255)                                                              cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGt                                                                                        <  1:1011432/151‑1 (MQ=255)                                                                  gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACCGACCGCGGCCGG‑GGGTTTTCCCCACTGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAACACCGCGCACCCCCCGtt                                                                                       >  2:358329/1‑148 (MQ=255)                                                                      gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                  <  2:564772/149‑1 (MQ=255)                                                                      gTAAGCAACGTTGGGCATCCGGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCgg                                                                                 <  1:841326/150‑1 (MQ=255)                                                                                 tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTGACCGAGGCCGA‑GGGTTTCCGCCCCTGTCAGCCCTTCCCCCAGACCTAACTCCTGATAAACCCCGCGCCCCCCCCGTTCAGCGGCACGGCGAg                                                                        >  2:996009/1‑148 (MQ=255)                                                                                      aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACCGAACGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTTCACCAGACTTAAATCCTTATAAAACCCGCGCACCCCCCCTTCATCGGCACGGCGAGGCa                                                                     >  2:686318/1‑146 (MQ=255)                                                                                          tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGGAGA‑GGGTTTCCGCCCCTGTCAGCCCCTCCCCCCGACTTTACTCCTGATAAAACGCGCGCCCCCCCCGTTCAGCGGCACGGCGAGGCAGCGTAGcct                                                            >  2:574555/1‑149 (MQ=255)                                                                                          tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACCGACCCCGGGCGA‑GGGTTTTCGCCCCCGTCAGCCCCTCCCCCAGAACTAACTCCTCATAAACCCCCCGCACCCCCCGTTCAGCGGCACGGCGGGTCAGCGTAGCtt                                                            >  1:404266/1‑151 (MQ=255)                                                                                             gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑GCACCGTCCGCGGCAGA‑GGCGTTTCGCCCCTTTCAGCCCTTCCCCCCAACCTTACTCCCGGTAAAACGCGCCCACCCCCAGCTCAGCGGCACGGGGCGCCAGccgagc                                                              >  1:254790/1‑141 (MQ=255)                                                                                             gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCg                                                          <  2:605438/150‑1 (MQ=255)                                                                                               cATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑GCACTGACCGAGGGAGA‑GGTTTTCCCCCCCTGTCAGCCCTTCCCCCAGACTTAAATCCTGATAAAACGCGCGCACCCCCAGGTCAGCGGCACGGCGAGTCCGc                                                                   >  2:360785/1‑139 (MQ=255)                                                                                                   tGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTTACCGCGGGAGA‑GGGTTTTCCCCACTGTCAGCCCTTCCACCCGAATTAACTCCTGATAAACCGCGCGCACCCCCCGGTCCGCGGCACGGCGAGTCAGCGGAGCTTCGTCAcc                                                     >  1:210653/1‑149 (MQ=255)                                                                                                   tGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCCCTGACCGGAGGAGA‑GGATTTTCGCCCCTTTTTGCCCTTCCACCACACTTTACTTCTTATATACCGCGCGCACCCCCCCTTTTCCGGCACGGCGAGGCTGCGTACCTTCGCCAcc                                                     >  1:502939/1‑149 (MQ=255)                                                                                                           atcaatcaTGCCGTTGCGGCGG‑GG‑GCCCTGACCCAGGCCGAGGGTTTTCCCCCCCT‑TCAGCCCTTCCACCCGACTTAACTCCTGATAAAACCCGCGCACCCCCAGTTCAGCGGCACGGCGAGGCAGCGTAGCCTCGTCACCTTCCAGGcg                                            >  2:596905/1‑150 (MQ=255)                                                                                                                       gTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATc                               <  2:388419/151‑1 (MQ=255)                                                                                                                       gTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATc                               <  2:699369/151‑1 (MQ=255)                                                                                                                                g‑gg‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                       >  2:544940/1‑150 (MQ=255)                                                                                                                                     ccACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcc                    >  1:933023/1‑150 (MQ=255)                                                                                                                                     ccACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcc                    <  1:603631/150‑1 (MQ=255)                                                                                                                                     ccACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCg                   <  2:503760/151‑1 (MQ=255)                                                                                                                                          gACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc              >  1:1079412/1‑151 (MQ=255)                                                                                                                                            ccGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg            <  2:82369/151‑1 (MQ=255)                                                                                                                                              gAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg          >  2:864730/1‑151 (MQ=255)                                                                                                                                                gccAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCCTAGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAACGCTATCCGGGCg        >  1:414555/3‑151 (MQ=255)                                                                                                                                                    ga‑gGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg    >  1:762256/1‑151 (MQ=255)                                                                                                                                                     a‑gGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg    <  1:192691/150‑1 (MQ=255)                                                                                                                                                       ggTTTTCCGCCACTGTCAGCCCTTCCGCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGc   <  1:301024/150‑1 (MQ=255)                                                                                                                                                       ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  >  1:103596/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                     CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGG‑GG‑CCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677890‑2678187

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