Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,869,637	C→A	L87L (CTC→CTA)	cbrA →	colicin M resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,869,637	0	C	A	90.5%	87.8 / 2.9	42	L87L (CTC→CTA)	cbrA	colicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (17/21);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

ATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCAT  >  NC_000913/3869398‑3869883                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       attattGGCCTCGGCTCGGCAGGATCGGAGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCAACAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGAAGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGGCGCCGCATCGCTAAAGCGTAACt                                                                                                                                                                                                                                              >  1:339359/1‑250 (MQ=255)                  ggCAGGGTCGGCGTTGGCACGAAAGTTATCAGGCAAAATGCAGGTGATCGCTCTGGATAAAAAGCACCAGTGTTGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCGGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTCATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGGCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCa                                                                                                                                                                                                                            <  2:269431/251‑1 (MQ=255)                   gCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCCCCGCATCGCTAACGCGTAACCACCAGCGAAGCTATATCa                                                                                                                                                                                                                            >  2:75688/1‑250 (MQ=255)                    cAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAAGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAat                                                                                                                                                                                                                          <  1:348989/251‑1 (MQ=255)                       gATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATAtt                                                                                                                                                                                                                       >  1:151936/1‑251 (MQ=255)                                tGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCa                                                                                                                                                                                                               <  1:53099/250‑1 (MQ=255)                                      gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTTGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTc                                                                                                                                                                                                        <  1:384706/251‑1 (MQ=255)                                       aaaGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTCCTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCAACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCg                                                                                                                                                                                                       >  1:358767/1‑251 (MQ=255)                                          gTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACt                                                                                                                                                                                                    <  1:277188/251‑1 (MQ=255)                                                   cAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGa                                                                                                                                                                                           <  1:347549/251‑1 (MQ=255)                                                     aaaTGCAGGTGAACGCTCTGGATCAAAGGGTCCAGTGTGGGACTTAAGGTTTTAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCTCTGCGATCTTTTATTCGCGATGGGCTGACTCTGCCTGTCGATGTGTTCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCAGCGCTAACGCGTAACTACCAGCGAAGCTATATCAAGATTAATCGCCATGCTTTCGACTTGTGGATGaaa                                                                                                                                                                                         <  2:87494/251‑1 (MQ=255)                                                       aTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTGTATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATc                                                                                                                                                                                       <  1:209279/251‑1 (MQ=255)                                                                                                tGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGttt                                                                                                                                                             >  1:299451/1‑236 (MQ=255)                                                                                                ggAAGGTTTCACCACACCTTGCGGCGGTCGGCTGGCTCCGGACGCGCAGCGATCTTTTATGCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCGGATCCCCGCCAGCGTTGAGGttt                                                                                                                                                             <  2:299451/235‑1 (MQ=255)                                                                                                      tttCAGCAACCCCTGCGGCGGTCTGCTGGCAACGGACGAGCAACGATCTTTTATTCGCGATGGGCTGACGATTCCTGACGATGTGATCGCTAATCCGCAGATTTTCAGCGTCCAAACCGTCGACGTCGCCGCATCGCTAACTCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCCGCGGTGAGGTTATTCACAACCGCCTGTGCCGa                                                                                                                                         >  2:248109/1‑250 (MQ=255)                                                                                                        tcagcacccccTGCGGCGGCCTGCTGGCACCGGACGAGGAGCGATCTTTTATTCGCGATGGGCTGACTCTTCCGGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAACCCGTCGACGCCGCCGCATCGCTAACGCGTAACCCCCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGaaa                                                                                                                                       <  2:27236/242‑1 (MQ=255)                                                                                                                                gggCACCGGACGCGCAGCGATTTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACTTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                              <  1:22887/250‑1 (MQ=255)                                                                                                                                 ggCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTTACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAACCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATg                                                                                                               <  1:329120/249‑1 (MQ=255)                                                                                                                                 agcccctgactccctgctctCTTTTATGCTCGTTGGGCTGACGCTTCCTTTCTATTTAATCGCTAATCCCAAAATGTTCAGCGCCAAACCCGTCGCTGCCGCCGCATCGCTAACGCGTACCTCCCAGCGAAGATATATCAATATAAATCGCCCTGCTTTCGACTTGTGGATGAACTCGCTGATCCCCGCTAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATg                                                                                                               <  2:14538/230‑1 (MQ=255)                                                                                                                                       ggACTCGCAACGATCTATTATTCGCGATGGGCTGACGCTTCCTGTCGATGAGATAGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCACCGCTAACGCGTAACTACCAGCGAAGCTATAACAATATTAATCGCCAAGCTTTCGACTTGTGGATGAAATTGATGATCCTCGACATCGTTTAGGTTTATAACGCCAGCCTGTGCCGAACAATCTGGCGTGAGTATGATAAATGGCATgacc                                                                                                       >  2:50512/1‑248 (MQ=255)                                                                                                                                            cgcAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAttttt                                                                                                  >  1:447391/1‑251 (MQ=255)                                                                                                                                              cAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCg                                                                                                >  1:88798/1‑251 (MQ=255)                                                                                                                                              cAGCGATCTTTTATTCGCGATGGGCTGACGCTGCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAA‑CCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGt                                                                                               <  1:190559/251‑1 (MQ=255)                                                                                                                                                gCGATTTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTg                                                                                              <  1:75688/251‑1 (MQ=255)                                                                                                                                                                    ggCTGACGCATCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGCGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACa                                                                           >  2:87631/1‑250 (MQ=255)                                                                                                                                                                               ccTGTCGATGTGATCGCTAATCCGCAGTTTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCACTATTAATCGCCATGCTTTCGACTTGAGGATGAAATCGCTGATCCCCACAAGCGTTGCGGTTTATCCCAACAGCCTGTGCCGAAAAATCTGGCGGGAGGGTGATAAATGGCATGTCATTTTTCGTGCAGCCGGCTGGAGGCAACAAATTAccgcccc                                                               >  2:330560/1‑250 (MQ=255)                                                                                                                                                                               ccTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgcccg                                                               >  1:302351/1‑251 (MQ=255)                                                                                                                                                                                               cgcatccctagtgtttcggccacAAAACCGTCGGCGCCGCAGCATCCCTAACGCTGAACTACCAGCTAAGCTACATCAATATTAAGCGCCTTGTTTTCGTGTTGGGGATAAAAGCGCGTACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGGGAGGATGATAAAGGGCATGTCATTTTTCGTGCAGTCGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGc                                                <  2:444046/228‑1 (MQ=255)                                                                                                                                                                                                    ccGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTACCGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATgg                                          <  1:75846/251‑1 (MQ=255)                                                                                                                                                                                                      gCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGtt                                                                                                                                                                    >  1:423197/1‑127 (MQ=255)                                                                                                                                                                                                      gCAGAGTTTCAGCGTCAAAACCGCCGACGTCGCCGCATCCCTAACGCGTAACTTCCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTGGTGGATGAAATCGCTGACCCCCGCCTGCGtt                                                                                                                                                                    <  2:423197/127‑1 (MQ=255)                                                                                                                                                                                                                  cGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGCTgg                            >  1:151145/1‑251 (MQ=255)                                                                                                                                                                                                                     aaaaaCCGTCGACGCCGCCGCATCGCTACCGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCAATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTgc                         <  1:425399/250‑1 (MQ=255)                                                                                                                                                                                                                      aaaaCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTgcg                        >  1:308786/1‑251 (MQ=255)                                                                                                                                                                                                                          ccGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTgcggcgg                    >  1:68979/1‑250 (MQ=255)                                                                                                                                                                                                                           cGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGCTGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATg                             >  2:45731/1‑241 (MQ=255)                                                                                                                                                                                                                           cGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATg                             <  1:45731/241‑1 (MQ=255)                                                                                                                                                                                                                           cGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGAc                   >  2:154451/1‑251 (MQ=255)                                                                                                                                                                                                                             tCGACGGCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATACTTTCGACTTGTGGATGAAATCGCTGATCCCCGACAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAACATCTGGCGTGAGGATGCTAAATGGCATGTAATTTTTCGAGCAGACGGCTGGGAGCAACATATTACCGCGCGCTATCTGGTCGGGGCCGATGGCTCAAACTCGATGTTgcggcgcaat                 >  2:142319/1‑247 (MQ=255)                                                                                                                                                                                                                                    cgccgcATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTAc           <  1:248109/250‑1 (MQ=255)                                                                                                                                                                                                                                            cGCTAACGCGTACCTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGatca   <  1:288260/250‑1 (MQ=255)                                                                                                                                                                                                                                            cGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAACCGCCATGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGatcat  <  1:154451/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       ATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCAT  >  NC_000913/3869398‑3869883

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