Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,028,546	G→A	A156V (GCC→GTC)	ygfS ←	putative electron transport protein YgfS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,028,546	0	G	A	95.7%	144.9 / ‑5.5	47	A156V (GCC→GTC)	ygfS	putative electron transport protein YgfS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/18);  minor base C (2/0);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCA  >  NC_000913/3028300‑3028792                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        gCTATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGAcaac                                                                                                                                                                                                                                                     >  2:441188/1‑251 (MQ=255)   tATGGCTGGACTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGCATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGacaaca                                                                                                                                                                                                                                                    >  2:39419/1‑250 (MQ=255)             cTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGctgctg                                                                                                                                                                                                                                         >  2:11036/1‑251 (MQ=255)                              cATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGtt                                                                                                                                                                                                                        >  1:379285/1‑251 (MQ=255)                              cATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCCTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGtt                                                                                                                                                                                                                        <  2:9384/251‑1 (MQ=255)                              cATCGTTCGGACTTGGTATATCGGTACCGTTATTATGCTGGGCTTTGGGAAACTGGTTGTAGGCCGCGTCTTTGCACGCGAATGGCAACAGCTTAATATTGGTACGTTGAGCATTACTGCCGCCCTGGTCGCATTTTTCGCGGGGGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGTAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAAAAATACGCTACTCGCGGACTCGCCTGatt                                                                                                                                                                                                                         <  2:316159/250‑3 (MQ=255)                                      tGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTc                                                                                                                                                                                                                 >  1:468935/1‑250 (MQ=255)                                        aCTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGt                                                                                                                                                                                                              >  1:35464/1‑251 (MQ=255)                                              aaTATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCaa                                                                                                                                                                                                        <  2:35464/251‑1 (MQ=255)                                              aaTATCGTTACCGGTATTATGCGTGGTGTTGTGACACTGGTCGTAGGCCGCGTCTGTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGTGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTGCAGTAAGCaa                                                                                                                                                                                                        <  2:50313/251‑1 (MQ=255)                                                   cGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCa                                                                                                                                                                                                   <  2:1704/251‑1 (MQ=255)                                                    gTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGACTATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGc                                                                                                                                                                                                    >  1:544069/1‑249 (MQ=255)                                                                     gggCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCa                                                                                                                                                                                   >  2:230963/1‑249 (MQ=255)                                                                                           ggCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCctc                                                                                                                                                          >  2:159245/1‑251 (MQ=255)                                                                                                 gTCTTTGCAAGCGAATGGCAAAAGCTGAATATTGGTACGGCGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGAGTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGCCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGTTTCAACGCAG‑GCTGGCCCCTCTTCCCg                                                                                                                                                    >  2:187326/1‑251 (MQ=255)                                                                                                  tCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGc                                                                                                                                                   >  1:87414/1‑251 (MQ=255)                                                                                                  tCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACAGCAGCGCCTGTTTTGGGCAAGATTCAACGCCG‑GCTGGCCCCTCTTCCCGc                                                                                                                                                   >  2:417414/1‑251 (MQ=255)                                                                                                         ccgcgAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTcac                                                                                                                                            <  1:365623/250‑1 (MQ=255)                                                                                                              catagCAAAAGCGGAATATAGGTAGGATGATAATTACTGCAGCAATGGCCGCATTTTACGCGGGTGGTTGGGCAATCTAATGGATTCTCCGCGCGTTTAAGCGAGGAGATTTCTTTCAGAGAGGATTCTCACCGGTGACAACAATACGCTGCTGGCGGAATCGCATGAGTACTCTTTCAGTAACAAACTGCAGCGCCTGTGTTGGGGGAGCTTCACCG‑AGAGCTGGCCCCTCTTGCCGCTGTTCACag                                                                                                                                          >  2:62812/5‑248 (MQ=255)                                                                                                                  gCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGCCGCATTTTACGCGGGTGGTTAGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAAAAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑CCTGGCCCCTCTTCCCGCTGTTCACAGCGATCGc                                                                                                                                   >  1:274951/1‑251 (MQ=255)                                                                                                                         cTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGAAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTcac                                                                                                                            >  1:226590/1‑251 (MQ=255)                                                                                                                           gAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACaa                                                                                                                          <  1:441188/251‑1 (MQ=255)                                                                                                                           gAATATTGGTAAGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATCGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGACCCTATTCACGCTGTTCACAGAGATCgactttcaaa                                                                                                                           >  1:435170/1‑241 (MQ=255)                                                                                                                                 tGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACg                                                                                                                                                                                                                                               <  1:436108/129‑1 (MQ=255)                                                                                                                                 tGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACg                                                                                                                                                                                                                                               >  2:436108/1‑129 (MQ=255)                                                                                                                                     aCGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTc                                                                                                                 <  1:230963/250‑1 (MQ=255)                                                                                                                                       ggTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATGTCTTTCAGAGAGGATTCTCACCGCTTACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGTCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCggg                                                                                                              <  1:378282/251‑1 (MQ=255)                                                                                                                                                                   ttttACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGAt                                                                                  >  1:557386/1‑251 (MQ=255)                                                                                                                                                                     tttcgcgGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGCTGATTCTCACCGCTGACAACAATACGCTGTTGGCGGACTCGTCTGTGTTCTCTTTCAGTAAGCAACTGCAGTGCCTGTGTTTGTCAAGATTCAACTCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCTCATTTCACAATTTGTTGTCGGGTATCCCCCGGCATCGACTGAATGGTGATc                                                                                 <  1:86329/247‑1 (MQ=255)                                                                                                                                                                     ttACGCGGGAGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGCTTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGAAG‑GCTGGCCCCTATTCCCGCTGGTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCCACGACTGAATGGAGATc                                                                                 >  2:589227/1‑250 (MQ=255)                                                                                                                                                                        cgcgGGAGGCTGGACACTCTAATAGATTCTCCGCGCCTTTCAGCGAGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACACCAACACGCTGCTGGCCGACTCGCCTGCGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGCCTGAATGGTGATCCTTc                                                                             >  2:102657/1‑251 (MQ=255)                                                                                                                                                                         gcgGGTGGTTGGGCAATCTAAGAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTcc                                                                            <  2:11291/251‑1 (MQ=255)                                                                                                                                                                             gtggtTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAAc                                                                        <  1:176273/251‑1 (MQ=255)                                                                                                                                                                                      cAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGgcacgc                                                                >  1:177506/1‑250 (MQ=255)                                                                                                                                                                                            aaTAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGAc                                                          <  2:37499/250‑1 (MQ=255)                                                                                                                                                                                                 aTTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAAct                                                    >  1:412530/1‑251 (MQ=255)                                                                                                                                                                                                  ttCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCAACGCTGACAACAATACGCTGCTGGCGGACTCGACTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCCTTTCACAATTTGTTGTCGGGTATCCCCCGGCCACGACTGAATGGTGATCATTCCATACGGGCACGCACTGACGCAActc                                                   >  1:31915/1‑251 (MQ=255)                                                                                                                                                                                                            cTTTCAGTGCGAACATTTCTTTCAGCGAGGATTCTCCCCGCTGCCAACAATCCGCTGCTGGCGGCCTCCACTGTGATCTCTTTGAGTACGCAACTGCAGCGCCTGTGTTGGGCCACATTACACGCAG‑GCTGCCCCATGTTCACGCTGTTCACAGAAATCGCATTTCACACTTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGCTCTATCTAAACGGTCACGCACTgctgacactatgacgaccac                                         >  2:31234/1‑232 (MQ=255)                                                                                                                                                                                                              ttCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAAt                                        <  1:11036/250‑1 (MQ=255)                                                                                                                                                                                                                cAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACACTACGCTGCTGGCGGAATCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCACGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGGTCATTCCAACCGGGCACGCACGGACGCAACTATGAACGCCAATAc                                      >  1:21936/1‑250 (MQ=255)                                                                                                                                                                                                                       gagaTTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGAACGCACTGACGCAACTCTGACAGCAAATACAACGGGc                               >  2:21805/1‑250 (MQ=255)                                                                                                                                                                                                                       gagaTTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGGGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGc                               >  2:476630/1‑250 (MQ=255)                                                                                                                                                                                                                              ctttcAGAGAGGATTCTCACCGCTGACAACACTACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTc                                                                                                                 >  1:495198/1‑161 (MQ=255)                                                                                                                                                                                                                              ctttcAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCTGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTc                                                                                                                 <  2:495198/161‑1 (MQ=255)                                                                                                                                                                                                                                  cAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTg                   <  1:278553/251‑1 (MQ=255)                                                                                                                                                                                                                                       agGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGACTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGcacc               >  1:537763/1‑250 (MQ=255)                                                                                                                                                                                                                                                aCCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGc     >  1:8183/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCa  <  1:516451/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GCTATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAG‑GCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCTCGCCCA  >  NC_000913/3028300‑3028792

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