Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,174,041	A→T	E320D (GAA→GAT)	birA →	DNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,174,041	0	A	T	100.0%	160.2 / NA	42	E320D (GAA→GAT)	birA	DNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (24/18);  total (24/18)

CGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGG  >  NC_000913/4173797‑4174286                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cggcCAAGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGataaataa                                                                                                                                                                                                                                                 <  2:151991/251‑1 (MQ=255)   gcCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGa                                                                                                                                                                                                                                               <  2:568155/251‑1 (MQ=255)            aTACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGt                                                                                                                                                                                                                                      >  1:592595/1‑251 (MQ=255)                 tGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTc                                                                                                                                                                                                                                  <  1:266054/250‑1 (MQ=255)                 ngAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTc                                                                                                                                                                                                                                  >  2:266054/2‑250 (MQ=255)                  gAATTACGTGCTGCGTTGGCACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGc                                                                                                                                                                                                                                >  1:347504/1‑251 (MQ=255)                              gcGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAat                                                                                                                                                                                                                    >  2:43525/1‑251 (MQ=255)                                     aaCTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGGATTCGCTCCCCTGCAAATTATTTGc                                                                                                                                                                                                             <  2:319333/251‑1 (MQ=255)                                         ctTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGt                                                                                                                                                                                                           <  1:347173/249‑1 (MQ=255)                                          ttCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCATGGATGGGCGGTGAAATATCCCTGCGTAGCGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGt                                                                                                                                                                                                        >  2:159460/1‑251 (MQ=255)                                               aCAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTACAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGAc                                                                                                                                                                                                   >  1:286721/1‑251 (MQ=255)                                                 aagaagGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACct                                                                                                                                                                                                 >  1:210704/1‑251 (MQ=255)                                                              aCCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATg                                                                                                                                                                                    >  1:125882/1‑251 (MQ=255)                                                                  tATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTCCTACCGCATGAt                                                                                                                                                                                  <  1:427118/249‑1 (MQ=255)                                                                       gTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGc                                                                                                                                                                              <  2:347504/248‑1 (MQ=255)                                                                          gcgcTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTACAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACtttt                                                                                                                                                                        <  2:286721/251‑1 (MQ=255)                                                                                 gAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAgg                                                                                                                                                                 >  2:609616/1‑251 (MQ=255)                                                                                   aaaGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGAt                                                                                                                                                               >  2:96297/1‑251 (MQ=255)                                                                                    aaGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAAt                                                                                                                                                                                                                                                   <  1:429618/166‑1 (MQ=255)                                                                                    aaGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAAt                                                                                                                                                                                                                                                   >  2:429618/1‑166 (MQ=255)                                                                                    aaGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGAtt                                                                                                                                                              >  1:453906/1‑251 (MQ=255)                                                                                      gCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTa                                                                                                                                                             >  1:511310/1‑250 (MQ=255)                                                                                             aaTTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCACTCCCCTGCAAATGATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCTCTTTTCGTCAGTATTAAACTGgcg                                                                                                                                                     <  2:322283/251‑1 (MQ=255)                                                                                                         cGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgt                                                                                                                                          <  2:440475/250‑1 (MQ=255)                                                                                                          gCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtag                                                                                                                                        <  1:448659/251‑1 (MQ=255)                                                                                                               gTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGa                                                                                                                                   <  2:592595/251‑1 (MQ=255)                                                                                                                                               tGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtctc                                                                                                   >  2:226202/1‑251 (MQ=255)                                                                                                                                                 gCATTTCACGCGTAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACGGCATGATTCGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCCGACAGTTGATCTCtt                                                                                                 >  1:470446/1‑251 (MQ=255)                                                                                                                                                         cgcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAAt                                                                                         <  1:497586/251‑1 (MQ=255)                                                                                                                                                         cgcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAAt                                                                                         >  1:344707/1‑251 (MQ=255)                                                                                                                                                           cgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGt                                                                                       >  1:319990/1‑251 (MQ=255)                                                                                                                                                                                   attaCTTGTGCAGGATGGAATAATAAAACACTGGATGGGCGGTGAAATATCCCTGCGTAGAGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCCCGAGTAGGTAGAATATTTTGCTTTAAGATCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTACTCGcttcttc                                                               >  2:328884/1‑251 (MQ=255)                                                                                                                                                                                      aCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCttt                                                            >  1:266137/1‑251 (MQ=255)                                                                                                                                                                                      aCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCttt                                                            <  2:453906/251‑1 (MQ=255)                                                                                                                                                                                               ggATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAAttt                                                   <  2:210704/251‑1 (MQ=255)                                                                                                                                                                                                 aTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGt                                                                                                                                                                     >  1:86369/1‑135 (MQ=255)                                                                                                                                                                                                 aTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGt                                                                                                                                                                     <  2:86369/135‑1 (MQ=255)                                                                                                                                                                                                   ggAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTcgc                                                <  1:335322/250‑1 (MQ=255)                                                                                                                                                                                                         ataaAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTAACGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGtt                                          >  1:563850/1‑250 (MQ=255)                                                                                                                                                                                                               cccTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGTGTAATCGCTCCCCTGCAAATTATTTGAGTAGTCTGACCTATTATACCGCATGATTAGCACTTTTCGCCGGGATTAAACTGGCGCGCTCCCGAGGCGGTAGATTATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTACTTTCGCGTAGTTTTGaac                                    >  2:87319/1‑249 (MQ=255)                                                                                                                                                                                                                     aTGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGa                              >  1:531127/1‑250 (MQ=255)                                                                                                                                                                                                                                           gTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTATTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCtt       >  2:456577/1‑251 (MQ=255)                                                                                                                                                                                                                                                gCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCgg  <  1:43525/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGG  >  NC_000913/4173797‑4174286

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