Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,990,451	T→C	18.3%	E72E (GAA→GAG)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,990,451	0	T	C	18.3%	163.3 / 25.7	71	E72E (GAA→GAG)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base C (9/4);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGC  >  NC_002947/5990305‑5990596                                                                                                                                                   |                                                                                                                                                  gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTc                                                                                                                                                 >  7:223699/1‑149 (MQ=255) gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTc                                                                                                                                                 >  2:377355/1‑149 (MQ=255)   cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                                                                                                               >  3:245317/1‑149 (MQ=255)   cACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                                                                                                               <  1:72648/149‑1 (MQ=255)             tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCgg                                                                                                                                     >  5:77607/1‑149 (MQ=255)              cGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCggg                                                                                                                                    >  3:386495/1‑149 (MQ=255)                   acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAtt                                                                                                                               >  6:324430/1‑149 (MQ=17)                    caacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                                                                                              >  1:25163/1‑149 (MQ=255)                     aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttcg                                                                                                                             >  3:152974/1‑147 (MQ=17)                     aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTgg                                                                                                                             >  4:383305/1‑149 (MQ=255)                      acacCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                                                                                                                            <  6:220060/149‑1 (MQ=255)                        gccggCAATGATCCGCGGCTTCTTGCCGGTCGCCTTCTCCGCCTTGACGATATCCTCGTTGTCCGAGGCGGTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTGCAGCTCGGGCACCAGCGGGCCTTTcgg                                                                                                                            <  4:295587/143‑1 (MQ=14)                          cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCt                                                                                                                        <  7:76632/149‑1 (MQ=17)                          cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCAGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCt                                                                                                                        >  8:359423/1‑149 (MQ=255)                          cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGAGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCt                                                                                                                        >  4:325053/1‑149 (MQ=255)                              aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                    >  8:320704/1‑149 (MQ=14)                                  aTCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                                                                                                <  4:257665/149‑1 (MQ=255)                                     aGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                             >  7:134456/1‑149 (MQ=14)                                     aGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                             <  3:171109/149‑1 (MQ=14)                                        tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCAGCCGGGAACATATCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggt                                                                                                          >  7:77892/1‑149 (MQ=14)                                           ttcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTggtggt                                                                                                       >  2:240079/1‑149 (MQ=11)                                           atcttGCCGGTCGCCTTCACCCCCTTCACGAATTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGCCAATGTACGGCGCATCCGGGAACAGGCCTTTCAGTTCCGGCACCAGCCGGCCGTTGGTGCCTAGCTCGAAGCTggtggt                                                                                                       >  4:260478/2‑149 (MQ=17)                                            tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTggaggt                                                                                                       >  5:30779/1‑148 (MQ=255)                                                gCCGGTCGCCTTCCCCGCCTTGACGAAATCCTCGTTGGCCCAGGCGTTGATCTGGCCAGGGCTGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                     <  2:153216/146‑1 (MQ=11)                                                gCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                     <  2:334273/146‑1 (MQ=255)                                                gCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGGACGGCGCCTCCGGGAACATCTCTTTCAGCTCGGGCACCCGCGGGCCATTCTGGCCTTGCTCGAAGCTGGGGGTGAGGa                                                                                                  >  8:368906/1‑149 (MQ=11)                                                 ccGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                                 <  4:386495/149‑1 (MQ=255)                                                   ggTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATgg                                                                                               <  1:377355/149‑1 (MQ=255)                                                    gTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                              >  1:316964/1‑149 (MQ=11)                                                    gTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGGTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                              <  6:30779/149‑1 (MQ=255)                                                       gCCTTCACCGCCTTGACGAAATCCTCGTGGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCgg                                                                                           <  3:259737/149‑1 (MQ=255)                                                         cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACACGCCTAACAGTTCAGGCACCAGCGGGCCGATGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCgg                                                                                           >  1:322985/1‑147 (MQ=17)                                                            cACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGt                                                                                      >  6:335820/1‑149 (MQ=255)                                                                gCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaa                                                                                  <  5:335820/149‑1 (MQ=255)                                                                 ccTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTg                                                                                    >  7:364214/1‑146 (MQ=255)                                                                     gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAAc                                                                              <  2:25163/148‑1 (MQ=255)                                                                       cGATATCCTCGTTGGCGCTGCCGTTCATCTGGCCAGGGCTGGCACTGTACGGCGCATCCGTGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGGCTTGCTGGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGATGAACTTg                                                                           <  4:249704/149‑1 (MQ=18)                                                                        gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTgg                                                                          >  1:333474/1‑149 (MQ=255)                                                                        gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTgg                                                                          >  1:381087/1‑149 (MQ=255)                                                                        gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCACGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTgg                                                                          >  3:175994/1‑149 (MQ=255)                                                                              cTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTTGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAgg                                                                     <  4:245317/148‑1 (MQ=255)                                                                                cGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGt                                                                    >  5:313100/1‑147 (MQ=255)                                                                                 gTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTca                                                                  <  6:313100/148‑1 (MQ=255)                                                                                      cccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGCCCAGGTcagccagc                                                            >  5:31719/1‑149 (MQ=255)                                                                                       ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt                                                                            >  5:111913/1‑132 (MQ=255)                                                                                       ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt                                                                            <  6:111913/132‑1 (MQ=255)                                                                                       ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcg                                                           >  7:88732/1‑149 (MQ=255)                                                                                             gTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCa                                                     >  8:182845/1‑149 (MQ=255)                                                                                                gATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgt                                                  >  4:141164/1‑149 (MQ=255)                                                                                                 aTCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCt                                                                                                                    >  5:349300/1‑82 (MQ=255)                                                                                                 aTCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCt                                                                                                                    <  6:349300/82‑1 (MQ=255)                                                                                                   cTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttg                                                <  7:359423/148‑1 (MQ=255)                                                                                                      gCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttgtt                                              <  3:383305/147‑1 (MQ=255)                                                                                                      gCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCt                                            >  5:342720/1‑149 (MQ=255)                                                                                                        cAGGGCGGGCAATGTACGGCGCATCCGGGAAGCGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTGGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAAGTTGGCGAGGTCAGCCAGCGCCAGCACGTTGTTCtt                                           >  4:122731/1‑148 (MQ=17)                                                                                                         aGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGa                                         <  7:335374/149‑1 (MQ=255)                                                                                                             cgggcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACtt                                     <  6:342720/149‑1 (MQ=255)                                                                                                               ggcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCAAGAAcaacg                                   <  4:175994/149‑5 (MQ=255)                                                                                                                gcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTc                                  >  1:204342/1‑149 (MQ=255)                                                                                                                 cAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                                    >  1:275495/1‑82 (MQ=255)                                                                                                                 cAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                                    <  2:275495/82‑1 (MQ=255)                                                                                                                 cAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCt                                 <  2:381087/149‑1 (MQ=255)                                                                                                                  aaTGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTg                                <  3:260478/149‑1 (MQ=255)                                                                                                                   aTGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTgg                               >  3:117240/1‑149 (MQ=255)                                                                                                                            gcATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATc                      <  2:333474/149‑1 (MQ=255)                                                                                                                            gcATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATc                      >  1:111085/1‑149 (MQ=255)                                                                                                                                  gggAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCAc                >  1:197760/1‑149 (MQ=255)                                                                                                                                         gggCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGAc         <  3:325053/149‑1 (MQ=255)                                                                                                                                               tCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttg                                                <  2:248571/104‑1 (MQ=255)                                                                                                                                                  gTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcgcca                                                        <  6:331208/93‑1 (MQ=255)                                                                                                                                                  gTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGc  >  2:91802/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGC  >  NC_002947/5990305‑5990596

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