Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,013,216	G→A	6.5%	G151E (GGG→GAG)	PP_2634 →	cellulose synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,013,216	0	G	A	6.5%	227.6 / 3.6	76	G151E (GGG→GAG)	PP_2634	cellulose synthase
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (2/3);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

AGGTCGAGCTGGCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGCGACGAAGCCTTCCTC  >  NC_002947/3013078‑3013364                                                                                                                                            |                                                                                                                                                      aGGTCGAGCTGGCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCg                                                                                                                                            <  1:88111/149‑1 (MQ=255)  ggTCGAGCTGGCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGc                                                                                                                                           <  1:407463/149‑1 (MQ=255)    tCGAGCTGGCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTc                                                                                                                                         <  3:57955/149‑1 (MQ=255)      gAGCTGGCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGAGCGGACGTcgcccga                                                                                                                                        <  4:96675/149‑7 (MQ=255)            gCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGAGCGGACGT‑CGCCCGAGCagca                                                                                                                                 <  1:9768/149‑4 (MQ=255)            gCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGAGCGGACGT‑CGCCCGAGCagc                                                                                                                                  <  3:321178/148‑3 (MQ=255)             ccGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGAGCGGACGT‑CGCCCGAGCagcag                                                                                                                                >  3:146221/1‑145 (MQ=255)              cGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATAACACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTc                                                                                                                                >  4:361184/1‑148 (MQ=255)                   aTTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACAtg                                                                                                                          >  7:177036/1‑149 (MQ=255)                       cGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACAtgctgc                                                                                                                      >  4:71645/1‑149 (MQ=255)                        gACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACAtgctgc                                                                                                                      <  4:69864/148‑1 (MQ=255)                         aCTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGctgctg                                                                                                                    <  3:329177/149‑1 (MQ=255)                          cTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTgg                                                                                                                   <  2:358406/149‑1 (MQ=255)                             gCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACg                                                                                                                                                >  8:175865/1‑117 (MQ=255)                             gCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACg                                                                                                                                                <  7:175865/117‑1 (MQ=255)                             gCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTc                                                                                                                 >  1:15370/1‑148 (MQ=255)                              cTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCaa                                                                                                               <  1:238236/149‑1 (MQ=255)                               tgttgtACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAAc                                                                                                              >  5:92433/1‑149 (MQ=255)                                gttgtACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCGGCTggccgggt                                                                                                             >  2:47081/1‑143 (MQ=255)                                    tACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTa                                                                                                         <  8:168298/149‑1 (MQ=255)                                      cGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATg                                                                                                       <  2:73160/149‑1 (MQ=255)                                          cTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGATGGTCAACCTCTATGAccc                                                                                                   >  8:79456/1‑149 (MQ=255)                                             ccGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGc                                                                                                <  3:198788/149‑1 (MQ=255)                                                      ccGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCtt                                                                                       >  5:358694/1‑149 (MQ=255)                                                      ccGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCtt                                                                                       >  3:381870/1‑149 (MQ=255)                                                           cGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCa                                                                                    >  6:226030/1‑147 (MQ=255)                                                               accaccGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCg                                                                                                                                        <  6:428400/91‑1 (MQ=255)                                                               accaccGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCg                                                                                                                                        >  5:428400/1‑91 (MQ=255)                                                                 caccGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAATGGGACAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGAc                                                                             >  1:119821/1‑149 (MQ=255)                                                                 caccGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGAcc                                                                            >  2:51000/1‑149 (MQ=255)                                                                 caccGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGAcc                                                                            >  6:371767/1‑149 (MQ=255)                                                                  accGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCt                                                                           >  3:347824/1‑149 (MQ=255)                                                                    cGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGt                                                                         <  1:51000/149‑1 (MQ=255)                                                                      gCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCg                                                                                                            <  5:380878/112‑1 (MQ=255)                                                                      gCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCg                                                                                                            >  6:380878/1‑112 (MQ=255)                                                                      gCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGtt                                                                        >  5:284864/1‑148 (MQ=255)                                                                      gCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTg                                                                       <  1:51807/149‑1 (MQ=255)                                                                       cAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGa                                                                      >  4:401669/1‑149 (MQ=255)                                                                            gCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGgt                                                                  <  8:93344/148‑1 (MQ=255)                                                                                 gagaCTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGc                                                             <  2:15370/148‑1 (MQ=255)                                                                                 gagaCTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGc                                                             <  2:323595/148‑1 (MQ=255)                                                                                 gagaCTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGc                                                             <  6:92433/148‑1 (MQ=255)                                                                                 gagaCTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGc                                                             <  5:371767/148‑1 (MQ=255)                                                                                  agaCTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGt                                                           <  6:284864/149‑1 (MQ=255)                                                                                     cTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGAc                                                        <  4:347824/149‑1 (MQ=255)                                                                                      tGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATg                                                                    >  2:321054/1‑136 (MQ=255)                                                                                      tGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACt                                                       >  2:282808/1‑149 (MQ=255)                                                                                      tGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCGGTTGacg                                                                    >  2:315281/1‑134 (MQ=255)                                                                                          gTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGTCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGtt                                                   >  7:184258/1‑149 (MQ=255)                                                                                           tACGGGTACTGGTGGCAGAAGCCGAAAT‑GGTCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGttt                                                  <  8:184258/149‑1 (MQ=255)                                                                                           tACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGttt                                                  <  2:426675/149‑1 (MQ=255)                                                                                             cGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTc                                                 <  6:358694/148‑1 (MQ=255)                                                                                             cGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCg                                                <  4:381870/149‑1 (MQ=255)                                                                                             cGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCg                                                <  3:361184/149‑1 (MQ=255)                                                                                              gggTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGt                                               <  3:401669/149‑1 (MQ=255)                                                                                                gTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGa                                             >  8:46840/1‑149 (MQ=255)                                                                                                      gtggCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCtt                                       >  8:263669/1‑149 (MQ=255)                                                                                                      gtggCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCtt                                       >  7:327218/1‑149 (MQ=255)                                                                                                       tggCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTg                                      >  6:296641/1‑149 (MQ=255)                                                                                                       tggCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTg                                      >  2:9342/1‑149 (MQ=255)                                                                                                        ggCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGc                                     >  6:298082/1‑149 (MQ=255)                                                                                                           aGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGa                                  >  1:339545/1‑149 (MQ=255)                                                                                                               gCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATc                                >  5:357546/1‑147 (MQ=255)                                                                                                                   aaaT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAggtg                          <  5:226030/149‑1 (MQ=255)                                                                                                                    aaT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAggtgg                         <  5:296641/149‑1 (MQ=255)                                                                                                                    aaT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAggtgg                         <  1:9342/149‑1 (MQ=255)                                                                                                                     aT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAggtggt                        <  7:46840/149‑1 (MQ=255)                                                                                                                         gccagccaCCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCa                     >  4:327317/1‑149 (MQ=255)                                                                                                                         gccagccaCCTGCAGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTcc                      <  2:389743/148‑1 (MQ=255)                                                                                                                            agccaCCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATcg                  <  5:298082/149‑1 (MQ=255)                                                                                                                            agccaCCTGCTGCTCGGGA‑GACGTCCGCGCGAGCACTTCGACATGCTGCTGGTCAACCGCTAAGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTTGTCCATcg                  <  7:79456/149‑1 (MQ=255)                                                                                                                                aCCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGcgac              >  7:90202/1‑149 (MQ=255)                                                                                                                                   tgctgcTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGT‑ATCGCTTGCCGATCCAGGTGGTCCATCGCGACGAAg          >  8:340279/1‑149 (MQ=255)                                                                                                                                     gtgctCGAGCGGACGT‑CGCCCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGCGACGAAGc         >  4:33392/2‑149 (MQ=255)                                                                                                                                     ctgctCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGCGACGaa           <  6:179971/147‑1 (MQ=255)                                                                                                                                        ctCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGCGACGAAGCCtt      >  2:345694/1‑149 (MQ=255)                                                                                                                                            ggC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGCGACGAAGCCTTCCTc  >  5:65357/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                      AGGTCGAGCTGGCCGAGCATTTCGACTGGCTGTTGTACGACCTGCCGCAATTGCCGGGCGATACCACCGGCAATGGCATCGAGACTGATGTACGGGTACTGGTGGCAGAAGCCGAAAT‑GGCCAGCCACCTGCTGCTCGGGC‑GACGTCCGCTCGAGCACTTCGACATGCTGCTGGTCAACCGCTATGACCCGGCCAGCCGCTTGCAGAGCGACCTGTTGATGGTGTGGCGTGACTTGTTTCGTGATCGCTTGCCGATCCAGGTGGTCCATCGCGACGAAGCCTTCCTC  >  NC_002947/3013078‑3013364

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