Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,552,453	G→A	A77V (GCA→GTA)	yddM ←	DNA‑binding transcriptional regulator YddM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,552,453	0	G	A	98.1%	185.0 / ‑4.9	53	A77V (GCA→GTA)	yddM	DNA‑binding transcriptional regulator YddM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/22);  minor base T (0/1);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTCATTAGTGTCTTTA  >  NC_000913/1552211‑1552693                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 gTCATTTTCATATGGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTGTTCTGCCTCTACCAGACTc                                                                                                                                                                                                                                          <  2:302173/251‑1 (MQ=255)    aTTTTCATATTGTAGACAACGTAGGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCTCTTTATAAAGTTCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCCCTGTTTTTTCTGCCTCTACCAGACTCCAg                                                                                                                                                                                                                                       <  1:520753/251‑1 (MQ=255)       ttCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGc                                                                                                                                                                                                                                     <  2:67261/250‑1 (MQ=255)          atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGttt                                                                                                                                                                                                                                 >  1:450173/1‑251 (MQ=255)            atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTGTCTGCCTCTACCAGACGCCAGGCGtttt                                                                                                                                                                                                                                <  2:308608/250‑1 (MQ=255)                           ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACAt                                                                                                                                                                                                                >  1:409751/1‑251 (MQ=255)                                 gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcgg                                                                                                                                                                                                          <  1:514533/251‑1 (MQ=255)                                                            ttATCCGGCATGAAAACCCTTCAAATACCATAGATTGCCGTGCACGTGTAGGCCCGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCCTCAGTTCTAAATGGCGCTTTCTAAAGTGACATTTTTTTATTGCGTAACCAGACGGCGTACTCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCCCATTTGCGGCGAACTGCCGATCACCACGGAAAGttt                                                                                                                                                                               >  2:234571/1‑251 (MQ=255)                                                                 cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCg                                                                                                                                                                          >  2:513828/1‑251 (MQ=255)                                                                  ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGc                                                                                                                                                                         >  2:386554/1‑251 (MQ=255)                                                                        aaaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTc                                                                                                                                                                   <  2:147565/251‑1 (MQ=255)                                                                         aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCt                                                                                                                                                                  <  1:463762/251‑1 (MQ=255)                                                                           aCCCTTCAAATCCAATCGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTgg                                                                                                                                                                >  1:434791/1‑251 (MQ=255)                                                                                  aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCCTAAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCaa                                                                                                                                                          >  1:590946/1‑250 (MQ=255)                                                                                          tAGATTGCAGAGACCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCATCCGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGACCTGCCGATCACCACGGCAAGTTTAATCGCCTTTTCTGGCGTCAAAGCTGCttt                                                                                                                                                 >  2:192451/1‑251 (MQ=39)                                                                                                   gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGtt                                                                                                                                                                                >  1:228393/1‑211 (MQ=255)                                                                                                   gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGtt                                                                                                                                                                                <  2:228393/211‑1 (MQ=255)                                                                                                    tGAACGTGTAGGCCTGATCAGCGTAGCGCATCAGGCCATTATGCGTTTGCCATCAGTTCTAAATGGCCCTTTATAAAGTGCCATTTTTTTATTGCGTCACCAGACGGCGTAATCGCGCCCCATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCTACTTTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGc                                                                                                                                       >  1:372074/1‑251 (MQ=255)                                                                                                      aaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCa                                                                                                                                      >  2:128795/1‑250 (MQ=255)                                                                                                            tAGGCATGATAAGCGTCGCGAATCAGGCACTGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCAATTTTTTTATTGCGTAACCCGACGGCGTAATCGCGACACATACACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTCTGCAGATTCAGCAACATTTGCGGCGAACTGCCGATCACCACGGCAAGCTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCTTGTCAGCatgcggct                                                                                                                               >  2:307793/1‑244 (MQ=255)                                                                                                             aGGCCTGATAAGCGTAGAGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCtt                                                                                                                              >  1:508475/1‑251 (MQ=255)                                                                                                              ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTg                                                                                                                             <  1:513828/251‑1 (MQ=255)                                                                                                                ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGcc                                                                                                                           <  1:345032/251‑1 (MQ=255)                                                                                                                   gATAAGCGTAGCGCATCAGGCAAGGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTGTCTTCCTCTACCAGACTCCAGGCGTGTTGCAGATTCAGCCACATTTGCGGCTAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGTCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGtt                                                                                                                        <  2:375045/251‑1 (MQ=255)                                                                                                                            aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTCAATGGCGCTTTATAAAGTGACATTTTTTTCTTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGc                                                                                                                >  2:559985/1‑250 (MQ=255)                                                                                                                             gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCaa                                                                                                              <  2:110245/251‑1 (MQ=255)                                                                                                                                           gTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTg                                                                                                <  1:156671/251‑1 (MQ=255)                                                                                                                                             tGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACGGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTgg                                                                                               <  2:198694/250‑1 (MQ=255)                                                                                                                                              gCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGc                                                                                              >  2:65317/1‑250 (MQ=255)                                                                                                                                                                           tataAAGTGCCATTTTTTTATTGCGTATCCAGACGGCATAATCGCGACACATCCACTGTTTTTTCTTCCTCTACCAGTCTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACGGCCGATCACCACGGATAGTTTAGGCTCCGTTTCTGGCGGCTATTCTGCTTTTCCTGTCAGCAAGCGTCTTGCCGTTGATGGATCAAGTTCCGTTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGt                                                                 <  1:79048/250‑1 (MQ=255)                                                                                                                                                                            ataAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTcc                                                               >  1:508617/1‑251 (MQ=255)                                                                                                                                                                             taAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTcc                                                               >  1:528132/1‑250 (MQ=255)                                                                                                                                                                                aGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGAGTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAAGCGCCATTTCGGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTGGCCGTTGAGGGAGCAATGTCCATTGCTCTGGCAAACTCGCGCAGGCGGACATTGAGTTCGTCCAGTg                                                           <  2:590946/251‑1 (MQ=255)                                                                                                                                                                                  tGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGAt                                                         >  2:272992/1‑251 (MQ=255)                                                                                                                                                                                  tGCCATTTTTTTATTGAGTAACCAGACGGCGTATTCGCGACCCATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGCTTCAGCCACATTTGCGGAGAACTGCCGCTCACCACGGAAAGATTAATTGCCATTTCTGGCGTCAATGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGATGGAGCAATTTCGATTGCTCGGGGACACTCGCGCTGGCTGCCATAAAGTTCGTCCAGTGAt                                                         >  2:281785/1‑251 (MQ=255)                                                                                                                                                                                      aTTTTTTTATTGCGTACCCAGACGGCGTAAGAGCGACACACCCACTGTTTTTGCTGCCTCTACCAGACTCCAGGCTTTTTGCAGATTCAGCCACATGTGCGGCGAACTGCCGATCACCACGGAAAGTTTAAGCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCTCAGGCTGACATTAAGTTCGTCCAGTGATTCCt                                                     <  2:13427/251‑1 (MQ=255)                                                                                                                                                                                              aTTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCTTTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATg                                             >  2:264844/1‑251 (MQ=255)                                                                                                                                                                                                tGCGTAACCTGACGACGTCATCGCGAAACATCCAATGTTTTTTCTGACTCTACGCGCCTCCAGGAGCTTTGCCGATTCAGCCACCTTTGCCGCGAACTGCCGATCACCACGGCAAGTTTAATCGCCCTTTCTGGCGTCCAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATACTGTc                                           >  2:565471/1‑251 (MQ=255)                                                                                                                                                                                                tGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGAGTCCTGAATAATGt                                            <  1:456723/250‑1 (MQ=255)                                                                                                                                                                                                 gCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTc                                           <  2:448613/250‑1 (MQ=255)                                                                                                                                                                                                     aaCCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGACTCCACCAGAATCCAGGCGTTTTGCAGATTCCGCCACATTTGCGGCGCACTGCCGATCCCCACGGAAAGTTTCATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCcg                                      >  1:284847/1‑251 (MQ=255)                                                                                                                                                                                                           aCGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGCTTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGCAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCcgggcggg                                >  2:498491/1‑251 (MQ=255)                                                                                                                                                                                                              gCGTAATCGCGACACACCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCCCGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGgat                             <  1:164289/251‑1 (MQ=255)                                                                                                                                                                                                                gTAATCGCGACACCTCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCAAGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGGGATTCCTGAATAATGTCCCCCGGGCGGGgatgg                           >  2:90281/1‑250 (MQ=255)                                                                                                                                                                                                                gTAATCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGgatg                            >  1:421814/1‑250 (MQ=255)                                                                                                                                                                                                                  aaTCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGGGATTCCTGAATAATGTCCCCCGGGCGGGgatgat                          <  2:107828/250‑1 (MQ=255)                                                                                                                                                                                                                   aTCGCGACACATCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTg                        >  1:294606/1‑251 (MQ=255)                                                                                                                                                                                                                       cgacacctcccgatgtttttCTGCCTCTTCCAGCCACCAGGCGGTTTGGAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTGTCATCGCAATGTCTGTCGTCCAAGCTGCTTTTCCTTGCAGCAATCGGCTTGCCGTTGTGGGATCAATTTCCATTGCTCTGTCGAACTCGCGCATGCTGACATTAATGTCGTCCAGTGTTTCCTAAAGAATGACGCCCGTGCGGGGAGGATTGGCCa                    <  2:241994/236‑1 (MQ=255)                                                                                                                                                                                                                            caTCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTc               <  2:450173/251‑1 (MQ=255)                                                                                                                                                                                                                              tCCACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTcct             >  2:18307/1‑249 (MQ=255)                                                                                                                                                                                                                                cACTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACCTTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTCAtt            >  1:38489/1‑250 (MQ=255)                                                                                                                                                                                                                                 aCTGTTTTTTCTGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCCGCCAAATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTCAtt            >  1:250713/1‑249 (MQ=255)                                                                                                                                                                                                                                     ttttttATGCCTCTACCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTACTCGCCCTTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGACTAATGTCCCCCGGGCGGGGATGCTTGGCCATTTTCATTAgtgtt      >  2:480603/1‑250 (MQ=255)                                                                                                                                                                                                                                          tctgcctctaccAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTCATTAGTGTCTTTa  <  2:294606/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 GTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTCTGGCAAACTCGCGCAGGCTGACATTAAGTTCGTCCAGTGATTCCTGAATAATGTCCCCCGGGCGGGGATGATTGGCCATTTTCATTAGTGTCTTTA  >  NC_000913/1552211‑1552693

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