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	98.0%	184.9 / NA	49	E320D (GAA→GAT)	birA	DNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/24);  minor base G (1/0);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGGA  >  NC_000913/4173798‑4174287                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ggcCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAg                                                                                                                                                                                                                                                 <  2:490388/251‑1 (MQ=255)             aCGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTAt                                                                                                                                                                                                                                     >  2:132818/1‑251 (MQ=255)                                    aaCTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGc                                                                                                                                                                                                              <  2:275139/251‑1 (MQ=255)                                               cAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGAcc                                                                                                                                                                                                   <  2:432450/251‑1 (MQ=255)                                                      gATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTcccc                                                                                                                                                                                                                            <  1:177604/219‑1 (MQ=255)                                                      gATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTcccc                                                                                                                                                                                                                            >  2:177604/1‑219 (MQ=255)                                                         tGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCtctt                                                                                                                                                                                               <  1:219804/245‑1 (MQ=255)                                                         tGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCtctt                                                                                                                                                                                               >  2:219804/1‑245 (MQ=255)                                                                  aTCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGCAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGAtt                                                                                                                                                                                  >  1:84526/1‑249 (MQ=255)                                                                    cTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGc                                                                                                                                                                               >  2:475225/1‑250 (MQ=255)                                                                           gcTGGGAAAAGCTGGATAATTGTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCg                                                                                                                                                                       <  1:232580/251‑1 (MQ=255)                                                                              gggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACctc                                                                                                                                                                                                 <  1:102882/222‑1 (MQ=255)                                                                              gggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTc                                                                                                                                                                     >  1:248080/1‑250 (MQ=255)                                                                              gggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACACTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGATCCCCTGCAAATTATTTGCGTAGTCTGACctc                                                                                                                                                                                                 >  2:102882/1‑222 (MQ=255)                                                                                      cTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATCAAGAAATATTTGGCATTTCAAGCGGAATTGACAAAAAGGGGGCTTTCTTACTAGAGAAGGATGGAATCATAAAACCCTGGCTGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCACTGCAAATTATTTGAGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTaa                                                                                                                                                             >  1:282804/1‑250 (MQ=255)                                                                                                  attaATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTc                                                                                                                                                 >  1:264969/1‑250 (MQ=255)                                                                                                    ttaTCGCCCAGTGAAACCTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATTGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCAc                                                                                                                                               <  2:439095/248‑1 (MQ=255)                                                                                                        cGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGGAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgta                                                                                                                                          <  2:264969/251‑1 (MQ=255)                                                                                                           ccAGTGAAACTTATAATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggg                                                                                                                                       >  2:124788/1‑250 (MQ=255)                                                                                                              gTGAAACTTATAATAGGTGATACAGAAATATTTGGCATTTCACGCGGTCTAGACCAACAGGGGGCTTTATAACTTGAGCAGGATGGAATAATAAAAACCTGGATGGGCGATGAAATATCCCTGCGTAGTGCAGAGAACTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTCGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGCgtaggt                                                                                                                                       >  2:252629/1‑248 (MQ=255)                                                                                                                    cTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCa                                                                                                                                                                    >  2:144639/1‑213 (MQ=255)                                                                                                                    cTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACGTTTCGTCa                                                                                                                                                                    <  1:144639/213‑1 (MQ=255)                                                                                                                                      gAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGcc                                                                                                              >  2:326951/1‑249 (MQ=255)                                                                                                                                               ggCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATctct                                                                                                   >  2:309114/1‑251 (MQ=255)                                                                                                                                                gCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCtt                                                                                                  <  2:335822/251‑1 (MQ=255)                                                                                                                                                         gcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACt                                                                                                                                                                            <  1:124549/168‑1 (MQ=255)                                                                                                                                                         gcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACt                                                                                                                                                                            >  2:124549/1‑168 (MQ=255)                                                                                                                                                            gAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTGCCACAATGTCa                                                                                      <  1:132818/251‑1 (MQ=255)                                                                                                                                                                  aCAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATGTCTGTCCACAAGGTCATGGCAg                                                                                <  1:475225/251‑1 (MQ=255)                                                                                                                                                                    aaaCAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTc                                                                              >  1:239477/1‑251 (MQ=255)                                                                                                                                                                     aaCAGGGGTCTTTATTACGTTAGCTGGATGGAATAATAAAACCCTGGATGGGCGGTTAAATATCCCTGCGTAGTTCAGATAAATAAGAAAGGTGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTCCTACCGCATTATTAGCACTGTTCGGCATGATTAAACTTGCGCGCTCACGAGTAAGTAGAATATTTTGCGTTAAGTGCAGCCAGTTGATCTCGTTCCACAATGTCAGGGGAGTCt                                                                             <  2:267562/251‑1 (MQ=255)                                                                                                                                                                             gCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGct                                                                     >  1:187134/1‑251 (MQ=255)                                                                                                                                                                                 tattaCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAACTATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGcttctt                                                                 >  2:536693/1‑251 (MQ=255)                                                                                                                                                                                  attaCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCTCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTGTCCACAATGTCATGGCAGTCTTAATCGcttcttc                                                                <  2:84526/251‑1 (MQ=255)                                                                                                                                                                                     aCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCttt                                                             >  2:39304/1‑251 (MQ=255)                                                                                                                                                                                          atcAGGATGGAATAATAAAACCCGGGATGTGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCGACCGCATGATTAGCACTTTTCGTCAGGATTAAACTTGCGCGCTCACGAGTAGGTAGAATATTTTTCTTTAATTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTa                                                        <  1:610661/249‑1 (MQ=255)                                                                                                                                                                                            cAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAAt                                                      >  2:232734/1‑251 (MQ=255)                                                                                                                                                                                               gATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTc                                                   <  1:225932/251‑1 (MQ=255)                                                                                                                                                                                                          aaaaCCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATg                                        <  1:526493/251‑1 (MQ=255)                                                                                                                                                                                                           aaaCCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGa                                       <  2:124634/251‑1 (MQ=255)                                                                                                                                                                                                                 tGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAAt                                   >  2:575721/1‑249 (MQ=255)                                                                                                                                                                                                                       ggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATcc                           <  2:137039/251‑1 (MQ=255)                                                                                                                                                                                                                       ggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGATCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATcc                           >  1:42543/1‑251 (MQ=255)                                                                                                                                                                                                                           gTGAAATATCCCTGCGTTGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGGCTTAAACGCGTCTTCTCTAGGTAATTTCGCGTAGTTATGAAAATAGGAATCCgggg                       <  2:71593/251‑2 (MQ=255)                                                                                                                                                                                                                                 tatCCCTGCGTAGTGAAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAATATAGGAATCCGGGTCGGTa                  >  2:79758/1‑250 (MQ=255)                                                                                                                                                                                                                                         cGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCCAATTATTTGAGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCt         >  1:41692/1‑251 (MQ=255)                                                                                                                                                                                                                                            aGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTcg      <  2:187134/251‑1 (MQ=255)                                                                                                                                                                                                                                               gCAGATAAATAAGAAAGGGGAGTATTAGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCgg   >  2:284024/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGAGTAGTGATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGGa  <  2:12029/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCGCGGA  >  NC_000913/4173798‑4174287

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