Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,257,872	G→A	5.9%	intergenic (‑71/+73)	PP_4634 ← / ← PP_4635	alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,257,872	0	G	A	5.9%	252.0 / 3.3	85	intergenic (‑71/+73)	PP_4634/PP_4635	alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base A (2/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

CAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCC  >  NC_002947/5257725‑5258019                                                                                                                                                    |                                                                                                                                                          cagcagAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGa                                                                                                                                                          <  7:333243/149‑1 (MQ=255)  agcagAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGaa                                                                                                                                                         >  5:173182/1‑149 (MQ=255)      gAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCg                                                                                                                                                     >  2:343389/1‑149 (MQ=255)        aGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCt                                                                                                                                                   <  6:272343/149‑1 (MQ=255)         gccgccGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGATAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTc                                                                                                                                                  >  8:206062/1‑149 (MQ=255)           cgccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGCAAGCCCGCTc‑‑cc                                                                                                                                              <  8:115045/149‑2 (MQ=255)            gccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑t                                                                                                                                            >  2:191609/1‑148 (MQ=255)             ccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TAc                                                                                                                                          <  2:238525/149‑1 (MQ=255)              cgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTc‑‑cca                                                                                                                                             >  3:348415/1‑145 (MQ=255)                aGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAgg                                                                                                                                       <  2:314870/149‑1 (MQ=255)                  gCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAgggg                                                                                                                                     <  5:22485/149‑1 (MQ=255)                  gCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAgggg                                                                                                                                     <  3:200699/149‑1 (MQ=255)                   cTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGATAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAggggg                                                                                                                                    <  7:206062/149‑1 (MQ=255)                    tGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑Ta                                                                                                                                           >  1:239693/1‑141 (MQ=255)                    tGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑Ta                                                                                                                                           <  2:239693/141‑1 (MQ=255)                      tATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGgag                                                                                                                                  >  4:148543/1‑148 (MQ=255)                      tATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGga                                                                                                                                   <  7:259373/147‑1 (MQ=255)                        tCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAAc                                                                                                                               <  4:406477/149‑1 (MQ=255)                         ccGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGACTCTTCGCGGGTGAACCCGCTC‑‑C‑‑AACAGGGGGAGAAcc                                                                                                                              >  8:51487/1‑149 (MQ=255)                          cGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑t                                                                                                                           >  5:83150/1‑149 (MQ=255)                          cGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCAACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑t                                                                                                                           >  1:306526/1‑149 (MQ=255)                               gCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTgg                                                                                                                       >  2:153618/1‑148 (MQ=255)                                  gcgGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGc                                                                                                                   <  8:399010/149‑1 (MQ=255)                                    ggggTGGCGAAATGGGCTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGGGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCgg                                                                                                                 <  1:287036/149‑1 (MQ=255)                                     gggTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCgg                                                                                                                 <  1:374940/148‑1 (MQ=255)                                       gTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCt                                                                                                              <  6:306163/149‑1 (MQ=255)                                          gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCCAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTga                                                                                                           >  1:126410/1‑148 (MQ=255)                                          gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTg                                                                                                            <  4:154131/148‑1 (MQ=255)                                          gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGc                                                                                                           >  6:229003/1‑149 (MQ=255)                                             aaaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCcg                                                                                                        >  2:224413/1‑149 (MQ=255)                                              aaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCtt                                                                                                             >  5:26867/1‑143 (MQ=255)                                              aaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCtt                                                                                                             <  6:26867/143‑1 (MQ=255)                                                    gTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAGCTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                 >  6:326808/1‑149 (MQ=255)                                                     tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTC‑‑C‑‑CACAGGTTCTCCCCC‑‑TGTAGGAGCGGGTTCACCCGCGAagag                                                                                                 <  3:303817/148‑1 (MQ=14)                                                     tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTC‑‑C‑‑CACAGGTTCTCCCCC‑‑TGTAGGAGCGGGTTCACCCGCGAagag                                                                                                 >  4:303817/1‑148 (MQ=14)                                                       gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAg                                                                                                                    <  5:273805/127‑1 (MQ=255)                                                       gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAg                                                                                                                    >  6:273805/1‑127 (MQ=255)                                                         gggTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTC‑‑C‑‑CACAGGTTCTCCCCC‑‑TGTAGGAGCGGGTTCACCCGCGAAGAGGccag                                                                                            <  4:348415/149‑3 (MQ=14)                                                                aTGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCTCCAAAACTAGCAAAGATTTTGTAGCGTGTACTGGCCTGTGCGCGGGTGAACCCGCTC‑‑C‑‑TACTGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGc                                                                                     <  4:407597/149‑1 (MQ=39)                                                                        aCTCAAGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAGGAGGCCGGTACAGGCTTACCCGa                                                                             <  3:95461/149‑1 (MQ=255)                                                                         cTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTATCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTC‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTAccc                                                                               <  1:343239/148‑1 (MQ=255)                                                                          tCAAGGCCCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGCCCTCATCGCGGGTGAACCCCCTC‑‑C‑‑TACACGGGGAGAACC‑‑TGTGGGAGCGCGCTTGCCCGCGAAGAGGCAGGTACAGGCTTACCCGAt                                                                            >  6:322057/1‑148 (MQ=255)                                                                            aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                          >  1:203313/1‑148 (MQ=255)                                                                                cgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCt                                                                     >  5:375399/1‑149 (MQ=255)                                                                                cgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCt                                                                     <  5:229003/149‑1 (MQ=255)                                                                                 gcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTc                                                                    >  4:87436/1‑149 (MQ=255)                                                                                     gTCCTGTTGTTGTGCGCCAAAAGTCGCAAAGATTTTGTAGCCTTTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGAc                                                                 <  7:51487/148‑1 (MQ=255)                                                                                         tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGt                                                            >  8:183983/1‑149 (MQ=255)                                                                                         tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑A‑‑TACAGGGGGAGAACC‑‑TGTGTGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACACGATTACAGCTCGACGCAg                                                             >  7:279139/1‑148 (MQ=255)                                                                                         tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTAATGGCCTCTTACCGGGTGCAACCCCGC‑‑C‑‑TACAGGGGGAGAAAA‑‑TGTGGGATCGGCCATGCCACACAAGAGCCCGGAACATGCATACCCGATTACATCTATACGCAGt                                                            >  8:215018/1‑149 (MQ=255)                                                                                           ttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTc                                                           <  6:375399/148‑1 (MQ=255)                                                                                              ttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGaa                                                        <  3:253494/148‑1 (MQ=255)                                                                                                gtgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACg                                                                <  5:211068/138‑1 (MQ=255)                                                                                                gtgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACg                                                                >  6:211068/1‑138 (MQ=255)                                                                                                 tgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACtt                                                     >  2:144340/1‑148 (MQ=255)                                                                                                  gcgcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑CGGTACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACtt                                                     >  1:125529/1‑149 (MQ=255)                                                                                                   cgcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTca                                                   <  7:183983/148‑1 (MQ=255)                                                                                                    gcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTcac                                                  >  6:99141/1‑148 (MQ=255)                                                                                                     ccAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACAt                                                >  4:363119/1‑149 (MQ=255)                                                                                                       aaaaCTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑AACAGGGGGAGCACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGGACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCg                                              >  7:141517/1‑149 (MQ=255)                                                                                                          ccttgCAATGATTTTGTAGCCTTTACTGGCCGATTCGCGGGTGAACCCGCTC‑‑C‑‑TACCGGGGGATAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCTATTACATCTCGACGCAGTCGAACTTCACATCGGTg                                           <  7:215018/145‑1 (MQ=38)                                                                                                            tAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGc                                         >  7:397636/1‑149 (MQ=255)                                                                                                               cAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCAc                                      <  1:224413/149‑1 (MQ=255)                                                                                                                aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACCTGTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCa                                       >  1:126038/1‑149 (MQ=255)                                                                                                                aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                     >  1:415501/1‑149 (MQ=255)                                                                                                                aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                     >  7:159754/1‑149 (MQ=255)                                                                                                                aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACCCAGTCGAACATCACATCGGTGGCCAcc                                     >  4:224884/1‑148 (MQ=255)                                                                                                                 aaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGCACTTCACATCGGTGGCCACGt                                    >  1:205446/1‑149 (MQ=255)                                                                                                                   gATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCt                                  <  5:312260/149‑1 (MQ=255)                                                                                                                   gATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCt                                  <  3:148543/149‑1 (MQ=255)                                                                                                                    aTTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCtt                                 >  3:337220/1‑149 (MQ=255)                                                                                                                     ttttGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                >  4:138215/1‑149 (MQ=255)                                                                                                                     ttttGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                <  1:191609/149‑1 (MQ=255)                                                                                                                         gTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTc                             <  6:83150/148‑1 (MQ=255)                                                                                                                            gCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTa                          <  3:138215/148‑1 (MQ=255)                                                                                                                            gCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAg                         >  4:270990/1‑149 (MQ=255)                                                                                                                            gCCTCTACTGTCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCAGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTa                          <  8:279139/148‑1 (MQ=255)                                                                                                                                 tACTGGCCTCTTCGGGGGTGATCCCGGTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCTGGAAGAGGCCGGTAGAGGCTTACCCGATTAGAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCaa                     <  5:322057/148‑1 (MQ=255)                                                                                                                                 tACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                    <  1:153618/149‑1 (MQ=255)                                                                                                                                   cTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACg                   <  3:87436/148‑1 (MQ=255)                                                                                                                                        ctctTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg             >  1:227049/1‑149 (MQ=255)                                                                                                                                         tctTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgc            >  6:311440/1‑149 (MQ=255)                                                                                                                                           ttCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcgc          >  3:186527/1‑149 (MQ=255)                                                                                                                                             cgcgGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcgc          <  7:32180/147‑1 (MQ=255)                                                                                                                                                   tGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc  <  2:227049/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                          CAGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCC  >  NC_002947/5257725‑5258019

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