Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,683,461	T→A	5.3%	F225I (TTC→ATC)	PP_3249 →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,683,461	0	T	A	5.3%	230.5 / 2.7	75	F225I (TTC→ATC)	PP_3249	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base A (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

AGTGGCCACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTC  >  NC_002947/3683314‑3683608                                                                                                                                                     |                                                                                                                                                    acgggCCACGACGAAGGCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCtt                                                                                                                                                    <  6:253155/146‑1 (MQ=255) aGTGGCCACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCtt                                                                                                                                                    <  2:9106/149‑1 (MQ=255)   tGGCCACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTc                                                                                                                                                   >  4:331211/1‑148 (MQ=255)     tccACGACGAAGTCGCGCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTcgac                                                                                                                                                <  1:153461/148‑1 (MQ=255)     gCCACGTCGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGTTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTcgac                                                                                                                                                <  5:201472/149‑1 (MQ=255)     gCCACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAAGTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTcga                                                                                                                                                 <  7:35293/148‑1 (MQ=255)     gCCACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTcga                                                                                                                                                 <  7:29795/148‑1 (MQ=255)       cACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTcgac                                                                                                                                                >  2:153461/1‑147 (MQ=255)        acgacgAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGAt                                                                                                                                             <  1:186831/149‑1 (MQ=255)          gacgaAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATgg                                                                                                                                           >  3:290296/1‑149 (MQ=255)             gaAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCATCGTCGAAGGccc                                                                                                                                        >  7:231777/1‑148 (MQ=255)              aaGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCgg                                                                                                                                       <  7:110445/149‑1 (MQ=255)               aGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCATCGTCGAAGGcccga                                                                                                                                      <  3:190885/149‑4 (MQ=255)                 tCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTg                                                                                                                                    >  1:189760/1‑149 (MQ=255)                   gCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGtt                                                                                                                                  >  3:346505/1‑149 (MQ=255)                   gCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCATCGTCGAAGGcccgattg                                                                                                                                   <  8:231777/148‑7 (MQ=255)                     ccGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCg                                                                                                                                >  8:233188/1‑149 (MQ=255)                      cgaaTCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCCCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGc                                                                                                                               <  8:67146/146‑1 (MQ=255)                            gTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCgg                                                                                                                           >  1:61598/1‑147 (MQ=255)                            gTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCgg                                                                                                                           <  2:61598/147‑1 (MQ=255)                             tGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCACCCGGGt                                                                                                                         >  7:277656/1‑148 (MQ=255)                                gcAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCaaa                                                                                                                     <  3:331211/149‑1 (MQ=255)                                         cGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGCGCCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAaggccaggc                                                                                                             >  1:264876/1‑149 (MQ=255)                                               aCTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCACCCGGGTCAAAGGCCAGGCCTTGCCg                                                                                                      <  8:277656/149‑1 (MQ=255)                                                 ttCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGCCGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCggg                                                                                                    >  6:92790/1‑149 (MQ=255)                                                 atCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CAGGGAACTGCTGCCGCTGTGCAAGGACAAGGGCGATGCGGAGTTCATCAATCGGGCGTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCggg                                                                                                    >  6:368621/2‑149 (MQ=255)                                                    aTCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTg                                                                                                 >  6:341233/1‑149 (MQ=255)                                                    aTCCAGATCAACTACTCGCTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTCCCAGGACAAGTGCGGTGCGGTGATGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGCGTCAAAGGCCAGGCCTTGCCGGGCTg                                                                                                 >  5:258345/1‑149 (MQ=255)                                                     tCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTgg                                                                                                <  4:290296/149‑1 (MQ=255)                                                     tCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTgg                                                                                                >  6:62768/1‑149 (MQ=255)                                                     tCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTgg                                                                                                >  5:284836/1‑149 (MQ=255)                                                        aGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTggg                                                                                               <  2:244037/147‑1 (MQ=255)                                                        aGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTggg                                                                                               >  1:244037/1‑147 (MQ=255)                                                          aTCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGgccgc                                                                                           <  2:264876/149‑1 (MQ=255)                                                              actactCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAaggccaggc                                                                                                             >  3:88652/1‑127 (MQ=255)                                                              actactCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAaggccaggc                                                                                                             <  4:88652/127‑1 (MQ=255)                                                                  ctCGGTCATTGCGCCGGACGCTGC‑CCGGGAACTGCTGCGGGTGTGCCAGGACTAGGGCGTTGCGGTGTTGTTCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAg                                                                                   <  6:258345/149‑1 (MQ=255)                                                                     ggTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGc                                                                                 <  7:293273/148‑1 (MQ=255)                                                                       tCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATc                                                                              >  2:86266/1‑149 (MQ=255)                                                                             cgcCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGATGCCTGCTTCGGcgg                                                                         <  5:255688/148‑3 (MQ=255)                                                                              gcCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCTGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCt                                                                       <  7:233188/149‑1 (MQ=255)                                                                              gcCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCt                                                                       >  3:313727/1‑149 (MQ=255)                                                                              gcCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCt                                                                       <  4:346505/149‑1 (MQ=255)                                                                               ccGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGCCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTg                                                                      >  6:82759/1‑149 (MQ=255)                                                                                cGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGGAGCAGGCATCGGCAGCTgg                                                                     <  2:189760/149‑1 (MQ=255)                                                                                cGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTgg                                                                     <  5:284971/149‑1 (MQ=255)                                                                                  gACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCa                                                                          >  7:85130/1‑142 (MQ=255)                                                                                  gACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCa                                                                          <  8:85130/142‑1 (MQ=255)                                                                                         c‑ccGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCt                                                             <  5:341233/148‑1 (MQ=255)                                                                                                aaCTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTtctc                                                      <  6:284836/149‑1 (MQ=255)                                                                                                  ctgctgCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCaa                                                    >  6:34137/1‑149 (MQ=255)                                                                                                   tgctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCgg                                                                            >  8:247923/1‑124 (MQ=255)                                                                                                   tgctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCgg                                                                            <  7:247923/124‑1 (MQ=255)                                                                                                    gctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCaa                                                    <  7:18356/147‑1 (MQ=255)                                                                                                    gctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCaa                                                    >  8:18356/1‑147 (MQ=255)                                                                                                     ctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGt                                                  <  4:313727/148‑1 (MQ=255)                                                                                                     ctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTAGGCACAGCTGTTTCTCAAGtt                                                 >  8:334389/1‑149 (MQ=255)                                                                                                      tgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTc                                                >  6:327088/1‑149 (MQ=255)                                                                                                        ccGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTc                                                <  5:261166/147‑1 (MQ=255)                                                                                                        ccGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGc                                              <  5:62768/149‑1 (MQ=255)                                                                                                              tgCCAGGTCATGGTCGTTGCGGTGTTGATCTATCCGTCCTTCGCCGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCTGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCAATTCTGTTTCTCATGTTGGCCATCa                                         <  4:133943/148‑1 (MQ=255)                                                                                                              tgCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCa                                         <  4:123514/148‑1 (MQ=255)                                                                                                               gCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGc                                       >  1:54398/1‑149 (MQ=255)                                                                                                                ccAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGcc                                      <  5:34137/149‑1 (MQ=255)                                                                                                                ccAGGACAAGGGCGTTGAGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTtct                                                       <  5:234459/132‑1 (MQ=255)                                                                                                                ccAGGACAAGGGCGTTGAGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTtct                                                       >  6:234459/1‑132 (MQ=255)                                                                                                                       aaGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAATGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCgg                                >  3:41976/1‑148 (MQ=255)                                                                                                                       aaGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTgg                                                                     >  8:97650/1‑111 (MQ=255)                                                                                                                       aaGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCAAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCgg                                >  5:105192/1‑148 (MQ=255)                                                                                                                                gCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAATGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCa                       <  4:41976/148‑1 (MQ=255)                                                                                                                                    tgtTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGt                   >  1:311511/1‑148 (MQ=255)                                                                                                                                        gATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATAc              >  4:14464/1‑149 (MQ=255)                                                                                                                                         aTCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAg                                                               >  7:322208/1‑99 (MQ=255)                                                                                                                                         aTCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAg                                                               <  8:322208/99‑1 (MQ=255)                                                                                                                                         aTCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATAcc             <  2:54398/149‑1 (MQ=255)                                                                                                                                               ccggCCATCGTCGAAGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTAc       >  2:248884/3‑149 (MQ=255)                                                                                                                                                    cTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTc  >  4:276615/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGTGGCCACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTC  >  NC_002947/3683314‑3683608

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