Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,473,613	G→A	6.8%	A147A (GCC→GCT)	bdcA ←	c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,473,613	0	G	A	6.8%	190.2 / 3.1	59	A147A (GCC→GCT)	bdcA	c‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base A (2/2);  total (31/28)
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.40e-01

ATGCAACATATCGCGCATTGGCCCGTTGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAA  >  NC_000913/4473469‑4473758                                                                                                                                                 |                                                                                                                                                  aTGCAACATATCGCGCATTGGCCCGTTGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCa                                                                                                                                                >  7:96082/1‑148 (MQ=255)        atatCGCGCATTGGCCCGTTGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAagcag                                                                                                                                         >  5:302629/1‑148 (MQ=255)           tCGCGCATTGGCCCGTTGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       <  6:302629/147‑1 (MQ=255)             gcgcATTGGCCCGTTGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATg                                                                                                                                    <  8:96082/148‑1 (MQ=255)                           tGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATaa                                                                                                                                             >  6:25816/1‑125 (MQ=255)                           tGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATaa                                                                                                                                             <  5:25816/125‑1 (MQ=255)                            ggCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCAt                                                                                                                      <  1:46781/147‑1 (MQ=255)                                    tAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       >  6:201589/1‑122 (MQ=255)                                    tAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       <  5:201589/122‑1 (MQ=255)                                           ggTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACg                                                                                                      <  4:319260/148‑1 (MQ=255)                                             tATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       <  7:386335/113‑1 (MQ=255)                                             tATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       >  8:386335/1‑113 (MQ=255)                                                cAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGcc                                                                                                 >  2:252836/1‑148 (MQ=255)                                                cAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGcc                                                                                                 >  3:87923/1‑148 (MQ=255)                                                        ccTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGAtt                                                                                          >  4:169522/1‑147 (MQ=255)                                                          tggctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       <  6:87509/100‑1 (MQ=255)                                                          tggctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       >  5:87509/1‑100 (MQ=255)                                                          tggctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGCGCCGATGAtt                                                                                          >  1:208624/1‑145 (MQ=255)                                                          tggctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGAtt                                                                                          >  2:117774/1‑145 (MQ=255)                                                           ggctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGAtt                                                                                          >  3:134801/1‑144 (MQ=255)                                                            gctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAg                                                                                       <  1:117774/146‑1 (MQ=255)                                                            gctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAg                                                                                       <  5:220115/146‑1 (MQ=255)                                                            gctggACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAg                                                                                       <  4:87923/146‑1 (MQ=255)                                                                  cgacgTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAg                                                                                       >  7:182618/1‑140 (MQ=255)                                                                  cgacgTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAg                                                                                       <  8:182618/140‑1 (MQ=255)                                                                       ttAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCtt                                                                          >  8:29055/1‑148 (MQ=255)                                                                        tAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTc                                                                         >  4:282394/1‑148 (MQ=255)                                                                        tAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTc                                                                         >  5:255992/1‑148 (MQ=255)                                                                            ggTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCggg                                                                      >  4:202704/1‑147 (MQ=255)                                                                             gTGATCCCCCGCGGGCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCa                                                                    >  6:127912/1‑148 (MQ=255)                                                                                 tCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTg                                                                >  5:250618/1‑148 (MQ=255)                                                                                  cccACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGc                                                               >  4:1554/1‑148 (MQ=255)                                                                                     aCGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTc                                                                                                         >  3:328042/1‑103 (MQ=255)                                                                                     aCGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTc                                                                                                         <  4:328042/103‑1 (MQ=255)                                                                                      cgcgGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCggg                                                           <  1:252836/148‑1 (MQ=255)                                                                                            ccAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       <  8:356599/66‑1 (MQ=255)                                                                                            ccAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGcc                                                                                                                                       >  7:356599/1‑66 (MQ=255)                                                                                            ccAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTaa                                                                                        <  6:184395/113‑1 (MQ=255)                                                                                            ccAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCtt                                                     >  3:187201/1‑148 (MQ=255)                                                                                            ccAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGCGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTaa                                                                                        >  5:184395/1‑113 (MQ=255)                                                                                             cAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCtt                                                     <  4:134801/147‑1 (MQ=255)                                                                                                aaTCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAAc                                                 <  6:250618/148‑1 (MQ=255)                                                                                                           aGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGAt                                      <  3:202704/148‑1 (MQ=255)                                                                                                            gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTGGCGGAGCCGATGATTAAGATGCGCCCGCGTGCGGGCCTCTGCCGGGCGGCTTCAGCTGGGGCATgg                                       >  2:334566/1‑145 (MQ=255)                                                                                                            gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATa                                     <  3:169522/148‑1 (MQ=255)                                                                                                            gCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATa                                     <  1:334566/148‑1 (MQ=255)                                                                                                               cgcgcgCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAg                                   >  2:130392/1‑147 (MQ=255)                                                                                                                 cgcgcCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAg                                >  1:29475/1‑148 (MQ=255)                                                                                                                      cATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGa                           <  2:29475/148‑1 (MQ=255)                                                                                                                        tGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAat                         <  6:255992/148‑1 (MQ=255)                                                                                                                         gCCTTGCAGGGCAGATTTGCTGGCAGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCAt                                                                   >  7:66336/1‑105 (MQ=255)                                                                                                                         gCCTTGCAGGGCAGATTTGCTGGCAGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCAt                                                                   >  7:72640/1‑105 (MQ=255)                                                                                                                         gCCTTGCAGGGCAGATTTGCTGGCAGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCAt                                                                   <  8:66336/105‑1 (MQ=255)                                                                                                                         gCCTTGCAGGGCAGATTTGCTGGCAGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCAt                                                                   <  8:72640/105‑1 (MQ=255)                                                                                                                              gCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCTCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTg                     <  5:127912/146‑1 (MQ=255)                                                                                                                                    cAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGaa             <  1:130392/148‑1 (MQ=255)                                                                                                                                    cAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGa              >  2:267233/1‑147 (MQ=255)                                                                                                                                     aGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGaaa            <  3:282394/148‑1 (MQ=255)                                                                                                                                        tttGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAgg         >  7:124380/1‑148 (MQ=255)                                                                                                                                              ggcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGaa             >  3:144416/1‑138 (MQ=255)                                                                                                                                              ggcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGaa             <  4:144416/138‑1 (MQ=255)                                                                                                                                               gcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCaa  >  2:180156/1‑148 (MQ=255)                                                                                                                                               gcggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCaa  >  3:54222/1‑148 (MQ=255)                                                                                                                                                cggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCaa  >  3:87486/1‑147 (MQ=255)                                                                                                                                                cggcATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGACGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCaa  >  7:18401/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ATGCAACATATCGCGCATTGGCCCGTTGGCGGGATTAGCGTCGGTATCAATTGGCCCTGGCTGGACGACGTTAATGGTGATCCCACGCGGTCCAAAATCACGGGCCAGCCCGCGCGCCATGCCTTGCAGGGCAGATTTGCTGGCGGCATAAGCAGCCATGCCTGCAACAGGCATACGATCGCCATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAA  >  NC_000913/4473469‑4473758

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