Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,100,070	G→A	C144Y (TGT→TAT)	ycdY →	chaperone protein YcdY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,100,070	0	G	A	97.9%	152.4 / NA	48	C144Y (TGT→TAT)	ycdY	chaperone protein YcdY
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (26/20);  minor base C (1/0);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATATGACGCAAGTCTCTGCCGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCT  >  NC_000913/1099803‑1100361                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    gtgATATGACGCAAGTCTCTGCCGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCAt                                                                                                                                                                                                                                                                      >  1:361267/1‑299 (MQ=255)                       cGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGTTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCGCTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCAGGCACCCACGCCTCTCTGGAGCAc                                                                                                                                                                                                                                               >  2:71964/1‑300 (MQ=255)                       cGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACAATGTTCAGTGAGTATCTGTTACCATGGTATGGTGTGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCAc                                                                                                                                                                                                                                               >  2:71987/1‑300 (MQ=255)                                 gCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCAcc                                                                                                                                                                                                                                      >  2:241080/1‑299 (MQ=255)                                 gCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCAAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTNTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCAc                                                                                                                                                                                                                                       >  2:487144/1‑298 (MQ=255)                                   gttgttTATCGGCGATGAATGTGCTGTGCCGCCATATTGTAGCGCATTGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCGTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCTGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCt                                                                                                                                                                                                                                   <  2:156582/300‑1 (MQ=255)                                      gttTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAAc                                                                                                                                                                                                                                >  1:289534/1‑300 (MQ=255)                                                        tgctgtgcCGCCATATCGTAGCGCATGGTTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCTAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                              <  2:512372/300‑1 (MQ=255)                                                        tgctgtgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGTAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATTAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGTTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCATGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                              <  2:312788/300‑1 (MQ=255)                                                        tgctgtgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGTAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                              <  2:313346/300‑1 (MQ=255)                                                        tgctgtgcCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGc                                                                                                                                                                                                              <  1:580421/300‑1 (MQ=255)                                                                 gccATATCGTAGCGCATAGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCAAATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTATGGCGCACCATGGCACCGCTAACCCGCGATGCCGTTAGTGCAATGTGGGa                                                                                                                                                                                                     >  2:330104/1‑300 (MQ=255)                                                                  ccATATCGTAGCGCACGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACACCTTTCTGGCGCCCCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGAc                                                                                                                                                                                                    >  2:259239/1‑300 (MQ=255)                                                                    atatCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAAACACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGa                                                                                                                                                                                                  >  2:527699/1‑300 (MQ=255)                                                                               cATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAAACCGCGATGCCATTAGTGCAATGTGGGACGAGCTGgaagaa                                                                                                                                                                                        >  2:90386/1‑299 (MQ=255)                                                                                                aCGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAgtg                                                                                                                                                                      <  2:364024/300‑1 (MQ=255)                                                                                                    aagcggaagTGCGCGCTTTTTTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCGGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTNGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGAt                                                                                                                                                                  <  2:185045/300‑1 (MQ=255)                                                                                                    aagcggaagTGCGCGCTTTTCTTTCTGAGCGATGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTNGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTTGCACACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTTAAGAATAAGTGTGAt                                                                                                                                                                  <  2:185514/300‑1 (MQ=255)                                                                                                                    ttttCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAAc                                                                                                                                                  <  1:422423/300‑1 (MQ=255)                                                                                                                                   aGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAttt                                                                                                                                   >  1:291266/1‑300 (MQ=255)                                                                                                                                          ccATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCAATCATAACTGCAACGGATAATttgttgtt                                                                                                                             >  1:228087/1‑299 (MQ=255)                                                                                                                                                  ggATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCAATCATAACTGCAACGGATAATTTGTTGTTGCATAAAAt                                                                                                                    >  1:264301/1‑300 (MQ=255)                                                                                                                                                   gATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtg                                                                                                                   <  2:616115/300‑1 (MQ=255)                                                                                                                                                      ctgcgagatgtgtataagagacagttGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAtgtgt                                                                                                                >  1:106103/25‑300 (MQ=255)                                                                                                                                                              cGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAGTGTGTGCTCGAtc                                                                                                        >  2:592175/1‑300 (MQ=255)                                                                                                                                                               gATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTATGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAAAGGATAATTTGTTGTTGCATAAAGTGTGTGCTCGAtct                                                                                                       >  2:529475/1‑300 (MQ=255)                                                                                                                                                                         gCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGc                                                                                             <  2:522489/300‑1 (MQ=255)                                                                                                                                                                                            ttCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATg                                                                          <  2:91128/300‑1 (MQ=255)                                                                                                                                                                                                                         aTGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGNCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc                                             <  1:439136/300‑1 (MQ=255)                                                                                                                                                                                                                         aTGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGNCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTGGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc                                             <  1:439074/300‑1 (MQ=255)                                                                                                                                                                                                                         aTGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTAGGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc                                             <  1:95914/300‑1 (MQ=255)                                                                                                                                                                                                                         aTGAGAGCGAAACACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGNCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGTACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc                                             <  1:439065/300‑1 (MQ=255)                                                                                                                                                                                                                           gagaGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCAGGCGCGGCATGAACATACTTCTCTCTATTGCaa                                           <  1:246477/300‑1 (MQ=255)                                                                                                                                                                                                                                       cTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAGTTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGTGTTTTCTGCTATCATGCGCGGCATGAAAATACCTCTCTATATTTCCCTCACAACAGGc                                >  2:2067/1‑299 (MQ=255)                                                                                                                                                                                                                                       cTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTATATTGCAATCACAACAGGCa                               >  2:2757/1‑300 (MQ=255)                                                                                                                                                                                                                                        tGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTTTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCa                               <  1:222438/299‑1 (MQ=255)                                                                                                                                                                                                                                           aaaCACTGTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATtct                           >  1:69076/1‑300 (MQ=255)                                                                                                                                                                                                                                                  gTTCAGTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGt                    >  1:460795/1‑300 (MQ=255)                                                                                                                                                                                                                                                       gTGAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCt                >  2:531586/1‑299 (MQ=255)                                                                                                                                                                                                                                                         gAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTTGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTggg             <  1:423156/300‑1 (MQ=255)                                                                                                                                                                                                                                                         gAGTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAAAATATTTATCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCt                >  1:428374/1‑297 (MQ=255)                                                                                                                                                                                                                                                          aGTATCTGTTACCCTGTTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTgggg            <  2:287540/300‑1 (MQ=255)                                                                                                                                                                                                                                                           gtgtttGTTACCCTGTTGTGGTGCGTTCCTTGGCAAATTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTNAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGAGAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTGTTCTGCTATCATGCGCGGCATGAACATACTTCGCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGa           <  2:200071/295‑1 (MQ=255)                                                                                                                                                                                                                                                           gTATCTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTACGGTATCTGGGGa           >  2:370126/1‑300 (MQ=255)                                                                                                                                                                                                                                                               cTGTTACCCTGGTCTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCAATCATAACTGCAANGGATAATTTGTTGTTGCATAAAGTGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTACGGTCTCTGGGGAt          >  2:196910/1‑297 (MQ=255)                                                                                                                                                                                                                                                               cTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAANGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCACTCTCTCCGGTATCTGGGGAt          >  2:197732/1‑297 (MQ=255)                                                                                                                                                                                                                                                               cTGTTACCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCACTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAANGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGAt          >  2:197125/1‑297 (MQ=255)                                                                                                                                                                                                                                                                    aCCCTGGTATGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCAATCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCt  >  2:518311/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    GTGATATGACGCAAGTCTCTGCCGATTACAATGCGTTGTTTATCGGCGATGAATGTGCTGTGCCGCCATATCGTAGCGCATGGGTTGAGGGCGCGACGGAAGCGGAAGTGCGCGCTTTTCTTTCTGAGCGAGGGATGCCATTAGCGGATACGCCAGCCGATCACATCGGCACATTGCTGCTCGCAGCTTCCTGGCTGGAAGATCAGTCAACGGAAGATGAGAGCGAAGCACTGGAAACACTGTTCAGTGAGTATCTGTTACCCTGGTGTGGTGCGTTCCTTGGCAAAGTGGAGGCCCATGCAACCACGCCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCT  >  NC_000913/1099803‑1100361

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