Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,062,856	G→T	7.2%	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.2%	201.4 / 3.0	64	G591W (GGG→TGG)	PP_1838	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base T (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACC  >  NC_002947/2062709‑2062994                                                                                                                                                    |                                                                                                                                           cTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTgggg                                                                                                                                           <  1:239367/149‑1 (MQ=255) cTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTgggg                                                                                                                                           >  2:78163/1‑149 (MQ=255)         cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAcc                                                                                                                                   >  1:38957/1‑149 (MQ=255)          atCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCGCGCCCCCAGCACTGGTGgcgagcctt                                                                                                                                  >  6:281146/2‑143 (MQ=255)            cAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCGCGCCCCCAGCACTGGTGgcgagcctt                                                                                                                                  <  5:281146/147‑7 (MQ=255)                gTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGCCATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCaa                                                                                                                            >  4:215362/1‑149 (MQ=255)                 tGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAAc                                                                                                                           <  4:106113/149‑1 (MQ=255)                          tgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGCACTGGTGgcgagc                                                                                                                                     >  5:292582/1‑127 (MQ=255)                          tgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGCACTGGTGgcgagc                                                                                                                                     <  6:292582/130‑4 (MQ=255)                          tgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCTAg                                                                                                                  >  2:262549/1‑149 (MQ=255)                               aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTg                                                                                                                        <  8:333329/138‑1 (MQ=255)                               aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTg                                                                                                                        >  7:333329/1‑138 (MQ=255)                                  gCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGAc                                                                                                            >  5:12630/1‑147 (MQ=255)                                     ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                                       >  8:43252/1‑149 (MQ=255)                                     ccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                                       >  7:102408/1‑149 (MQ=255)                                          cATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  <  7:43252/149‑1 (MQ=255)                                           aTTCAGGAGAAGTTGGTTGCTGTGAATCTCACCGTGGGTATCCAGGTCGAAATCGTGCCGAACTTCATTGGCCGCCTTGGTGGCAAGGCTCGCCACCATTGCTGGTGGCGCTACCTGCTCAACCTGCCCGAGTGTGACCAGGGCGTAg                                                                                                  <  1:262549/148‑1 (MQ=255)                                            ttCAGGAGAAGTTCGGTGCGGTGAATCACCCCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  >  2:178076/1‑147 (MQ=255)                                            ttCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  <  1:178076/147‑1 (MQ=255)                                               aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  >  4:163132/1‑144 (MQ=255)                                               aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  <  3:163132/144‑1 (MQ=255)                                               aGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACAAGGGCGTAGGCAt                                                                                              >  6:149881/1‑148 (MQ=255)                                                     aGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                                       <  8:151545/149‑1 (MQ=255)                                                       ttCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  <  1:204457/136‑1 (MQ=255)                                                       ttCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCATGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAg                                                                                                  >  2:204457/1‑136 (MQ=255)                                                               ggtgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCa                                                                             <  1:78163/149‑1 (MQ=255)                                                                 tgAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGc                                                                           >  7:35675/1‑149 (MQ=255)                                                                   aaTCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGTGCAGCg                                                                          <  6:194403/148‑1 (MQ=255)                                                                     tCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGc                                                                       <  4:66241/149‑1 (MQ=255)                                                                      cacaCCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAg                                                                        >  8:338993/1‑147 (MQ=255)                                                                          ccGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGAtt                                                                  <  6:12630/149‑1 (MQ=255)                                                                           cGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTg                                                                 >  5:166886/1‑149 (MQ=255)                                                                            gTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGt                                                                >  3:3204/1‑149 (MQ=255)                                                                            gTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGt                                                                <  4:3204/149‑1 (MQ=255)                                                                                 tACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGc                                                           <  1:237366/149‑1 (MQ=255)                                                                                    ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGa                                                                                            <  3:233777/113‑1 (MQ=255)                                                                                    ccAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGa                                                                                            >  4:233777/1‑113 (MQ=255)                                                                                          cGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                  <  3:215362/149‑1 (MQ=255)                                                                                           gACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCATCGAGCAGATTGTCGGGCTGGTGCAggg                                                 <  1:196409/149‑1 (MQ=255)                                                                                            aCATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                                       >  5:40588/1‑94 (MQ=255)                                                                                            aCATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAggg                                                                                                       <  6:40588/94‑1 (MQ=255)                                                                                              aTCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGa                                              >  1:243887/1‑149 (MQ=255)                                                                                                 gTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGAcc                                            <  8:102408/148‑1 (MQ=255)                                                                                                  tGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCg                                           >  6:148653/1‑148 (MQ=255)                                                                                                      gACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                  <  1:252281/137‑1 (MQ=255)                                                                                                      gACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                  >  2:252281/1‑137 (MQ=255)                                                                                                       aCCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATc                                       >  5:368756/1‑147 (MQ=255)                                                                                                       aCCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATc                                       <  6:368756/147‑1 (MQ=255)                                                                                                               gggCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCAGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGt                             >  7:274126/1‑149 (MQ=255)                                                                                                                ggCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGt                             <  8:35675/148‑1 (MQ=255)                                                                                                                 gccgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTc                            >  5:226104/1‑148 (MQ=255)                                                                                                                  ccgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTcc                           >  8:125703/1‑148 (MQ=255)                                                                                                                   cgccTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCaa                         <  5:148653/149‑1 (MQ=255)                                                                                                                            ggCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                                       <  5:210615/78‑1 (MQ=255)                                                                                                                            ggCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAg                                                                                       >  6:210615/1‑78 (MQ=255)                                                                                                                            ggCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAg                        <  3:198524/141‑1 (MQ=255)                                                                                                                            ggCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAg                        >  4:198524/1‑141 (MQ=255)                                                                                                                             gCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACa                    >  1:46803/1‑144 (MQ=255)                                                                                                                             gCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACa                    <  2:46803/144‑1 (MQ=255)                                                                                                                                  gCTCGCGACCTGTGGTGGGGGGGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGGGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCGTGATCGAGTCCAAGGACATGAGCCcgc           <  5:149881/148‑1 (MQ=255)                                                                                                                                   cTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                  >  6:143241/1‑108 (MQ=255)                                                                                                                                   cTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAgg                                                  <  5:143241/108‑1 (MQ=255)                                                                                                                                      gCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATg      >  8:222240/1‑149 (MQ=255)                                                                                                                                       ccaccaGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGAATTATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACAt                   <  5:20932/135‑1 (MQ=255)                                                                                                                                       ccaccaGTGCTGGGGGCCCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGCTGCAGGGTGCCCGTATCATGATCGAGTCCAAGCACAt                   >  6:20932/1‑135 (MQ=255)                                                                                                                                          ccaGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTAcc  >  5:227822/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                           CTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGGGCGTAGGCATGATCAAGCCGTCGGGCAGCGAGCAGATTGTCGGGCTGGTGCAGGGTGACCGTATCCTGATCGAGTCCAAGGACATGAGCCCGCGGATGTACC  >  NC_002947/2062709‑2062994

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