Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,973,874	C→T	A89A (GCC→GCT)	rffH →	glucose‑1‑phosphate thymidylyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,973,874	0	C	T	100.0%	108.2 / NA	36	A89A (GCC→GCT)	rffH	glucose‑1‑phosphate thymidylyltransferase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/20);  total (16/20)

CGGGCGGTTCCGGCACCCGATTGCATCCGATTACGCGCGGCGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTG  >  NC_000913/3973627‑3974117                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    cgggcggTTCCTGCACCCGCTTGCATCCGATTACGCGCGGCGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCAGGAAGATAAAGGCTATTTCCAGCGCCTGCTCGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCtttt                                                                                                                                                                                                                                                  >  1:190402/1‑251 (MQ=255)                    aTTGCATCCGATTACGCGCGGCGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAAccttcc                                                                                                                                                                                                                               <  1:73043/251‑1 (MQ=255)                      tGCATCCGATTACGCGCGGTGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTGACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCAGCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCTGTGAAAccttcct                                                                                                                                                                                                                              <  2:231708/250‑1 (MQ=255)                       gCATCCGATTACGCGCGGCGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAAccttcctt                                                                                                                                                                                                                             <  2:223526/250‑1 (MQ=255)                                    cgcgGCGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACcttc                                                                                                                                                                                                               >  2:147417/1‑251 (MQ=255)                                       ggcggatctaagctactttagCCGATTTACAATAAGCCAGTGACTTGCTATCCGCTGTCGGTGCTGATGAGTGCCGGTAGCCTCGAAATTCTCATCATCACTACGCCGGAAGAGAAAGGTTATTTCCAGCGCCTGCTGGGCGAGGTTAGTGAGTTCGGGATCCATCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGTTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTg                                                                                                                                                                                                            <  1:255964/231‑1 (MQ=255)                                                aaGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCACTGCCGGTGCTGATGCTGGCCGGAATCAGCGAAATTCTCATCCTCACTACGCCGGCAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACATTCTTGTCTGGTGCt                                                                                                                                                                                                   >  1:20373/1‑251 (MQ=255)                                                 aGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTg                                                                                                                                                                                                  <  1:162651/251‑1 (MQ=255)                                                 aGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGACGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCATTAATGGTGAACCTTCTTGTCTGGTGCTg                                                                                                                                                                                                  >  1:164828/1‑251 (MQ=255)                                                           gCCGATTTACGATAAGACAATGATATACTATCCGCTGTCGGTGCTGCTGCTGGCCGGTATCCGCGAAATTCTCATCATCTCTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGACGATGGTAGTGAGTTCGGTATACAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGTAGCTTCTTGTCTGGTGCGGGGCGATAAc                                                                                                                                                                                         >  2:190076/1‑250 (MQ=255)                                                             cGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGAGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGCCGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtc                                                                                                                                                                                      >  1:88462/1‑251 (MQ=255)                                                             cGATTTACGATAAGCCAATGATTTAATATCAGCAGTCGGTGCAGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTCTTTCCAGCGCCTGCTGGGCGATGGGAGTGAGTTCGGTAACCAGCTGGAATATGCCGAACAACCCAGCCCGGACGGTCTGGCGCAGACTTTTATCATCGGTGAAAACTTCCTTAATGGTGACACTTCTTCTCTGGGGCTGGGCGATAACAtc                                                                                                                                                                                      >  2:233020/1‑251 (MQ=255)                                                                            cAATGATTTACTCTCCTCTTGCGGTGCTGATGCTGGCCGGGATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTGCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCTAGCTGGAATATGCCGAACAGCCCAGCTCGGACGTGCTTGCGCAGGCTTTTATCATCGGTGAAACCGTCCTTAATGGTGAACCTGCTTGTCTGTTGCTGGTCGATAACATCTTCGGCGGGCAGGGc                                                                                                                                                                       <  2:238456/251‑1 (MQ=255)                                                                                                    tGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGtt                                                                                                                                               >  1:136904/1‑251 (MQ=255)                                                                                                             tGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGgcgcgc                                                                                                                                      >  2:94732/1‑251 (MQ=255)                                                                                                               gCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCAc                                                                                                                                    >  2:146418/1‑251 (MQ=255)                                                                                                                   gTATCCGCGAAATTCTCATCATCACTACGACGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGCTGGTAGCACGTTCGGTATCCAGATGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCACAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGAGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCGGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGACGCGCAcggc                                                                                                                                 >  1:113524/1‑247 (MQ=255)                                                                                                                            aaaTTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACg                                                                                                                       >  2:69360/1‑251 (MQ=255)                                                                                                                              aTTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTGCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGTTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTTCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGt                                                                                                                     <  1:262442/251‑1 (MQ=255)                                                                                                                                                  ggAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATgg                                                                                                 <  2:116826/251‑1 (MQ=255)                                                                                                                                                         aaaGGGTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTGCGGTAGCCAGCTGGAATATGCCGAACAGCCCAGCCCGTACGGTCGGGCGCAGGCTTTTATCATCGGTGAATCCTTCCGTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGa                                                                                          <  2:221859/251‑1 (MQ=255)                                                                                                                                                                   tCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAACCCTTCCTTAATGGTGAACCTTCTTGGCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCAGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTCTCAGGTGATGGCCCCGGAACGCTTTGGc                                                                                >  2:764/1‑251 (MQ=255)                                                                                                                                                                        cgcCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCgtggt                                                                           >  2:44338/1‑251 (MQ=255)                                                                                                                                                                                gggCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTg                                                                   >  1:53081/1‑251 (MQ=255)                                                                                                                                                                                   cGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTgacg                                                                >  1:114787/1‑251 (MQ=255)                                                                                                                                                                                       ggTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACaa                                                            <  1:82031/251‑1 (MQ=255)                                                                                                                                                                                       ggTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACa                                                             <  2:73134/250‑1 (MQ=255)                                                                                                                                                                                       ggTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTGGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACaa                                                            <  1:18079/251‑1 (MQ=255)                                                                                                                                                                                              aGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCcgc                                                     <  1:56805/251‑1 (MQ=255)                                                                                                                                                                                               gTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCAGCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCTGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCcgcg                                                    <  2:180795/251‑1 (MQ=255)                                                                                                                                                                                               gTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCACGGACGGTCTGGCGCAGGCTGTTATCATCGGTGAAACCTTACTTAATGGTGAACCTTCGTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGTCTTCAGTCCGAAGCTGCGTCATGGTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGATTTTGACGACAACTTCcgcg                                                    <  2:34222/251‑1 (MQ=255)                                                                                                                                                                                                                       cGAACAGCCCCGCCCGGCCGGTCTGGCGAAGCCTTTTATCGTCGGTGGCACCTTCCTTCATGGTGAACCTTATTGTCTGGAGCTAGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAACTGGCCACGGTTTTTGGCTATAAGGGGGTGGAATCGGAACGCTTTGGCGTGGTAGGGTTTGACGACAACTTCCGCGCTATCTCGATGGAAGAAAAGCCa                             >  2:98722/1‑250 (MQ=255)                                                                                                                                                                                                                          aCAGCCCAGCCCGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTGAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTGCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAAGTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAAc                         <  1:200292/251‑1 (MQ=255)                                                                                                                                                                                                                                 agcccGGACGTGCTGGCGCAGGCTTTTATCAGCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGGCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGTCTTCAGTCCGAAGCTGCGTCATGTGGCGTCGCGCACCGAAGGGGCGACGGTTTTTGGCTAGCAGGGGAGGGACCCGGAACGGTTTGGCGGGGTGGAGTTTGACGACAATTTCCGCGCTAGCTCGCTGTAAGAAAATCCATAACAGCCGaa                  <  2:124938/251‑1 (MQ=255)                                                                                                                                                                                                                                 agcccGGACGGTCTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGaa                  <  1:69360/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggcTGGCGCAGGCTTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTTTTGTCTGGGGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTAGCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTgg        <  1:241729/249‑1 (MQ=255)                                                                                                                                                                                                                                                 cgcAGGCTTTTATCATCGTTGAAACCTTCCTTAATGGTGAACCTTCTGGTCTGGTGCTGGGCGATAACATCTGCTTCGGTCAGGGCTTCAGTCCGAAGCTTCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTGGACGACAATTTCCTCGCTATCGCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTg  <  1:162066/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CGGGCGGTTCCGGCACCCGATTGCATCCGATTACGCGCGGCGTATCGAAGCAACTGTTGCCGATTTACGATAAGCCAATGATTTACTATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTG  >  NC_000913/3973627‑3974117

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