Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,064,246	C→T	E74E (GAG→GAA)	cobT ←	nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,064,246	0	C	T	97.9%	144.9 / ‑4.0	48	E74E (GAG→GAA)	cobT	nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/21);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGC  >  NC_000913/2064000‑2064485                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGTGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑ccc                                                                                                                                                                                                                                             <  1:166247/251‑1 (MQ=255) gACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑ccc                                                                                                                                                                                                                                             >  2:369331/1‑251 (MQ=255)      ccAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACg                                                                                                                                                                                                                                        <  1:522351/251‑1 (MQ=255)             cAACTTTTTAGCCTGACGTCGACTCATTGACGGAGCTGAAGCTATATTGCCGCTTCTTCGTGCGAAACGCATGGTGATAAGCTCGGGGATAGGCGAAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCCCACACTCGCCGGTTGTTCCACGTGTAATATTGTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGag                                                                                                                                                                                                                                 <  1:308253/251‑2 (MQ=255)                  tttCAGCCTGACGGCGACTCATTGCCGGAGCTGACGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGc                                                                                                                                                                                                                             >  2:65214/1‑250 (MQ=255)                  tttCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCgg                                                                                                                                                                                                                              <  1:400732/249‑1 (MQ=255)                     cAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGcacac                                                                                                                                                                                                                         <  1:165579/251‑1 (MQ=255)                       gtcTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCTTGCTACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCATCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTGGGAGAAATAGCGACCCCT‑CCTGCCCAGACGCCGTGATCGGcacacag                                                                                                                                                                                                                       <  2:152908/249‑2 (MQ=255)                              ggCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATATGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCCTGATCATCACACATAACCAGt                                                                                                                                                                                                                >  2:152760/1‑251 (MQ=255)                                ccaCTCATTGCCGGGGCTGAAACAACATTGCCTCTACCCCCTGCGCCACGCAAGTTGATAAGCCCGGGGATCGGCTCAGCAGTATCAATACCAACATCAATTCCGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGCCCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTAc                                                                                                                                                                                                              <  2:163957/249‑1 (MQ=255)                                                  gAAGCAATATTGCCGCTCCCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATTCCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACAt                                                                                                                                                                                             <  1:36827/250‑1 (MQ=255)                                                             gCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCATTCTT‑CCCAGACGCCGTGATCGGCACACATAACCCGTACCGCTTTTTTGCCCACATGCGGTATGCCAt                                                                                                                                                                                 >  2:478198/1‑251 (MQ=255)                                                               cGCTTCCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACTGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTGTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTAGGCCATTc                                                                                                                                                                               <  2:277040/251‑1 (MQ=255)                                                                             acacGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACCTGGACGTTAGCGCCCGCTTGGTCTGCCAGCCCACACACGCCGGTTGTTCCACTGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATTGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATATCCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCTGCATTcc                                                                                                                                                                 <  1:344254/251‑1 (MQ=255)                                                                               acGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTg                                                                                                                                                               >  1:104568/1‑251 (MQ=255)                                                                                          tAAGCTCGGGGATAGGCTAAGCATTATCAATACCAGCATGAATTACGTGGACGTTTGCGCCCGCTTGTTCTGCCAGCACACACTCGCCGGTTGTTCCACGGGTCATATTTTTAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGCCGCCGTGATCGGca                                                                                                                                                                                                                            <  1:581763/179‑1 (MQ=255)                                                                                          tAAGCCCGGGGACTGACTCTGCAGGATCACTACCAGCATCAATTCCGTGCACGGTAGCGCACGATTGTTCTGCCAGCACACACACGCAGGTCACTCCACTGGTACTATTTTCAGCATGTATGGCTGTCACTTCTTTTGGAGAAATATCGCCCCCTTCTT‑CCCAGACGCCGTGATCGGca                                                                                                                                                                                                                            >  2:581763/1‑179 (MQ=255)                                                                                               ccGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCa                                                                                                                                                 >  2:212349/1‑249 (MQ=255)                                                                                                        ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAg                                                                                                                                      <  1:65214/251‑1 (MQ=255)                                                                                                         gCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGc                                                                                                                                     >  1:62179/1‑251 (MQ=255)                                                                                                          cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACAAACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATCGCTGTCACTTCTTTTGGAGGAATAGCGACACCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCt                                                                                                                                    >  1:532677/1‑251 (MQ=255)                                                                                                              gcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTcc                                                                                                                                <  2:104568/251‑1 (MQ=255)                                                                                                                         aCCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCa                                                                                                                      <  1:498072/250‑1 (MQ=255)                                                                                                                                 cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGttt                                                                                                              <  1:369331/250‑1 (MQ=255)                                                                                                                                          ggACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTGTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCTCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATAGCATGCACCTCCAGCTTGCCCAGGCTACTAACAGGTTTGAGTAACCCg                                                                                                    <  2:486140/251‑1 (MQ=255)                                                                                                                                           gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTATT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCACACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGATTGCGTAACCCGt                                                                                                   >  2:169066/1‑251 (MQ=255)                                                                                                                                           gACGTTAGCGCCCGCTTGTACTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACCTGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCACCAGGTTTGAGTacaccgt                                                                                                   >  2:337277/1‑245 (MQ=255)                                                                                                                                              gTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCCCACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCa                                                                                                 >  1:214310/1‑250 (MQ=255)                                                                                                                                               ttAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGCATCCTATTTTCAGCCTGAATGGCTACCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CACAGACGCCGCGATCGGCACACATAACCCGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCACCCCCGGAATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCACCAGGTTTGAGAAACCCGTCAat                                                                                               >  1:49195/1‑251 (MQ=255)                                                                                                                                                     cccGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTc                                                                                                                                         >  1:100768/1‑203 (MQ=255)                                                                                                                                                     cccGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCg                                                                                         >  1:227355/1‑251 (MQ=255)                                                                                                                                                        gCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGCCATATTTTCAGCCTGTATAGCTGTCACTTCTTTTGGAGAAATAGCGCCCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGAATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTtg                                                                                      >  2:587210/1‑251 (MQ=255)                                                                                                                                                                tGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGCCGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCACAAGGTTTGAGTAACCCGTCAATATGCAGTTGTGCAcgcg                                                                              >  2:375970/1‑251 (MQ=255)                                                                                                                                                                            cacGCTGGTTATTCCACGGGTAATATTTTCAGCCTGTATGGCTGTCACTTGTTTTGGAGAAATAGCGACACATTCTT‑ACCAGACGCCGTGATCGGCAAACATCACCAGTACCGCTTTTTTGCCCACATGCGGTATGCAATTCAACCCCGGCATTCCTGCCAGATGTAAGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCCCGACACATAGCGGt                                                                    >  2:382411/1‑249 (MQ=255)                                                                                                                                                                                  ggTTGTTCAACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCg                                                             >  1:496203/1‑250 (MQ=255)                                                                                                                                                                                       ttCCAAGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCt                                                       >  2:37163/1‑251 (MQ=255)                                                                                                                                                                                         ccACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCAGTCCTGCCAGTTGTATGGCAAGCACCGCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTgg                                                     <  2:59761/251‑1 (MQ=255)                                                                                                                                                                                            cGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGAt                                                  <  2:426291/251‑1 (MQ=255)                                                                                                                                                                                                      tttCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAg                                         <  1:152760/250‑1 (MQ=255)                                                                                                                                                                                                       ttCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGt                                        >  1:172879/1‑250 (MQ=255)                                                                                                                                                                                                        tCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTa                                       <  1:73986/250‑1 (MQ=255)                                                                                                                                                                                                              tGTATGTCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGc                                >  2:73614/1‑251 (MQ=255)                                                                                                                                                                                                                       gTCACTTCTTTTGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGATTTCCAAGGCTACCAACAGGGTTGGGTAACCCGTCAATATGCCGTTGTGCACGCGACATCGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTg                       >  1:376981/1‑251 (MQ=255)                                                                                                                                                                                                                               ttttGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCAc               >  2:273497/1‑251 (MQ=255)                                                                                                                                                                                                                                 ttGGAGAAATAGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATGCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAAGCGGCAAGTATTTGCATCTCACg              <  1:85907/250‑1 (MQ=255)                                                                                                                                                                                                                                           aGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCCCCTCCATCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGTTATAATCGATGTCTGGGACCGTATTCAGTAAATCGGCAAGTAGTTGCATCTCTCGTCTCTTtata    <  2:376981/250‑1 (MQ=255)                                                                                                                                                                                                                                           aGCGACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTtat     >  1:536083/1‑249 (MQ=255)                                                                                                                                                                                                                                              gACCCCTTCTT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGc  >  1:508162/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 GACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCT‑CCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGTATGGCAAGCACCTCCAGCTTTCCCAGGCTACCAACAGGTTTGAGTAACCCGTCAATATGCCGTTGTGCACGCGACATAGCGGTAGAATCGATGGCTGGGATCGTATTCAGTAAATCGGCAAGTATTTGCATCTCACGTCTCTTTATAGC  >  NC_000913/2064000‑2064485

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