Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,647,927	G→A	Y265Y (TAC→TAT)	oprB‑II ←	carbohydrate‑selective porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,647,927	0	G	A	94.9%	173.9 / ‑1.4	59	Y265Y (TAC→TAT)	oprB‑II	carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (29/27);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCG  >  NC_002947/1647783‑1648068                                                                                                                                                  |                                                                                                                                               gTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAATCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaatagt                                                                                                                                            <  4:226417/149‑1 (MQ=255) gTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaata                                                                                                                                              <  2:216043/147‑1 (MQ=255)  tCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaata                                                                                                                                              >  3:56837/1‑146 (MQ=255)  tCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaata                                                                                                                                              <  4:56837/146‑1 (MQ=255)    ggACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCTGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGc                                                                                                                                         <  1:313748/149‑1 (MQ=255)    ggACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGTACATCTTGTGCCTTGGCATTACTGTAATAGTAGc                                                                                                                                         <  1:254583/149‑1 (MQ=255)     gACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGATACTGCTGCGGTAGGCGGCGCCACTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGc                                                                                                                                         >  1:300273/1‑149 (MQ=255)     gACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGcc                                                                                                                                        <  3:160075/149‑1 (MQ=255)      aCGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCg                                                                                                                                       >  2:296558/1‑149 (MQ=255)         ccAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCAGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCa                                                                                                                                    <  2:69169/149‑1 (MQ=255)          cAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCAc                                                                                                                                   <  8:198165/149‑1 (MQ=255)                     ccTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCt                                                                                                                        <  4:19328/149‑1 (MQ=255)                      ctgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCtt                                                                                                                       >  6:359256/1‑149 (MQ=255)                            ctgGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCACTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑Caa                                                                                                                  <  2:300273/147‑1 (MQ=255)                                gcgcCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTg                                                                                                            <  5:210009/149‑1 (MQ=255)                                  gcCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGa                                                                                                          >  7:5505/1‑149 (MQ=255)                                       tCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGt                                                                                                     <  7:278461/149‑1 (MQ=255)                                        ccACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTg                                                                                                    >  6:334680/1‑149 (MQ=255)                                          acaAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTgg                                                                                                           <  5:83281/140‑1 (MQ=255)                                          acaAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTgg                                                                                                           >  6:83281/1‑140 (MQ=255)                                          acaAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTggg                                                                                                  >  2:370006/1‑149 (MQ=255)                                             aGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTc                                                                                               <  7:186336/149‑1 (MQ=255)                                             aGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTc                                                                                               <  1:296558/149‑1 (MQ=255)                                                  tgcttgctCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATAc                                                                                          <  3:262133/149‑1 (MQ=255)                                                  tgcctgctCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATAc                                                                                          <  5:359256/149‑1 (MQ=255)                                                        ctCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTc                                                                                    >  7:32117/1‑149 (MQ=255)                                                         tCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTc                                                                                    >  8:72393/1‑148 (MQ=255)                                                                tgctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACgg                                                                              >  3:71847/1‑148 (MQ=255)                                                                tgctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACgg                                                                              >  3:181401/1‑148 (MQ=255)                                                                    gcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTa                                                                         >  8:191630/1‑149 (MQ=255)                                                                    gcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTa                                                                         >  7:167438/1‑149 (MQ=255)                                                                          ggcggcGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCAc                                                                   >  1:360647/1‑149 (MQ=255)                                                                            cgtcgccgcTGAGGCCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCAccc                                                                 <  1:245033/146‑1 (MQ=255)                                                                                     tGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CACGCCTTGGATACGTGGGGTCCATACCAGTTCGCACGGGATTCCCGCGCCCTGGGtgccg                                                        >  8:147956/1‑149 (MQ=255)                                                                                     tGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATATGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGtgccg                                                        >  4:217465/1‑149 (MQ=255)                                                                                         ggcggcCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGtgccgctgc                                                    >  4:185426/1‑149 (MQ=255)                                                                                         ggcggcCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGtgccgctgc                                                    <  4:342273/149‑1 (MQ=255)                                                                                           cggccggcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGACTTGGATACGTGGGGTGCATACCAGTTCGAACGGCATTACCGCACCCTGGGtgccgctgcgg                                                  >  4:270913/1‑147 (MQ=255)                                                                                                 gcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGc                                            <  5:334680/149‑1 (MQ=255)                                                                                                    gACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCtt                                          >  6:186123/1‑148 (MQ=255)                                                                                                       cGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGaaa                                      <  4:71847/149‑1 (MQ=255)                                                                                                        gTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAAc                                     >  8:275912/1‑149 (MQ=255)                                                                                                           gCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCa                                   <  8:5505/148‑1 (MQ=255)                                                                                                               tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGcc                              <  7:191630/149‑1 (MQ=255)                                                                                                               tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGcc                              >  4:306093/1‑149 (MQ=255)                                                                                                               tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGcc                              >  7:156190/1‑149 (MQ=255)                                                                                                                 cTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCgg                            <  7:275912/149‑1 (MQ=255)                                                                                                                  ttGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                           >  8:83343/1‑149 (MQ=255)                                                                                                                  ttGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                           >  1:241718/1‑149 (MQ=255)                                                                                                                      gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCg                       <  3:306093/149‑1 (MQ=255)                                                                                                                       aaCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTc                        >  3:258175/1‑147 (MQ=255)                                                                                                                       aaCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGa                      <  8:32117/149‑1 (MQ=255)                                                                                                                        aCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTAGAAGCGTTGTATACGTGGGGGCATTACCAGTTCGAACGGCATTACCGCACCAGGGTTGCCGCTGCCGCTGAGCTTGAAACCATTGCGGTATTc                        >  1:101066/1‑147 (MQ=255)                                                                                                                          aTCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGt                   <  6:333716/149‑1 (MQ=255)                                                                                                                            cTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTg                 >  2:174033/1‑149 (MQ=255)                                                                                                                                   cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCAtt                                 >  6:180111/1‑126 (MQ=255)                                                                                                                                   cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCAtt                                 <  5:180111/126‑1 (MQ=255)                                                                                                                                      ggaaTTACTGTAATAGTAGCCTGCACGATATTCCCCTTT‑CAAGCCTTGGATATGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTTAAACCATTGCCGGATTCGAGGTTGGAAGGGTTc        <  3:217465/145‑1 (MQ=255)                                                                                                                                      ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCt       >  5:166613/1‑149 (MQ=255)                                                                                                                                       gCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTg      >  1:361722/1‑149 (MQ=255)                                                                                                                                        cATTAATGTAATACTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTg      <  2:101066/148‑1 (MQ=255)                                                                                                                                         aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCt    <  3:185426/149‑1 (MQ=255)                                                                                                                                              tGTAATAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGACGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCg  >  7:206926/1‑146 (MQ=255)                                                                                                                                                  |                                                                                                                                               GTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGA‑ACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTT‑CAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCG  >  NC_002947/1647783‑1648068

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