Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,437,893	C→A	5.7%	T45K (ACG→AAG)	PP_3057 →	pectinesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,437,893	0	C	A	5.7%	227.8 / 3.7	70	T45K (ACG→AAG)	PP_3057	pectinesterase
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base A (2/2);  total (37/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 = 5.00e-01

ATAAGACCCTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATACGCTCGCC  >  NC_002947/3437750‑3438035                                                                                                                                                 |                                                                                                                                               aTAAGACCCTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTc                                                                                                                                            >  5:113021/1‑148 (MQ=255)      aCCCTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTc                                                                                                                                      <  6:113021/149‑1 (MQ=255)       cccTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTc                                                                                                                                            <  5:314679/142‑1 (MQ=255)       cccTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTc                                                                                                                                            >  6:314679/1‑142 (MQ=255)        ccTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCaa                                                                                                                                    >  5:205125/1‑149 (MQ=255)          tATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTc                                                                                                                                            <  5:761/139‑1 (MQ=255)          tATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTc                                                                                                                                            >  6:761/1‑139 (MQ=255)            tGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAAcc                                                                                                                                  <  1:4226/147‑1 (MQ=255)            tGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAg                                                                                                                                <  6:125315/149‑1 (MQ=255)             gAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGt                                                                                                                               >  4:61365/1‑149 (MQ=255)             gAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGt                                                                                                                               >  1:289821/1‑149 (MQ=255)              aGTCTCGAAACTGAACTTGCAGGCCGCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGt                                                                                                                               <  4:279536/148‑1 (MQ=255)              aGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGtt                                                                                                                              >  4:260367/1‑149 (MQ=255)                tctcGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAAAGTTCTACGTc                                                                                                                                      >  4:170642/1‑139 (MQ=255)                tctcGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAAAGTTCTACGTc                                                                                                                                      >  4:175130/1‑139 (MQ=255)                tctcGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAAAGTTCTACGTc                                                                                                                                      <  3:170642/139‑1 (MQ=255)                tctcGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAAAGTTCTACGTc                                                                                                                                      <  3:175130/139‑1 (MQ=255)                  tcGAAACTGAAATTGCAGGCCTCTCCAGCAGAGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGAt                                                                                                                           >  1:203161/1‑149 (MQ=255)                     aaaCTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGc                                                                                                                       <  5:164236/149‑1 (MQ=255)                      aaCTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGa                                                                                                                            <  2:305211/143‑1 (MQ=255)                      aaCTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGa                                                                                                                            >  1:305211/1‑143 (MQ=255)                       agtGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCAACAGGTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCGGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGa                                                                                                                     <  7:64141/147‑1 (MQ=255)                       aCTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGa                                                                                                                     <  3:61365/149‑1 (MQ=255)                       aCTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGa                                                                                                                     >  6:66603/1‑149 (MQ=255)                       aCTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGa                                                                                                                     <  2:172380/149‑1 (MQ=255)                            aaTTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACa                                                                                                                <  8:169134/149‑1 (MQ=255)                                  aGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATgcgc                                                                                                          >  1:19371/1‑149 (MQ=255)                                    gCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCtt                                                                                                        <  7:138245/149‑1 (MQ=255)                                           cAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAAtt                                                                                                 >  6:147741/1‑149 (MQ=255)                                            aGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTc                                                                                                <  5:66603/149‑1 (MQ=255)                                             gCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCg                                                                                               <  5:207257/149‑1 (MQ=255)                                              cA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCgg                                                                                              >  5:323997/1‑149 (MQ=255)                                                    tCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGc                                                                                         >  2:47795/1‑149 (MQ=255)                                                       ccTCGCTGATCGACTACGTCACCATCTTCAAAACTTCTAGTGCTAAAGCAATCTCCGAGGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCa                                                                                      <  5:303856/149‑1 (MQ=255)                                                       ccTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCa                                                                                      <  6:19152/149‑1 (MQ=255)                                                           gCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTc                                                                                   <  3:260367/148‑1 (MQ=255)                                                              gATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGccc                                                                                       >  8:76392/1‑141 (MQ=255)                                                              gATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGccc                                                                                       <  7:76392/141‑1 (MQ=255)                                                              gATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCaaaa                                                                               >  4:155996/1‑149 (MQ=255)                                                                  gACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCa                                                                            <  2:19371/148‑1 (MQ=255)                                                                   aCTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGc                                                                                         <  8:26529/134‑1 (MQ=255)                                                                   aCTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGc                                                                                         >  7:26529/1‑134 (MQ=255)                                                                          cACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAg                                                                   >  2:279034/1‑149 (MQ=255)                                                                             cAGCTTCAAACCTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCtt                                                                                                        <  7:176555/109‑1 (MQ=255)                                                                             cAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCtt                                                                                                        >  8:176555/1‑109 (MQ=255)                                                                             cAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCa                                                                >  6:49305/1‑149 (MQ=255)                                                                                 ttCAAAACTTCTAGTGCTAAAGCAATCGACGATGCGGTGGCCGCTTCGCAGCCCATCTCGCGAACGTTCTACGTCAAACAGTTCATCGGCGATGACAATGAGCTAGGTAATTCGGATAGCCCACTCAAAACCATTGCACGAGCCATcgc                                                            >  5:206523/1‑149 (MQ=255)                                                                                         ttCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGc                                                    >  5:12825/1‑149 (MQ=255)                                                                                               tGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATg                                              >  8:197662/1‑149 (MQ=255)                                                                                               tGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATg                                              <  2:203161/149‑1 (MQ=255)                                                                                                 cTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGc                                            >  2:291551/1‑149 (MQ=255)                                                                                                    aaGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCggg                                         >  3:22493/1‑149 (MQ=255)                                                                                                    aaGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCg                                           <  5:49305/147‑1 (MQ=255)                                                                                                    aaGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCg                                           <  6:205125/147‑1 (MQ=255)                                                                                                        aaTCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAg                                     >  8:107925/1‑149 (MQ=255)                                                                                                        aaTCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAg                                     >  1:153891/1‑149 (MQ=255)                                                                                                          tCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCg                                   <  1:291551/149‑1 (MQ=255)                                                                                                            gCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGa                                 >  5:45221/1‑149 (MQ=255)                                                                                                                   cGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCAtt                          <  7:107925/149‑1 (MQ=255)                                                                                                                            cTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGaa                 <  5:147741/149‑1 (MQ=255)                                                                                                                               cgcCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGa               <  7:79220/148‑1 (MQ=255)                                                                                                                                gcCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGAc              >  1:135923/1‑148 (MQ=255)                                                                                                                                 ccGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACta            >  5:301653/1‑149 (MQ=255)                                                                                                                                    gCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATAc         <  2:153891/149‑1 (MQ=255)                                                                                                                                       aTCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATACGCt      >  5:137507/1‑149 (MQ=255)                                                                                                                                          tcGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATACGCTCGc   >  4:246141/1‑149 (MQ=255)                                                                                                                                            gcgAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATACGCTCGcc  >  3:192984/1‑148 (MQ=255)                                                                                                                                            gcgAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATACGCTCGcc  >  4:229101/1‑148 (MQ=255)                                                                                                                                             cgAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTgg                      >  6:148977/1‑127 (MQ=255)                                                                                                                                             cgAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTgg                      <  5:148977/127‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               ATAAGACCCTATGAGTCTCGAAACTGAAATTGCAGGCCTCTCCAGCA‑AGGTCACCTCGCTGATCGACTACGTCACCAGCTTCAAAACTTCTAGTGCTAAAGCAATCGCCGATGCGGTGGCCGCTGCGCCGGCCATCTCGCGAACGTTCTACGTCAACCAGTTGATCGGCGATGACAATGCGCTTGGTAATTCGGATAGCCCACTCAAAACCATTGCCCGAGCCATCGCCGCTACGCCGAATGGCGGGTTAGCGGACATCATTCTGGTCGAAGACTATACGCTCGCC  >  NC_002947/3437750‑3438035

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