Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,222,743	C→T	intergenic (‑256/+61)	arpA ← / ← iclR	regulator of acetyl CoA synthetase/DNA‑binding transcriptional repressor IclR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,222,743	0	C	T	97.2%	109.6 / ‑6.1	36	intergenic (‑256/+61)	arpA/iclR	regulator of acetyl CoA synthetase/DNA‑binding transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/11);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACCCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGTGC  >  NC_000913/4222498‑4222980                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               cATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑At                                                                                                                                                                                                                                          >  1:264013/1‑251 (MQ=255)          ttGATGTCTTGAATTTTAACTCTCAGCACAGTGAAAAAATTTAATATGATTACAAATATAGCAATATCATAAATCGCTCAATCTCATAATGCAGCAGTATAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAGTTAAGCGGATTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGAGAAGTGGATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGAAAATCTT‑ATTCTATTGc                                                                                                                                                                                                                                  >  2:565699/1‑250 (MQ=255)            gATGTCTTGAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCa                                                                                                                                                                                                                                <  2:246152/250‑1 (MQ=255)                    gAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGGGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTa                                                                                                                                                                                                                        <  1:213409/250‑1 (MQ=255)                                               aaTTTAATATGATTACAACTAAAGAAATAACATAAATCGCACAATCTAATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGAGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAATTAATGCTGAGTTCATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCAACTCAGGTATGATTGGAAGAATATTGCCTCTGCCCGc                                                                                                                                                                                            >  2:292311/1‑251 (MQ=255)                                                   tAATATGATTACAACTAGCGAATTATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGCGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGAGCAGAATATTGCCTCTGCCCGcc                                                                                                                                                                                           >  2:289783/1‑248 (MQ=255)                                                      tatGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTACGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTCTGATGGGCAGATTATTGCCTCTGCCCGCCAGaca                                                                                                                                                                                      >  2:56794/1‑248 (MQ=255)                                                        tGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGaaac                                                                                                                                                                                     >  2:164113/1‑248 (MQ=255)                                                             aCAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAgc                                                                                                                                                                              <  1:164113/251‑1 (MQ=255)                                                                     aGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGCGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGCATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTcc                                                                                                                                                                       >  1:50099/1‑250 (MQ=255)                                                                           atCATAAATCGCTCAAGCACATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTa                                                                                                                                                                 <  1:296814/250‑1 (MQ=255)                                                                              aTAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGcc                                                                                                                                                             <  1:371198/251‑1 (MQ=255)                                                                                                 aaTGCAGCCTTAAAAGTTTCGGTGGAATTAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAAGAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTcgccg                                                                                                                                          <  1:463341/251‑1 (MQ=255)                                                                                                  aTGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTcgccgc                                                                                                                                         >  1:153732/1‑251 (MQ=255)                                                                                                   tGCAGCCGTCAAAGTTTCGGAGGAATGGGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGATGAGTTAATAAACGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCAACCGTAAGCCAGCGTCCCTTCCTTCGCCGCt                                                                                                                                        >  2:498832/1‑251 (MQ=255)                                                                                                          gTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGCTTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGAGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCAACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGAGTCCCTTCCTTCGCAGCTTTAATc                                                                                                                                  >  2:383669/1‑250 (MQ=255)                                                                                                          gTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTACTAATTAAGTGGTTTGTTGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCACTTTAATc                                                                                                                                  >  2:119030/1‑250 (MQ=255)                                                                                                                gTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTAGGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATcg                                                                                                                           <  1:486265/251‑1 (MQ=255)                                                                                                                  ttCGGGGGAATGATAT‑CTTGCGATTTTGTTAATAATTAAGGGGTTTGATGTTTAAAAATAATGCGGAGTTAATAAAAGTGTGTGAAGTGTAGGACATTAAGTTATTTGGGCGTCAATGCGATTAACAGAC‑CCTTTGATTCTATGGCCACTCAGGTATGATGGGGAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGGACGCCAGCGTCACTTCCTTCGCGGCTTTAATCACCATcgc                                                                                                                          <  2:264013/250‑1 (MQ=255)                                                                                                                   ttgttGGAATGAGATTCTT‑CGATTGTGTTAATAATTAATTGGTTTTATGTTTAAAATAAATTCTGAGTAAATAAAAGGGTGTGAAGGGTAGGACATTAAGTTATGTTGGCGTCAATGCAATAAACAGACACTCTT‑ATTCTGTTTCCACTCAGGGATGATGGGCAGAATATGGCCTCGGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCGTTCGCCGCTTTAATCGCCATcgcg                                                                                                                         <  2:278662/246‑1 (MQ=255)                                                                                                                                         aTTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCAt                                                                                                   <  1:383669/251‑1 (MQ=255)                                                                                                                                                         aGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAAt                                                                                   >  1:362044/1‑251 (MQ=255)                                                                                                                                                              tttGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCCTCGGTCATACGTGAAATCGGTc                                                                              >  2:462149/1‑251 (MQ=255)                                                                                                                                                                tGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGATATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCATCGTCACTTCCTTCGCCGCTTTAATCCCCATCGCGCCAAACTCGGTCACGCGGTCATCGGTCATACGTGAAATCGGTCCg                                                                            >  2:134383/1‑251 (MQ=255)                                                                                                                                                                                tGCTGAGTTAATAAAAGTGTGTGACGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCgg                                                            >  1:471302/1‑251 (MQ=255)                                                                                                                                                                                             aaaGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCACACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAAAGGTTCACg                                               >  2:262124/1‑251 (MQ=255)                                                                                                                                                                                               aGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACggg                                             >  2:155852/1‑250 (MQ=255)                                                                                                                                                                                                      gAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGTTGCTCATc                                      >  1:300765/1‑251 (MQ=255)                                                                                                                                                                                                             aTGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGCCACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGAAATACGTGAAATCGGTCCGGAAATAGAAATTGCGACAAACGGTTCACGGTACTCATCGAAAAt                                >  2:423196/1‑250 (MQ=255)                                                                                                                                                                                                             aTGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATa                               <  1:155852/251‑1 (MQ=255)                                                                                                                                                                                                                    tAAGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGACATCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTAGGTCACGCGGTCATCGGTAATACGTGAACTCGGTCCGGAACTAGAAATTGCGACAACCGGTTCACGGTGCTCATctaacatcaactctg                        >  2:100412/1‑237 (MQ=255)                                                                                                                                                                                                                      aGTTATTTTGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGc                       >  1:18945/1‑250 (MQ=255)                                                                                                                                                                                                                            tttGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAATGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCg                                                                            >  2:265418/1‑191 (MQ=255)                                                                                                                                                                                                                            tttGGCGTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACCCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATGCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCg                                                                            <  1:265418/191‑1 (MQ=255)                                                                                                                                                                                                                                  gTCAATGCGATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAACCGGTTCACGGTGCTCGTCGAAAATACACGCTGCAAGGCAACGTAGc          >  1:15710/1‑251 (MQ=255)                                                                                                                                                                                                                                          gATTAACAGACACTCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGTGc  >  1:33771/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CATTGTTTATTGATGTCTTGAATTTTAACTATCAGAACAGTGAAAAAATTTAATATGATTACAACTAAAGAAATATCATAAATCGCTCAATCTCATAATGCAGCCGTAAAAGTTTCGGTGGAATGAGAT‑CTTGCGATTTTGTTAATAATTAAGTGGTTTGATGTTTAAAAATAATGCTGAGTTAATAAAAGTGTGTGAAGTGTATGACATTAAGTTATTTTGGCGTCAATGCGATTAACAGACACCCTT‑ATTCTATTGCCACTCAGGTATGATGGGCAGAATATTGCCTCTGCCCGCCAGAAAAAGTCAGCGCATTCCACCGTACGCCAGCGTCACTTCCTTCGCCGCTTTAATCACCATCGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCAGTGC  >  NC_000913/4222498‑4222980

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