Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,103	G→C	intergenic (+71/‑531)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,103	0	G	C	95.3%	136.5 / ‑1.4	43	intergenic (+71/‑531)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base C (14/27);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCTCGACTTCTTACAGTGT‑CCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGAT  >  NC_000913/88861‑89342                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                gCTCGACTTCTTACAGTGT‑CCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGTCATCAATATCCTTAAAGt                                                                                                                                                                                                                                         <  1:340214/251‑1 (MQ=255) gCTCGACTTCTTACAG‑GTGCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTaaa                                                                                                                                                                                                                                           <  2:120843/249‑1 (MQ=255)  cgcGACTTCTTACAGGGT‑CCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTa                                                                                                                                                                                                                                        <  2:193678/249‑1 (MQ=255)      aCTTCTTACAGTGT‑CCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGAGTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGGAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATa                                                                                                                                                                                                                                    <  2:19027/251‑1 (MQ=255)            tactgtgt‑CCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGTACATCAATATCCTTAAAGTAAATAAGa                                                                                                                                                                                                                                 <  2:133245/244‑1 (MQ=255)                     cgggTGCTGGCAGTGGCTCAACGTCACAGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCGGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAAATAGACTATTCCTGAc                                                                                                                                                                                                                      <  2:313047/249‑1 (MQ=255)                         ggCTGGCAGTGGTTCAACTTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCGGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACtatt                                                                                                                                                                                                                  <  2:96791/250‑1 (MQ=255)                          gCTGGCAGTGGCTCAACGTCTCCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTTGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACtatt                                                                                                                                                                                                                  <  1:315085/250‑1 (MQ=255)                                 gTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaa                                                                                                                                                                                                          <  1:261498/251‑1 (MQ=255)                                               ccGCGATGTCGCAGAGCGTGTTCTGGAGATTGTCCTGGCAAACCTGGACGAACCGCGTAAGCCAAAACCTGTTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCTAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCTCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAccc                                                                                                                                                                                            <  1:326744/251‑1 (MQ=255)                                                cgcgATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAcccg                                                                                                                                                                                           <  1:201874/251‑1 (MQ=255)                                                      gTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGt                                                                                                                                                                                      <  1:234524/250‑1 (MQ=255)                                                           agggcttgtgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCTTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAAGACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGTACATAAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGttaatta                                                                                                                                                                                <  2:213883/245‑1 (MQ=255)                                                                   tgCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAcc                                                                                                                                                                         >  1:246101/1‑250 (MQ=255)                                                                                 cTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATACGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTcaca                                                                                                                                                           >  2:2098/1‑250 (MQ=255)                                                                                  tGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaat                                                                                                                                                         >  1:85940/1‑251 (MQ=255)                                                                                     caaaccTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTAAAACCCGCACGTTAATTAACACACCAGCTGAAATTTACaataat                                                                                                                                                      >  2:194169/5‑251 (MQ=255)                                                                                     cAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaatact                                                                                                                                                      >  1:283587/1‑249 (MQ=255)                                                                                                    cgcgTAAGCCAAAACCTGGTTTAACGCGAATTAAACGTAATCTCTACCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACa                                                                                                                                       <  2:196249/251‑1 (MQ=255)                                                                                                        tAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAgcgc                                                                                                                                   >  1:290522/1‑251 (MQ=255)                                                                                                                               gcATTAAACGAAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAATAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGCCTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCCGCTGAAATTCACAATAATTAAGTGATATCGACAGCCAGTTTTTGCATTATTTGGTTACAt                                                                                                            >  1:59063/1‑251 (MQ=255)                                                                                                                                 aTTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATCTCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTACGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgc                                                                                                          >  1:420946/1‑251 (MQ=255)                                                                                                                                    aaaCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAACCGCGAACAAAACTACGAGCCAGGTGCATTTCCCTCTGGCGCGTAGAGGACGTGAATGGACATCAATATCCTTAAAGAAAATAAACCTATTCCTGACTATTATTGATATGTTCTTTTAAACCAGCCCGGGAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATAGCGCCCTCGCGTTTTTGCATTATTTTGTAACagtc                                                                                                          >  2:35983/1‑245 (MQ=255)                                                                                                                                          aaTCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGACGCGAAAAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGa                                                                                                  <  2:430056/250‑1 (MQ=255)                                                                                                                                            ccgctATCGCCGCGGCGTGTTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAATTCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATGCAGAGATGAAtt                                                                                               <  2:238568/248‑1 (MQ=255)                                                                                                                                                 aTCGCCGCGGCGTGCTTAGCCGTATCTAAGCCGTGAACAAAAATACGCGTCAGGTTAATGTTCCTCTGGCGCGTAGAGTACGTGAATGGACTGCAATATCCGTAAATTAAATAAGACTATTCCTGACTATTATTGAGAAGTTCTATTAAACCCGCCCGTTAAGTAACTCATCAGCTGAAATTCACAATAATTAATGGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGcc                                                                                            <  2:336737/249‑1 (MQ=255)                                                                                                                                                    gccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGCACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATGCAGAGATGAATTGCCGAttt                                                                                       >  1:282891/1‑251 (MQ=255)                                                                                                                                                                 cAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACttt                                                                           <  1:309885/250‑1 (MQ=255)                                                                                                                                                                      gTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTCTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTg                                                                     <  1:363889/251‑1 (MQ=255)                                                                                                                                                                        aGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGTCGATTTAACAAACACTTTTCTTTGCt                                                                   <  2:3501/251‑1 (MQ=255)                                                                                                                                                                         gCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCtt                                                                  >  2:159098/1‑251 (MQ=255)                                                                                                                                                                          cTAAGCCGCGCACACAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGAACTTTTCTTTGCttt                                                                 >  2:36265/1‑251 (MQ=255)                                                                                                                                                                                     aCAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAccc                                                      <  2:214742/251‑1 (MQ=255)                                                                                                                                                                                          aaTACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTgg                                                 <  2:290522/251‑1 (MQ=255)                                                                                                                                                                                             aCGCGCCATGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTGAAAGTAAATAATACTATTCCTGACGATTATTGATAAGTTCGTTGAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATGAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCATAGATGAAGTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCAt                                              <  2:282891/251‑1 (MQ=255)                                                                                                                                                                                                            ttCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAg                                <  2:150198/250‑1 (MQ=255)                                                                                                                                                                                                                      ggcgTAGAGTACGTGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACaggg                     <  2:354659/250‑4 (MQ=255)                                                                                                                                                                                                                            gagTACGGGAATGGAAATCAATATCCTCAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGtt               <  2:326671/251‑1 (MQ=255)                                                                                                                                                                                                                            gagTACAGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAAGTAACTCACCAGCTGAAATTCTCAATAATTAAGTGATAGCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAAGTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGtt               <  2:393297/251‑1 (MQ=255)                                                                                                                                                                                                                                 cGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAACAGCTGAAATTCCCAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGCGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCg          >  2:402020/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGa   <  2:46942/251‑1 (MQ=255)                                                                                                                                                                                                                                         gACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGATTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGat  <  1:205287/251‑1 (MQ=255)                                                                                                                                                                                                                                          aCATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCAGCCCGTTATTTAACTCACCAGCTGAACTTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAAAAAAAACTTTTCTTTGCTTTTACGCAAACCCGCTGGCATAACACGCCACACAGACGCAACAGGGACTGTTACACGGGGAAGat  >  2:187326/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                GCTCGACTTCTTACAGTGT‑CCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTGGTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGAT  >  NC_000913/88861‑89342

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