Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,009,508	G→A	6.6%	intergenic (‑414/‑1254)	PP_2628 ← / → PP_2631	ABC transporter ATP‑binding protein/cellulose biosynthesis protein BcsF/YhjT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,009,508	0	G	A	6.6%	180.9 / 2.1	60	intergenic (‑414/‑1254)	PP_2628/PP_2631	ABC transporter ATP‑binding protein/cellulose biosynthesis protein BcsF/YhjT
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base A (2/2);  total (31/29)
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.18e-01

ACGCCTCCTTGCTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAATGCGCTGCGCCAGGCATT  >  NC_002947/3009365‑3009650                                                                                                                                                     |                                                                                                                                               aCGCCTCCTTGCTATGCCGCCTGGTAGTCGAAGTCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAAc                                                                                                                                           <  2:202279/149‑1 (MQ=255)       cctTGCTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCaa                                                                                                                                     <  1:345788/149‑1 (MQ=255)       cctTGCTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGCCACTGTATGAGGCCAACTCGc                                                                                                                                       >  2:358586/1‑147 (MQ=255)        ctTGCTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCaa                                                                                                                                     <  8:286015/148‑1 (MQ=255)           gCTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGGGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGc                                                                                                                                   <  2:209127/147‑1 (MQ=255)            cTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCCCTGTATGAGGCCAACTCGCAAGCCAc                                                                                                                                <  1:209414/149‑1 (MQ=255)              aTGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCAccc                                                                                                                              >  3:186862/1‑149 (MQ=255)               tGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCAccc                                                                                                                              <  4:186862/148‑1 (MQ=255)                 ccgccAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGt                                                                                                                            >  5:75835/1‑148 (MQ=255)                   gccAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTAt                                                                                                                         >  2:283319/1‑149 (MQ=255)                    ccAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATc                                                                                                                        <  1:283319/149‑1 (MQ=255)                      aGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTg                                                                                                                      <  4:186525/149‑1 (MQ=255)                       ggTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTgg                                                                                                                     <  4:272225/149‑1 (MQ=255)                         tAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGtt                                                                                                                   <  7:343066/149‑1 (MQ=255)                              gAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGa                                                                                                              >  1:148219/1‑149 (MQ=255)                               aaGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGAc                                                                                                             >  5:30201/1‑149 (MQ=255)                                 ggCATGCCCGATAATGCCAGGGCATCCCTGATCGGCGATGCGCAAGCCATGGTTGTCGCCTGCTGGGTGCGCCTCCCGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGaa                                                                                                                <  1:121882/149‑1 (MQ=255)                                   cATGCGCGCAAATGCCAGCGCAGCCCTGATCGCCGCTGCGCAAGCCATGGATGACG‑‑GGCTGGGTGCG‑CG‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCAACCGTTATCTGGTTGAAGACCTGc                                                                                                         <  2:318570/149‑1 (MQ=255)                                    aTGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGcc                                                                                                        >  6:166670/1‑149 (MQ=255)                                        ccGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACACGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGt                                                                                                    >  7:261637/1‑149 (MQ=255)                                         cGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTc                                                                                                   <  4:30884/149‑1 (MQ=255)                                         cGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTc                                                                                                   >  4:70707/1‑149 (MQ=255)                                          gATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCt                                                                                                  >  5:133094/1‑149 (MQ=255)                                              aTGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGAt                                                                                              <  6:133094/149‑1 (MQ=255)                                                gCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGc                                                                                            >  3:304235/1‑149 (MQ=255)                                                gCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACGGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGc                                                                                            >  4:74508/1‑149 (MQ=255)                                                 ccAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGc                                                                                            <  1:358586/148‑1 (MQ=255)                                                     ggCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGc                                                                                                         >  4:351369/1‑131 (MQ=255)                                                     ggCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGc                                                                                                         <  3:351369/131‑1 (MQ=255)                                                          cccTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CTACGCGACGCGCGGTCCGGCCGATCTGGCACTGTATGAGGCAAACTCGCAAGCAAACAGTTATATGGTTGAAGCCCTGCCTCGTCTTGATGCAAGTGGGCgt                                                                                  >  6:186240/1‑149 (MQ=255)                                                          cccTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑AT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCgt                                                                                  >  6:191859/1‑149 (MQ=255)                                                          cccGGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCGGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCgt                                                                                  <  7:260662/149‑1 (MQ=255)                                                               aTCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGtgtg                                                                               <  6:75835/147‑1 (MQ=255)                                                                       tGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGtgctgc                                                                     <  2:148219/149‑1 (MQ=255)                                                                          gcAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTgg                                                                  <  8:192804/149‑1 (MQ=255)                                                                          gcAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTgg                                                                  >  5:99381/1‑149 (MQ=255)                                                                          gcAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGACCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCgt                                                                                  <  6:47571/133‑1 (MQ=255)                                                                          gcAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGACCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCgt                                                                                  <  6:41945/133‑1 (MQ=255)                                                                          gcAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGACCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCgt                                                                                  >  5:47571/1‑133 (MQ=255)                                                                          gcAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGACCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCgt                                                                                  >  5:41945/1‑133 (MQ=255)                                                                               ccATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCg                                                             <  3:70707/149‑1 (MQ=255)                                                                                cATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACACGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCgg                                                            <  8:261637/149‑1 (MQ=255)                                                                                   ggATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTTCTGGTGCTGCTGGCACCGGcac                                                         >  1:388723/1‑149 (MQ=255)                                                                                       gACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTTTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCAc                                                               >  1:230626/1‑139 (MQ=255)                                                                                       gACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCAc                                                               <  2:230626/139‑1 (MQ=255)                                                                                          g‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTg                                                  <  5:166670/149‑1 (MQ=255)                                                                                                         cT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCg                                       >  5:392407/1‑148 (MQ=255)                                                                                                          t‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCt                                                                                                           <  5:112270/79‑1 (MQ=255)                                                                                                          t‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCt                                                                                                           >  6:112270/1‑79 (MQ=255)                                                                                                          t‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCg                                       >  3:174860/1‑147 (MQ=255)                                                                                                               acgcgacgcgCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTc                                  >  8:69273/1‑149 (MQ=255)                                                                                                                  cgacgcgCGGTCCGGCCGAGCTGGCACTGAATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAg                               >  6:152509/1‑149 (MQ=255)                                                                                                                               ggCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAAtg                  >  8:119518/1‑149 (MQ=255)                                                                                                                                gCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGaa                                     >  3:159276/1‑129 (MQ=255)                                                                                                                                gCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAAtgc                 <  6:392407/149‑1 (MQ=255)                                                                                                                                        ggCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCAAACGCCTGGGCGGATACCGAAGTCGACCACTGGAGCAAtgcgcagcgc          >  4:303841/1‑148 (MQ=255)                                                                                                                                          cACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAATGCGCTGCGCCAg       >  1:86953/1‑149 (MQ=255)                                                                                                                                           tcTGTATGAGGCCTACTGGCATGCCAGCCGTTATCTGGTTGATGACCTGCGTCGTTTTTATGCCCGTGGGCGTGTGCTGGTGGTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAATGCGCTGCGCCa        <  8:166081/146‑1 (MQ=255)                                                                                                                                           aCTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAATGCGCTGCGCCa        >  7:166081/1‑147 (MQ=255)                                                                                                                                               tATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAATGCGCTGCGCCAGGCAtt  >  8:309090/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               ACGCCTCCTTGCTATGCCGCCAGGTAGTCGAAGGCATGCCCGATAATGCCAGGGCAGCCCTGATCGCCGATGCGCAAGCCATGGATGACG‑‑TGCTGGGTGCG‑CT‑‑CGACGCGACGCGCGGTCCGGCCGAGCTGGCACTGTATGAGGCCAACTCGCAAGCCACCCGTTATCTGGTTGAAGACCTGCCTCGTCTTGATGCCCGTGGGCGTGTGCTGGTGCTGCTGGCACCGGCACACGCCTGGGCGGATACCGAAGTCGAGCACTGGTGCAATGCGCTGCGCCAGGCATT  >  NC_002947/3009365‑3009650

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: ‑