Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	65,820	G→A	6.4%	G293D (GGT→GAT)	betA‑I →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	65,820	0	G	A	6.4%	188.2 / 2.3	62	G293D (GGT→GAT)	betA‑I	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base A (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCGGGCCGCTGGCCATGG  >  NC_002947/65677‑65966                                                                                                                                                |                                                                                                                                                   gCATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCTGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATTCCAGGCGTtggtgg                                                                                                                                                 <  5:96906/147‑1 (MQ=255) gCATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTtggtgg                                                                                                                                                 <  5:92303/147‑1 (MQ=255) gCATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGc                                                                                                                                                >  5:817/1‑148 (MQ=255) gCATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGc                                                                                                                                                <  2:191699/148‑1 (MQ=255)  cATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCaa                                                                                                                                              >  3:295960/1‑149 (MQ=255)   aTGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAAc                                                                                                                                             >  7:190030/1‑149 (MQ=255)   aTGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAAc                                                                                                                                             >  2:3713/1‑149 (MQ=255)    tGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACt                                                                                                                                            >  5:116725/1‑149 (MQ=255)          tGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAg                                                                                                                                       <  6:817/148‑1 (MQ=255)           gcggcgCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGa                                                                                                                                     >  2:21797/1‑149 (MQ=255)               cgcgcCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGaccac                                                                                                                                 <  5:118875/149‑1 (MQ=255)                gcgcCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGaccacc                                                                                                                                <  1:152410/149‑1 (MQ=255)                   gcgcgAGAGCTTCCTGTGTGCTGGGGCGGTAGGTTCAGCGGGCATCCTCCAGCGCTCCGGCATCGGCGCGCGCCAGTTGGTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGct                                                                                                                                        <  6:261816/137‑2 (MQ=255)                   ccgcGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGCCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCa                                                                                                                                        >  5:261816/1‑138 (MQ=255)                    cgcgAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGc                                                                                                                             >  1:62875/1‑148 (MQ=255)                           tcatcCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTgcgc                                                                                                                     <  1:3713/149‑1 (MQ=255)                            catcCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGc                                                                                                                             <  5:190077/140‑1 (MQ=255)                            catcCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGc                                                                                                                             >  6:190077/1‑140 (MQ=255)                              tcCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACACGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTg                                                                                                                  <  4:300719/149‑1 (MQ=255)                                  ctGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGaccacc                                                                                                                                <  4:301684/131‑1 (MQ=255)                                  ctGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGaccacc                                                                                                                                >  3:301684/1‑131 (MQ=255)                                     tgCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTAc                                                                                                            <  6:116725/148‑1 (MQ=255)                                       cTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCa                                                                                                          <  1:21797/148‑1 (MQ=255)                                                         cGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTAcc                                                                                       >  4:268837/1‑149 (MQ=255)                                                              aTCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGaa                                                                                  <  6:231369/149‑1 (MQ=255)                                                               tcctccAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAAc                                                                                 >  5:209696/1‑149 (MQ=255)                                                                cctccAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCtt                                                                                     <  1:241264/144‑1 (MQ=255)                                                                cctccAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCtt                                                                                     >  2:241264/1‑144 (MQ=255)                                                                 ctccAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAAc                                                                                 >  4:190934/1‑147 (MQ=255)                                                                      gcgcTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATgg                                                                          <  3:268837/149‑1 (MQ=255)                                                                            cGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCaa                                                                      >  3:303045/1‑147 (MQ=255)                                                                             ggCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAg                                                                   >  7:156671/1‑149 (MQ=255)                                                                              gCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGc                                                                  >  4:283773/1‑149 (MQ=255)                                                                               cATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGcc                                                                 <  3:283773/149‑1 (MQ=255)                                                                                  cGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCtgt                                                              >  4:268918/1‑149 (MQ=255)                                                                                    gCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTgg                                                            >  7:10374/1‑149 (MQ=255)                                                                                        gcgcCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCa                                                        <  2:62875/149‑1 (MQ=255)                                                                                              gTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGatgg                                                  <  2:239458/149‑1 (MQ=255)                                                                                                 gCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGatgggc                                                <  4:295960/148‑1 (MQ=255)                                                                                                   tGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACa                                                                    <  1:213114/126‑1 (MQ=255)                                                                                                   tGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACa                                                                    >  2:213114/1‑126 (MQ=255)                                                                                                     gAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGatgggcatggg                                           >  1:324123/1‑149 (MQ=255)                                                                                                           ttGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCg                                     >  7:69159/1‑149 (MQ=255)                                                                                                           ttGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCg                                     >  1:182770/1‑149 (MQ=255)                                                                                                               gTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTAc                                 <  3:190934/149‑1 (MQ=255)                                                                                                                     gtgtTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTgcgc                                                                                                                     >  4:77456/1‑59 (MQ=255)                                                                                                                     gtgtTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTgcgc                                                                                                                     <  3:77456/59‑1 (MQ=255)                                                                                                                     gtgtTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGaa                                                                                  <  5:160185/94‑1 (MQ=255)                                                                                                                     gtgtTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGaa                                                                                  >  6:160185/1‑94 (MQ=255)                                                                                                                     gtgtTCGCAACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGAGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTAc                           >  4:241772/1‑149 (MQ=255)                                                                                                                           gCCACGACATGCCAGGCGTTGATGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCg                      >  7:263182/1‑148 (MQ=255)                                                                                                                           gCCACGACATGCCAGGCGTTGATGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGAGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCg                      >  7:258387/1‑148 (MQ=255)                                                                                                                            ccACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCg                      >  3:274537/1‑147 (MQ=255)                                                                                                                            ccACGACATGCCAGGCGTTGATGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTa                    <  8:258387/149‑1 (MQ=255)                                                                                                                            ccACGACATGCCAGGCGTTGATGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCAAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTa                    <  8:263182/149‑1 (MQ=255)                                                                                                                             cACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGatgggc                                                <  2:48752/120‑1 (MQ=255)                                                                                                                             cACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGatgggc                                                >  1:48752/1‑120 (MQ=255)                                                                                                                                gACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCgg                <  3:125999/149‑1 (MQ=255)                                                                                                                                gACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCgg                >  4:125999/1‑149 (MQ=255)                                                                                                                                       cAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGTGGCAAGATTGGCATGGGCCTGCGt                                    <  5:221068/122‑1 (MQ=255)                                                                                                                                       cAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGTGCATGGGCCTGCGt                                    >  6:221068/1‑122 (MQ=255)                                                                                                                                        aGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCGGGCCGCTgg        >  6:9222/1‑149 (MQ=255)                                                                                                                                        aGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCGGGCCGCTgg        >  6:40021/1‑149 (MQ=255)                                                                                                                                        aGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCGGGCCGCTgg        >  4:297922/1‑149 (MQ=255)                                                                                                                                               ggtggCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCGGGCCGCTGGCCATgg  <  3:241772/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GCATGAGTTTGCGGCGCGCCGCGAGATCATCCTCTGTGCTGGGGCGGTAGGTTCACCGGGCATCCTCCAGCGCTCCGGCATCGGCCCGCGCCAGTTGCTGGAAAGCTTGGGTATTGGTGTTCGCCACGACATGCCAGGCGTTGGTGGCAACTTGCAGGACCACCTGCAACTGCGCCTGATCTACCAGATCCGCAACACCCGTACCTTGAACCAGATGGCCAACAGCCTGTGGGGCAAGATGGGCATGGGCCTGCGTTACCTCTACGACCGTAGCGGGCCGCTGGCCATGG  >  NC_002947/65677‑65966

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