Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,062,856	G→T	7.3%	G591W (GGG→TGG)	PP_1838 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,062,856	0	G	T	7.3%	185.6 / 4.0	63	G591W (GGG→TGG)	PP_1838	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base T (2/3);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

CCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTAC  >  NC_002947/2062708‑2062999                                                                                                                                                     |                                                                                                                                                ccTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTggg                                                                                                                                                 <  8:123489/149‑1 (MQ=255)  cTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTAGGTGGCAAGGCTCGCCACCAGTGCTggg                                                                                                                                                 <  6:348353/148‑1 (MQ=255)     ggCCGGTTCAACGTGGCGCTGGTGCTGATTGCCCCGGGCATTGAGGAGAAGTTCGGGTCGGTGTATCACACCGGGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGCACTGGTGgcg                                                                                                                                             <  3:37191/149‑1 (MQ=255)         gCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                                                                                                                                         <  6:176557/149‑1 (MQ=255)               aCGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTc                                                                                                                                   <  6:23317/149‑1 (MQ=255)                cGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTc                                                                                                                                   >  7:214111/1‑148 (MQ=255)                      gctgctgctGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCACCCTGc                                                                                                                            >  6:335471/1‑149 (MQ=255)                          ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCGCGCCCCCAGCACTGGTGgcgagcctt                                                                                                                                       <  1:237640/134‑7 (MQ=255)                          ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCGCGCCCCCAGCACTGGTGgcgagcctt                                                                                                                                       >  2:237640/1‑128 (MQ=255)                             ctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAgg                                                                                                                      <  4:156314/148‑1 (MQ=255)                             ctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAgg                                                                                                                      <  3:293681/148‑1 (MQ=255)                              tgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGa                                                                                                                        >  7:313541/1‑145 (MQ=255)                              tgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGa                                                                                                                        <  8:313541/145‑1 (MQ=255)                                aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGa                                                                                                                  >  6:14492/1‑149 (MQ=255)                                aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGa                                                                                                                  >  5:355925/1‑149 (MQ=255)                                aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGa                                                                                                                  >  4:313966/1‑149 (MQ=255)                                     cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGCACTGGTGgcga                                                                                                                                            >  7:55730/1‑117 (MQ=255)                                     cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGCACTGGTGgcga                                                                                                                                            <  8:55730/118‑2 (MQ=255)                                           cATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                       <  6:194652/149‑1 (MQ=255)                                             ttCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGc                                                                                                     <  2:203473/149‑1 (MQ=255)                                                aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCAt                                                                                                   >  6:36076/1‑148 (MQ=255)                                                  gagaAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGAt                                                                                                >  1:33355/1‑149 (MQ=255)                                                     aaGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCaa                                                                                             <  1:300532/149‑1 (MQ=255)                                                      aGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                                            <  3:313966/149‑1 (MQ=255)                                                                ggtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCggg                                                                                    >  7:35304/1‑147 (MQ=255)                                                                ggtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCa                                                                                  <  4:304589/149‑1 (MQ=255)                                                                 gtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAg                                                                                 >  4:192109/1‑149 (MQ=255)                                                                      tCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAg                                                                             >  7:336405/1‑148 (MQ=255)                                                                      tCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAg                                                                             >  8:313521/1‑148 (MQ=255)                                                                              tGGGTACCCAGGTCGACATCGTTCGGAGCATCATGGGCCGCCTGGGTGGCAAGGCTGGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTc                                                                    <  6:299902/149‑1 (MQ=255)                                                                              tGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTc                                                                    <  3:72515/149‑1 (MQ=255)                                                                                ggTACCCAGGTCGACATCGTGCCGACCATCCTGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCgg                                                                  <  5:14492/149‑1 (MQ=255)                                                                                ggTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCgg                                                                  <  2:188793/149‑1 (MQ=255)                                                                                    cccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTg                                                              <  2:355423/149‑1 (MQ=255)                                                                                     ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTgg                                                             >  4:75690/1‑149 (MQ=255)                                                                                     ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATAAAGCCGTCGGGCAGCGAGCAGATTGTCGGCCTgg                                                             >  8:342727/1‑149 (MQ=255)                                                                                       aGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTg                                                           >  5:222567/1‑149 (MQ=255)                                                                                         gTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTg                                                           >  6:245133/1‑147 (MQ=255)                                                                                          tCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAg                                                        <  5:36076/149‑1 (MQ=255)                                                                                            gACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                       <  8:214111/148‑1 (MQ=255)                                                                                                tCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAc                                                  <  6:355925/149‑1 (MQ=255)                                                                                                tCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAc                                                  <  3:192109/149‑1 (MQ=255)                                                                                                 cGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCAGGAGCAGGGTGAcc                                                 >  6:209556/1‑149 (MQ=255)                                                                                                  gTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTTGGGGCGCGACCTGCTCAACCTTCCCGAGGGTGACCAGGGCGTAGGCATGAACAAGCCGTCGGGCAGCGTGCTGATTGTCGGGCTGGTGCAGGGTGAcc                                                 <  8:299287/148‑1 (MQ=255)                                                                                                   tGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGt                                                                     >  1:253189/1‑127 (MQ=255)                                                                                                   tGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGt                                                                     <  2:253189/127‑1 (MQ=255)                                                                                                    gCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTa                                              >  2:290274/1‑149 (MQ=255)                                                                                                     ccGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTa                                              <  1:290274/148‑1 (MQ=255)                                                                                                          catcatGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTAt                                             <  3:289765/144‑1 (MQ=255)                                                                                                          catcatGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTAt                                             >  4:289765/1‑144 (MQ=255)                                                                                                              atGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATc                                      >  6:114173/1‑147 (MQ=255)                                                                                                                  gcctgcTGGGTGGCAAGGCTCGCCACGAGTGCTGGGGGCGCGACCTTCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCTGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTcc                                <  7:342727/144‑1 (MQ=255)                                                                                                                      ccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAgg                            >  3:342041/1‑149 (MQ=255)                                                                                                                       cTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGa                           <  7:313521/149‑1 (MQ=255)                                                                                                                       cTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGa                           >  7:220249/1‑149 (MQ=255)                                                                                                                         gggTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGCCCGTAACCTGATCGAGTCCAAGGACa                         >  3:264627/1‑149 (MQ=255)                                                                                                                          ggTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGAc                          <  8:40367/147‑1 (MQ=255)                                                                                                                          ggTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACAt                        <  5:114173/149‑1 (MQ=255)                                                                                                                          ggTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACAt                        <  8:35304/149‑1 (MQ=255)                                                                                                                                aaGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGccc                  >  7:330967/1‑149 (MQ=255)                                                                                                                                       gCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGa             <  2:33355/147‑1 (MQ=255)                                                                                                                                           ccaGGGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGGGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTAcc       <  4:186901/149‑1 (MQ=255)                                                                                                                                               tGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTa   >  2:357646/1‑149 (MQ=255)                                                                                                                                                 cTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTAc  <  1:357646/148‑1 (MQ=255)                                                                                                                                                 cTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTAc  <  5:245133/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACCGCTAC  >  NC_002947/2062708‑2062999

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