Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,145,313	C→A	S683Y (TCC→TAC)	yegE →	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,145,313	0	C	A	96.7%	88.9 / NA	31	S683Y (TCC→TAC)	yegE	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (13/16);  minor base G (0/1);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAG  >  NC_000913/2145022‑2145576                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                        agaTGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGCACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTTCATTACAGTATTACGCCGTCAAGTACTCTGGACAGCAGCAATCTTGGTTCGGTTCTGCTGATTCATGCAGTTCCCGACTCCCGCACATCGCTGTCCCCGCTTAGCTACACTCCCACCCattatg                                                                                                                                                                                                                                                                 >  1:207547/1‑300 (MQ=255)                  cAGGAAGAAGGGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCACACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGcc                                                                                                                                                                                                                                                <  2:54490/300‑1 (MQ=255)                  aaGGAAGAAGCGTTAGGTGTTCCGCTCCGGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGcc                                                                                                                                                                                                                                                <  2:277408/299‑1 (MQ=255)                      aagaatTGTTTGGTGGTCCGGTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGTCACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTGCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACcg                                                                                                                                                                                                                                           <  1:55600/301‑1 (MQ=255)                               tccggcgtCCGTTCCGGACTGTGGTTCGTATGTCTTTTGGCGACAACGGACAATTACTGGTGAACCTTTAAAGTGCCGACCCCTCGCGTTCCGAGTTGGGGAGAGATGGTGTGGTGGGCTGGCGGGGCGGCGGCAGTTACGACGTGCGTTACAGTATTACGCCGTTAAGGACTCTGTACAGCAGCTATATGGTTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTAGCATGGTGCACTGACGCATCTCGCCAACCGCGCCAGtt                                                                                                                                                                                                                                   <  1:196393/293‑1 (MQ=255)                                              tGACGGTGTTGCACATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCg                                                                                                                                                                                                                   <  2:218469/301‑1 (MQ=255)                                                       ttcttatTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGACACAACTGCGTATCctgct                                                                                                                                                                                                          <  2:14346/297‑1 (MQ=255)                                                         attatTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTACTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCtgctgc                                                                                                                                                                                                        <  2:142814/299‑1 (MQ=255)                                                           tatTACTTTTGGCGACAACGGATCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGTTTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGTTGAGCTACAGCGCCTACCCAGATGCACCGACGCATCCCGACAACCGCGCCAGTTTTGAGAAACAACTGCGTAACctgacagaa                                                                                                                                                                                                      >  1:14345/1‑295 (MQ=255)                                                             ttACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAc                                                                                                                                                                                                    <  2:119016/301‑1 (MQ=255)                                                             tctcTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAAc                                                                                                                                                                                                    <  2:315103/298‑1 (MQ=255)                                                                cTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGATAAAGCGCCTACCATGATGCACTGACGCATCTCCCCAACCGCGCGCGATTTCAGAAACAACTGCGTTCCctgctcccatcaca                                                                                                                                                                                                 >  1:315095/1‑292 (MQ=255)                                                                                              aCAATTACAGTGACGACACCTCACGGTCCGCGATGGAGAAAGATGTGGTATTGCACTGTCGGAGCGGAGGCAGTTACGACGTGGATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACACATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGgt                                                                                                                                                                   <  1:30560/301‑1 (MQ=255)                                                                                                            gACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCCGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAACCAGTACACATCAGCGCCATGCCCTGGTGTTTATCGacatcg                                                                                                                                                     >  1:10196/1‑296 (MQ=255)                                                                                                                         ccGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGCCGGCAGCAATATTGGTTCGGTTCTGGTGATGCAGGACGTTACCGAATCACGCAAAATGCTGCGCCCGCTGAGCTACAGCGCCTACTATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCCGCTGCAAACGGCAAACAGTACACACCAGCGACATGCCCTGGTGTTTATCGATCTTGACCGCTTTaacg                                                                                                                                         >  1:161098/1‑298 (MQ=255)                                                                                                                                        aTGTGGTGTTGCACTGTCGCAGCGGCGGTAGTTACGACGTGCATTACCGTATTACGCCGTTCGCTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCCCGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCCACCATGATGCACTGACGCCTCTGCCCAACCGCGCCAGCTTTGACACACAACTGCGTATCCTGCTGCAAGCGGTAAACAGCACACAACAGCGACATGCCCGGGTGTTTGTTGGTTTTGATCGCTTTaagaggtggaagtaatgcg                                                                                                                          >  1:163286/1‑283 (MQ=255)                                                                                                                                             ctgtTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCGAGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGc                                                                                                                    <  2:219548/300‑1 (MQ=255)                                                                                                                                               gtTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCCCTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAGACAACTGCGTATCCTTCTGCCAACGGTAAACAGTACACATCAGCGACATTCCATGGTGGTTATCGATCTGGACCGCTGTAACGCGGTGATTGACAGCTCCGCGCag                                                                                                                  >  1:219541/1‑300 (MQ=255)                                                                                                                                                          ggAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCgcatgcccggcgctgt                                                                                                        <  2:50755/300‑16 (MQ=255)                                                                                                                                                             gcggcgGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAACCAGTACACATCAGCGCCATGCCCTGGCGTTTATCGATCTTGATCGCTTCAACGCGGTCAATGAAAGCGCCGGGCATGCGGCGGGttgcg                                                                                                     >  1:90800/1‑296 (MQ=255)                                                                                                                                                             gcggcgGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTTCACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACg                                                                                                     >  2:59680/1‑300 (MQ=255)                                                                                                                                                                         acgacgTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTCAGAAACAACTGCGTACCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGCTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGCGCCTGCGGCGGGTGCCGCTTTGCTGcccg                                                                                         >  2:12841/1‑300 (MQ=255)                                                                                                                                                                                    tACAGTATTACGCCGTTAAGTACGCTGTACGGCAGCTATATTGGGTCGGTTCTGGGGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTTCACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGcgtcgt                                                                             <  1:59678/301‑1 (MQ=255)                                                                                                                                                                                                       gTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGTGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCCGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGCACTGGCGTCGTTAATGCTGAGTATGCTGCg                                                          >  2:159381/1‑301 (MQ=255)                                                                                                                                                                                                             cgatctgcagcaATATTGGTTCGGTTCTTGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTACAg                                                    <  2:188322/295‑1 (MQ=255)                                                                                                                                                                                                               gACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCg                                                  >  2:135198/1‑301 (MQ=255)                                                                                                                                                                                                                   gcagcaATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGAGGCGGGTGACGCTTTGCTGCGCGAACTGGCGCCGGTAATGCTGATGATGATGCGTTcacgaacttg                                              >  1:188316/1‑292 (MQ=255)                                                                                                                                                                                                                           gggggtggggttcGGGGGAGGCAGGACGTTACGGAATCACGCGAAGTGGGGGGCCGGCTGGGCTACAGGGCCTACCATGAGGCACTGAGGCATCGCGGCAACGGCGCCAGTTTTGAGAAAGAACTGCGGATCCTGCTGGAAACGGTAGACAGTACACATCAGGGGGATGCGCGGGTGTTTATCGATCTTGATCGCGTTAAAGCGGGGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCGGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGgcg                                       <  1:12840/292‑1 (MQ=255)                                                                                                                                                                                                                           attggtagagtgctagggagtcgaggcgggaCCGAAGCACGCAAAATGCTGCGCAAGCTGAGCTACGGCGCCTGCCATGATGCACTGACGAATCGCGCCAACCGCACCAGTTATGAGAAACGACTGCGTATCCTGCTGCAAACGGTAAACAGTACACAGCAGCGACAAGCCCTGGGGTTTATCGATCTTGGTCGCTGTAAAGCGGTAAATGACAGCGCCGGGCATGCGGCAGGGGACGCTTTGCGGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTGCCAGCGACATGCTGgcgc                                      <  1:159378/271‑1 (MQ=255)                                                                                                                                                                                                                                ttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCACGGTGCACTGACGAATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACt                                  <  2:315057/300‑1 (MQ=255)                                                                                                                                                                                                                                                                cACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGATGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCTCGCGACGTGCTGCCGCGACTCGGTGGGGATGAATCTGGTCTTCTCCTAcctg  >  2:60006/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                        AGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAG  >  NC_000913/2145022‑2145576

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