Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,683,469	T→A	5.4%	D227E (GAT→GAA)	PP_3249 →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,683,469	0	T	A	5.4%	230.5 / 2.6	74	D227E (GAT→GAA)	PP_3249	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base A (2/2);  total (36/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.82e-01

ACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGCCTGAC  >  NC_002947/3683321‑3683616                                                                                                                                                      |                                                                                                                                                    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                                                                                                        >  6:62768/1‑149 (MQ=255)                                              tCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTgg                                                                                                        <  4:290296/149‑1 (MQ=255)                                              tCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTgg                                                                                                        >  5:284836/1‑149 (MQ=255)                                                 aGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTggg                                                                                                       >  1:244037/1‑147 (MQ=255)                                                 aGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTggg                                                                                                       <  2:244037/147‑1 (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                                                                                  <  8:85130/142‑1 (MQ=255)                                                                           gACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCa                                                                                  >  7:85130/1‑142 (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                                                                                    <  7:247923/124‑1 (MQ=255)                                                                                            tgctgcCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCgg                                                                                    >  8:247923/1‑124 (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                                                               >  6:234459/1‑132 (MQ=255)                                                                                                         ccAGGACAAGGGCGTTGAGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTtct                                                               <  5:234459/132‑1 (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)                                                                                                                                                 gacgaTGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGc       <  5:82759/148‑1 (MQ=255)                                                                                                                                                 gacgaTGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGc       <  5:368621/148‑1 (MQ=255)                                                                                                                                                  acgaTGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGCCt     <  7:267891/149‑1 (MQ=255)                                                                                                                                                  acgaTGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTAGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGCCt     <  7:334389/149‑1 (MQ=255)                                                                                                                                                    gaTGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGCCTGa   >  2:291751/1‑149 (MQ=255)                                                                                                                                                     aTGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGCCTGAc  <  6:105192/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACGACGAAGTCGCCCGTATCGTGCGCAGCGAGCCGCTGGACTTCATCCAGATCAACTACTCGGTCAATGCGCCGGACGCTGC‑CCGGGAACTGCTGCCGCTGTGCCAGGACAAGGGCGTTGCGGTGTTGATCAATCGGGCCTTCGACGATGGCCGGCTGTTCGCCCGGGTCAAAGGCCAGGCCTTGCCGGGCTGGGCCGCCGAAGCAGGCATCGGCAGCTGGGCACAGCTGTTTCTCAAGTTCGCCATCAGCCACCCGGCAGTGACCACGGTCATACCCGCTACCAGTCGGCCTGAC  >  NC_002947/3683321‑3683616

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