Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,567,688	A→C	8.1%	E82A (GAG→GCG)	PP_3150 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,567,688	0	A	C	8.1%	187.4 / 7.2	62	E82A (GAG→GCG)	PP_3150	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base C (2/3);  total (29/33)
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.32e-01

CTGCCTTTGCCCAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACGGCATTGCACTGCT  >  NC_002947/3567540‑3567835                                                                                                                                                     |                                                                                                                                                    cTGCCTTTGCCCAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGa                                                                                                                                                     <  1:17415/149‑1 (MQ=255)        tGCCCAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCTCCTGCGcgtcg                                                                                                                                               <  6:101094/148‑1 (MQ=255)           ccAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTc                                                                                                                                           >  6:52083/1‑149 (MQ=255)            cAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCTCCTGCGcgtcgg                                                                                                                                              >  4:71419/1‑144 (MQ=255)            cAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCTCCTGCGcgtcgg                                                                                                                                              <  3:71419/145‑2 (MQ=255)              gCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGcc                                                                                                                                        <  6:39973/149‑1 (MQ=255)                    aGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGc                                                                                                                                         >  6:242260/1‑142 (MQ=255)                    aGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGc                                                                                                                                         <  5:242260/142‑1 (MQ=255)                     ggCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCg                                                                                                                                       <  6:263373/143‑1 (MQ=255)                     ggCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCg                                                                                                                                       >  5:263373/1‑143 (MQ=255)                       ccGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCTTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTg                                                                                                                               <  1:298298/149‑1 (MQ=255)                       ccGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTg                                                                                                                               <  5:150395/149‑1 (MQ=255)                        cGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTgc                                                                                                                              >  8:181722/1‑149 (MQ=255)                              ccTCGCCAGCGCGGCCGTTACCGAGTTCTTCCCCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTg                                                                                                                               <  7:99276/142‑1 (MQ=255)                              ccTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTg                                                                                                                               >  8:99276/1‑142 (MQ=255)                                tCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGc                                                                                                                       >  6:132947/1‑148 (MQ=255)                                tCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGc                                                                                                                       >  6:75196/1‑148 (MQ=255)                                   ccAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCAAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAAc                                                                                                                   >  3:264846/1‑149 (MQ=255)                                    cAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCACCTCGGCCGACGCTCCTGCGcgtcgg                                                                                                                                              >  1:79209/1‑120 (MQ=255)                                    cAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCACCTCGGCCGACGCTCCTGCGcgtcgg                                                                                                                                              <  2:79209/121‑2 (MQ=255)                                            ccGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGaa                                                                                                          <  5:132947/149‑1 (MQ=255)                                              tttaactaGTTCTACACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTAGTGCAAAAACTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCTGAAATCGTGCGCATCGCCAACGAAGTCGAAg                                                                                                         <  8:194242/141‑1 (MQ=255)                                                  ccGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGc                                                                                                                       <  4:56330/130‑1 (MQ=255)                                                  ccGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGgtgt                                                                                                     <  1:223920/148‑1 (MQ=255)                                                   cGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAg                                                                                                         <  6:158896/143‑1 (MQ=255)                                                   cGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAg                                                                                                         >  5:158896/1‑143 (MQ=255)                                                   cGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCg                                                                                                   >  2:137821/1‑149 (MQ=255)                                                    gAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGc                                                                                                  <  5:52083/149‑1 (MQ=255)                                                    gAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGc                                                                                                  >  2:10402/1‑149 (MQ=255)                                                            caccacCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGgc                                                                                          >  3:33212/1‑149 (MQ=255)                                                             accaccGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGgc                                                                                          <  7:181722/148‑1 (MQ=255)                                                               caccGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCt                                                                                       <  4:33212/149‑1 (MQ=255)                                                                        tgttgtCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGGCAACGAAGTCGAAGGTGTCGccccc                                                                                              >  5:264189/1‑133 (MQ=255)                                                                        tgttgtCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGccccc                                                                                              <  6:264189/133‑1 (MQ=255)                                                                          ttgtCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCa                                                                             <  3:143995/148‑1 (MQ=255)                                                                            gtCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAgg                                                                           <  3:173404/148‑1 (MQ=255)                                                                                aTGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGca                                                                      <  5:75196/149‑1 (MQ=255)                                                                                aTGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGca                                                                      >  8:70787/1‑149 (MQ=255)                                                                                         ccGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAAtt                                                              <  3:89825/148‑1 (MQ=255)                                                                                                    gTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGc                                                  >  6:238426/1‑149 (MQ=255)                                                                                                       cAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGtt                                               >  5:59341/1‑149 (MQ=255)                                                                                                        aaaaCCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCTTTTCGCTGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTc                                              >  6:291757/1‑149 (MQ=255)                                                                                                          aaCCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGc                                            <  7:28600/149‑1 (MQ=255)                                                                                                             cTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCAt                                         <  6:4304/149‑1 (MQ=255)                                                                                                                ccACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATggg                                      >  5:148352/1‑149 (MQ=255)                                                                                                                ccACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATgg                                       <  7:70787/148‑1 (MQ=255)                                                                                                                ccACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATg                                        <  5:238426/147‑1 (MQ=255)                                                                                                                   cGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCt                                   <  6:148352/149‑1 (MQ=255)                                                                                                                     ccAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTgg                                 <  1:10402/149‑1 (MQ=255)                                                                                                                      cAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTggg                                >  1:209079/1‑149 (MQ=255)                                                                                                                       aGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTggg                                >  2:72511/1‑148 (MQ=255)                                                                                                                         ccTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGc                               <  4:264846/147‑1 (MQ=255)                                                                                                                                   ccTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGa                    >  8:149566/1‑148 (MQ=255)                                                                                                                                       ggCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGaa                                                    >  2:14925/1‑112 (MQ=255)                                                                                                                                       ggCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGaa                                                    <  1:14925/112‑1 (MQ=255)                                                                                                                                       ggCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACg               <  1:137821/149‑1 (MQ=255)                                                                                                                                        gCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACgg              >  2:191761/1‑149 (MQ=255)                                                                                                                                                gcAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATgg                                       <  7:109455/116‑1 (MQ=255)                                                                                                                                                gcAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATgg                                       >  8:109455/1‑116 (MQ=255)                                                                                                                                                  aGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACGGCATTGCActg    <  7:149566/149‑1 (MQ=255)                                                                                                                                                  aGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGACCAGGTAGCACAGCAAATCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACGGCATTGCActg    >  1:94551/1‑149 (MQ=255)                                                                                                                                                    gAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACGGCATTGCActgct  >  4:205655/1‑149 (MQ=255)                                                                                                                                                    gAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACGGCATTGCActgct  >  2:289712/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGCCTTTGCCCAGCGACCAGGCCGACCGCCTCGCCAGCGCGGCCGTTACCGAGTTCTTCACCACCGCCAGTGTTGTCGATGCCGCTACCGTTCTGTTCGTGCAAAACCTGCCACGCCAGCCTGCCACCGCCTCGGCCGACGCGCAGGAGCGTCGTTTCGCCGAAATCGTGCGCATCGCCAACGAAGTCGAAGGTGTCGCCCCCGGGCGCTTCCAGGCCCAGGTAGCACAGCAATTCGATCGCGAAGCGTTCGCCATGGGCCTGGGCATGCTGGGCGAAAACGGCATTGCACTGCT  >  NC_002947/3567540‑3567835

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