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%	161.1 / NA	41	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 (22/19);  total (22/19)

GCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGG  >  NC_000913/4173804‑4174286                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      gCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAggg                                                                                                                                                                                                                                          >  2:23510/1‑251 (MQ=255)         gTGAATTACGTGATGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTc                                                                                                                                                                                                                                  >  2:283458/1‑251 (MQ=255)           gAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGc                                                                                                                                                                                                                                >  1:656905/1‑251 (MQ=255)            aaTTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCt                                                                                                                                                                                                                               >  2:530629/1‑251 (MQ=255)                 cGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCt                                                                                                                                                                                                                          >  1:30808/1‑251 (MQ=255)                       gcGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAat                                                                                                                                                                                                                    >  1:75808/1‑251 (MQ=255)                                  ctTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGt                                                                                                                                                                                                           >  2:502893/1‑249 (MQ=255)                                          aagaagGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACct                                                                                                                                                                                                 >  1:236978/1‑251 (MQ=255)                                           agaagGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACctc                                                                                                                                                                                                >  1:1476/1‑251 (MQ=255)                                                 aTTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTAc                                                                                                                                                                                          <  1:149439/251‑1 (MQ=255)                                                           tATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGAtt                                                                                                                                                                                 >  1:153300/1‑250 (MQ=255)                                                                  cgcgCTGGGAAAAGATGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGCTTAGCACttt                                                                                                                                                                         >  2:63593/1‑251 (MQ=255)                                                                        gggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTc                                                                                                                                                                    <  1:283458/250‑1 (MQ=255)                                                                          gAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAgg                                                                                                                                                                 >  1:28694/1‑251 (MQ=255)                                                                               gCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTaa                                                                                                                                                            <  1:573208/251‑1 (MQ=255)                                                                                       aTTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcgc                                                                                                                                                    >  1:453111/1‑251 (MQ=255)                                                                                               aaTCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAg                                                                                                                                                                                                         >  1:485331/1‑190 (MQ=255)                                                                                               aaTCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAg                                                                                                                                                                                                         <  2:485331/190‑1 (MQ=255)                                                                                                       aGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTTGTCAGGATTAAACTGGCGCGCTCACGAgtaggtag                                                                                                                                    <  2:175592/251‑1 (MQ=255)                                                                                                           aaaCTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAat                                                                                                                                 >  2:28738/1‑250 (MQ=255)                                                                                                                        gTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGt                                                                                                                   >  1:500044/1‑251 (MQ=255)                                                                                                                                        tGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtct                                                                                                    <  1:223610/250‑1 (MQ=255)                                                                                                                                                            aCAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAg                                                                               <  2:453111/251‑1 (MQ=255)                                                                                                                                                             cAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTTATGGCAGt                                                                              <  2:656905/251‑1 (MQ=255)                                                                                                                                                                aCAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCtt                                                                           <  1:23510/251‑1 (MQ=255)                                                                                                                                                                    ggggCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATc                                                                       <  1:182980/251‑1 (MQ=255)                                                                                                                                                                              taCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGcttcttct                                                              <  2:311235/250‑1 (MQ=255)                                                                                                                                                                                      cAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACct                                                                                                                                                                                                 >  1:541202/1‑111 (MQ=255)                                                                                                                                                                                      cAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACct                                                                                                                                                                                                 <  2:541202/111‑1 (MQ=255)                                                                                                                                                                                        ggATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAAttt                                                   <  2:236978/251‑1 (MQ=255)                                                                                                                                                                                          aTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTcg                                                 <  1:412780/251‑1 (MQ=255)                                                                                                                                                                                          aTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTcg                                                 >  1:84443/1‑251 (MQ=255)                                                                                                                                                                                                         ccTGGCTGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGACAGGGGAGTATTCGCACCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGAGTAGTTATGCAAAt                                  >  2:497586/1‑251 (MQ=255)                                                                                                                                                                                                             gATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGa                              <  1:502893/251‑1 (MQ=255)                                                                                                                                                                                                                  gCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCTAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCg                         >  1:690145/1‑251 (MQ=255)                                                                                                                                                                                                                       gAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggta                                                                                                                                     >  1:570940/1‑138 (MQ=255)                                                                                                                                                                                                                       gAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggta                                                                                                                                     <  2:570940/138‑1 (MQ=255)                                                                                                                                                                                                                       gAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTc                     <  1:490066/250‑1 (MQ=255)                                                                                                                                                                                                                         aaTATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGt                  <  2:1476/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGcc           <  1:28738/249‑1 (MQ=255)                                                                                                                                                                                                                                          cAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCgg  >  2:698222/1‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      GCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGG  >  NC_000913/4173804‑4174286

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