Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,528,718	G→A	11.0%	V353V (GTG→GTA)	ftsA →	cell division protein FtsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,528,718	0	G	A	11.0%	165.2 / 9.5	64	V353V (GTG→GTA)	ftsA	cell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base A (3/4);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

TCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGAT  >  NC_002947/1528570‑1528865                                                                                                                                                     |                                                                                                                                                    tctTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTg                                                                                                                                                     >  7:181243/1‑149 (MQ=255) tctTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTCCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTg                                                                                                                                                     >  7:4663/1‑149 (MQ=255)      aCCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGca                                                                                                                                                >  5:16559/1‑149 (MQ=255)        ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                              <  3:276843/149‑1 (MQ=255)        ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCgg                                                                                                                                                 >  1:121038/1‑145 (MQ=255)        ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCgg                                                                                                                                                 <  2:121038/146‑2 (MQ=255)        cattATCCAGGACGAACTGCGCCGCAGTGGCTACGAAGAAATGGTGCCCGCCGGCATCGTCATCACCGGCGGCACCGCGAAGATGGAAGGCCCCGTCGAACTGCACGAGGAAATCTTCAACATGCCGCTACGCCTGCGAGTGCCGca                                                                                                                                                >  4:166110/3‑147 (MQ=255)          ttATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                            >  4:226223/1‑149 (MQ=255)                aGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCggc                                                                                                                                                <  4:82630/139‑3 (MQ=255)                aGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCggc                                                                                                                                                >  3:82630/1‑137 (MQ=255)                 ggCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCggg                                                                                                                                     >  8:255135/1‑149 (MQ=255)                 ggCCGAACTGCGCCGCAGTGGCGACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGGGCCGCACAGCGTTCggg                                                                                                                                     <  6:324993/149‑1 (MQ=255)                  gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                              >  6:32736/1‑139 (MQ=255)                  gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                              <  5:32736/139‑1 (MQ=255)                    cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                            >  7:273720/1‑139 (MQ=255)                    cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                            <  8:273720/139‑1 (MQ=255)                             ccgcAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACg                                                                                                                         >  7:284141/1‑149 (MQ=255)                                  gTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtggtg                                                                                                                    <  8:181243/149‑1 (MQ=255)                                   tGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                                   <  8:112741/149‑1 (MQ=255)                                    ggCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                                   <  3:106971/148‑1 (MQ=255)                                       tACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGCCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCaa                                                                                                               >  6:168449/1‑149 (MQ=255)                                             gACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                                         <  3:2663/149‑1 (MQ=255)                                                 tGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtg                                                                                                                       >  2:28036/1‑131 (MQ=255)                                                 tGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACgtg                                                                                                                       <  1:28036/131‑1 (MQ=255)                                                   gTGCCCGCCGGCATCGGCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTAt                                                                                                     <  5:41287/147‑1 (MQ=255)                                                     gCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCa                                                                                                 <  8:284141/149‑1 (MQ=255)                                                        cgccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTg                                                                                              >  5:82444/1‑149 (MQ=255)                                                        cgccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTg                                                                                              >  1:247152/1‑149 (MQ=255)                                                          ccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGgt                                                                                            <  6:210023/149‑1 (MQ=255)                                                             gCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGtgtg                                                                                         <  1:293789/149‑1 (MQ=255)                                                             gCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGtgtg                                                                                         <  3:226223/149‑1 (MQ=255)                                                                tCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgcg                                                                                                                  <  8:80761/121‑1 (MQ=255)                                                                tCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgcg                                                                                                                  >  7:80761/1‑121 (MQ=255)                                                                  gTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCtt                                                                                    >  8:15779/1‑149 (MQ=255)                                                                      tccccGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTc                                                                                >  5:265366/4‑149 (MQ=255)                                                                        aCCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTc                                                                                <  4:103994/147‑1 (MQ=255)                                                                              ggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCTGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACgg                                                                        <  8:96095/149‑1 (MQ=255)                                                                                 aCCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCct                                                                     <  6:82444/149‑1 (MQ=255)                                                                                  gcgcgAAGATGGAATGCTCCGTCGAACTGGCCGAGGAAATCTTGCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGTGTGAGCGACGTGGTGCGCATCCCGATGTATTCAACTGGTTTGGGCTTGCTCACCTACGGCctg                                                                    <  3:166110/148‑1 (MQ=255)                                                                                   cgcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgc                                                                   >  5:158672/1‑149 (MQ=255)                                                                                    gcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGTAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgct                                                                  >  4:308111/1‑149 (MQ=255)                                                                                      gAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgctg                                                                 >  3:30677/1‑148 (MQ=255)                                                                                      gAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGa                                                                <  4:30677/149‑1 (MQ=255)                                                                                       aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                               <  7:255135/149‑1 (MQ=255)                                                                                             gAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGt                                                          <  2:221569/148‑1 (MQ=255)                                                                                               aGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTcc                                                        >  1:308932/1‑148 (MQ=255)                                                                                                 gcgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAg                                                     <  6:16559/149‑1 (MQ=255)                                                                                                  cgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTTCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAgg                                                    <  3:200335/149‑1 (MQ=255)                                                                                                    ccGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                                  <  7:15779/149‑1 (MQ=255)                                                                                                         gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                               >  5:55947/1‑131 (MQ=255)                                                                                                         gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                               <  6:55947/131‑1 (MQ=255)                                                                                                         gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGc                                             >  8:104669/1‑149 (MQ=255)                                                                                                           aCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGccc                                           <  7:104669/149‑1 (MQ=255)                                                                                                               gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGTAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGa                                        <  3:308111/148‑1 (MQ=255)                                                                                                                 cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                                     >  4:149490/1‑149 (MQ=255)                                                                                                                 cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                                     >  1:311953/1‑149 (MQ=255)                                                                                                                  gAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                                    <  2:308932/149‑1 (MQ=255)                                                                                                                         tCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacag                              >  7:29208/1‑148 (MQ=255)                                                                                                                             ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagc                             <  6:158672/145‑1 (MQ=255)                                                                                                                             ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagc                             <  5:168449/145‑1 (MQ=255)                                                                                                                               acaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                                    <  6:151056/136‑1 (MQ=255)                                                                                                                               acaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                                    >  5:151056/1‑136 (MQ=255)                                                                                                                                ctgtgcGGCACGCCCAGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac                          >  1:97568/6‑145 (MQ=255)                                                                                                                                ctgtgcGGCACGCCCAGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac                          <  2:97568/140‑1 (MQ=255)                                                                                                                                          cGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACa             <  2:311953/148‑1 (MQ=255)                                                                                                                                             cTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTa         >  3:188393/1‑149 (MQ=255)                                                                                                                                                     gCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGAt  >  3:310873/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTATGGCGAT  >  NC_002947/1528570‑1528865

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