Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,709,378	G→C	5.6%	R1132R (CGG→CGC)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,709,378	0	G	C	5.6%	223.2 / 3.1	72	R1132R (CGG→CGC)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base C (2/2);  total (37/35)
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.58e-01

CGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGAC  >  NC_002947/5709231‑5709521                                                                                                                                                     |                                                                                                                                                 cgttGCAGCCTTCGAAAACGACATGCTGCTGTGCCGCCGAGTGCTGGAGCAGGTCGACCACCAGGCGCTTAGC‑GCACAGTCCATCTGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACgg                                                                                                                                                  <  7:331630/145‑1 (MQ=255) cGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACGACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACggg                                                                                                                                                 <  4:231339/149‑1 (MQ=255)      cAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                            <  4:217194/149‑1 (MQ=255)      cAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCcgc                                                                                                                                            <  3:71612/149‑1 (MQ=255)        gCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGa                                                                                                                                          <  5:338154/149‑1 (MQ=255)              gacaaCTACATGCTGAGGTGCCGACGCTTGCTTGAGAAGGTCGACCACCCTGCGCTGAGC‑GCACAGGCCATCGGGCTGACCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATAGCGAAAGGGTCCGCGATTcca                                                                                                                                     <  6:221303/145‑1 (MQ=255)               aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGGACCCGTTCGCGAtaatcg                                                                                                                                    >  4:24287/1‑143 (MQ=255)               aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGGACCCGTTCGCGAtaatcg                                                                                                                                    <  3:24287/148‑6 (MQ=255)               aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                   >  5:44853/1‑149 (MQ=255)               aaaaCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagc                                                                                                                                   >  4:13407/1‑149 (MQ=255)                  aCGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTccagccag                                                                                                                                <  3:13407/149‑1 (MQ=255)                    gACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAgccaggc                                                                                                                              <  1:220855/149‑1 (MQ=255)                    gACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGACCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAgccaggc                                                                                                                              >  2:135304/1‑149 (MQ=255)                         gctgctGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgcc                                                                                                                         >  2:310262/1‑149 (MQ=255)                             ctgTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGCGGCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgcc                                                                                                                      >  1:263664/1‑149 (MQ=255)                             ctgTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAgccgc                                                                                                                       <  1:135304/147‑1 (MQ=255)                                  ccgccgCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCAGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGc                                                                                                               >  6:26950/1‑149 (MQ=255)                                     ccgcTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTg                                                                                                            >  3:219935/1‑149 (MQ=255)                                            cTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCgtca                                                                                                      <  5:19681/148‑1 (MQ=255)                                               gAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCgtcagt                                                                                                    >  2:38102/1‑147 (MQ=255)                                                 gCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGc                                                                                                >  4:350942/1‑149 (MQ=255)                                                    ggTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGc                                                                                             <  6:44853/149‑1 (MQ=255)                                                     gTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGAGGCt                                                                                            >  5:152511/1‑149 (MQ=255)                                                       cGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTg                                                                                                            >  1:356663/1‑131 (MQ=255)                                                       cGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTg                                                                                                            <  2:356663/131‑1 (MQ=255)                                                             ccaGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTggatcgga                                                                                    >  6:315861/1‑149 (MQ=255)                                                                    cTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGc                                                                             <  4:219935/149‑1 (MQ=255)                                                                     tGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCg                                                                            >  1:205517/1‑149 (MQ=255)                                                                      gAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCg                                                                            >  5:259434/1‑148 (MQ=255)                                                                       aGC‑GCACAGTCTATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCgg                                                                           <  1:310262/148‑1 (MQ=255)                                                                                gTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGc                                                                  <  5:315861/149‑1 (MQ=255)                                                                                 tCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTg                                                                   >  4:50702/1‑147 (MQ=255)                                                                                    aTCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGa                                                              >  7:114411/1‑149 (MQ=255)                                                                                        ggCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCAGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGAcccc                                                          <  5:26950/149‑1 (MQ=255)                                                                                             ccccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGAcccc                                                          >  1:147489/1‑144 (MQ=255)                                                                                             ccccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGa                                                      <  2:205517/148‑1 (MQ=255)                                                                                               ccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTgg                                                                                          >  3:274572/1‑110 (MQ=255)                                                                                               ccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTgg                                                                                          <  4:274572/110‑1 (MQ=255)                                                                                               ccGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATc                                                    >  7:27022/1‑148 (MQ=255)                                                                                                       cAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTccc                                           <  8:27022/149‑1 (MQ=255)                                                                                                         gagaCTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCaa                                         <  3:50702/149‑1 (MQ=255)                                                                                                           gaCTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCaaa                                        >  3:344101/1‑147 (MQ=255)                                                                                                              tGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCa                                                                 <  1:319628/120‑1 (MQ=255)                                                                                                              tGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCa                                                                 >  2:319628/1‑120 (MQ=255)                                                                                                                ccGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCTGGACCCCAAGATCGATTTTCCCAATAAAGcc                                  <  1:38102/149‑1 (MQ=255)                                                                                                                ccGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAg                                    <  1:149020/147‑1 (MQ=255)                                                                                                                  gAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACCGCAGGACCCCAAGATCGAtttt                                              <  6:21225/135‑1 (MQ=255)                                                                                                                  gAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACCGCAGGACCCCAAGATCGAtttt                                              >  5:21225/1‑135 (MQ=255)                                                                                                                        aTACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAAt                                        <  1:369009/135‑1 (MQ=255)                                                                                                                        aTACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATACTGGCGGAAACACAGCAGGACCCCAAGAAAGATTTTACCaaa                                        >  2:369009/1‑134 (MQ=255)                                                                                                                           ctgCATTACCTCGATTATCGCTAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        <  6:259434/146‑1 (MQ=255)                                                                                                                            aGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        <  6:137473/147‑1 (MQ=255)                                                                                                                            aGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        >  5:137473/1‑147 (MQ=255)                                                                                                                              cAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGGCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACGCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCaa                    <  8:211018/149‑1 (MQ=255)                                                                                                                              cAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGcc                      >  3:281969/1‑147 (MQ=255)                                                                                                                                 tACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        <  8:307711/142‑1 (MQ=255)                                                                                                                                 tACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAg                        >  7:307711/1‑142 (MQ=255)                                                                                                                                   ccTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGAcc                <  4:281969/148‑1 (MQ=255)                                                                                                                                    cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGc                                                                  <  2:68584/97‑1 (MQ=255)                                                                                                                                    cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGc                                                                  >  1:68584/1‑97 (MQ=255)                                                                                                                                    cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAt              >  4:117035/1‑149 (MQ=255)                                                                                                                                    cTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAt              >  7:311820/1‑149 (MQ=255)                                                                                                                                     tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATaaa                                     >  1:246927/1‑125 (MQ=255)                                                                                                                                     tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATaaa                                     <  2:246927/125‑1 (MQ=255)                                                                                                                                     tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             >  4:126730/1‑149 (MQ=255)                                                                                                                                     tCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             >  2:260065/1‑149 (MQ=255)                                                                                                                                      cGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             >  1:229752/1‑148 (MQ=255)                                                                                                                                       gATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCCGGCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             >  1:71587/1‑148 (MQ=255)                                                                                                                                         ttATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAAc         <  3:350942/149‑1 (MQ=255)                                                                                                                                            tCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             >  2:315524/1‑142 (MQ=255)                                                                                                                                            tCGCGGACCCGTTCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCAtt             <  1:315524/142‑1 (MQ=255)                                                                                                                                            tCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCACC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACg        >  4:232127/1‑147 (MQ=255)                                                                                                                                               cgAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATATTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGa   <  1:260065/149‑1 (MQ=255)                                                                                                                                                gAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGAc  <  6:152511/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CGGTGCAGCCATCGAAAACGACATGCTGCTGTGCCGCCGCTTGCTGGAGCAGGTCGACCACCAGGCGCTGAGC‑GCACAGTCCATCGGGCTGCCCCGCGTTCCAGAGACTGCCGAGCAGATACAGCAATACCTCGATTATCGCGAACGGGTCCGCGATTCCAGCCAGGCCAGCCGCC‑GCCTGCTTGTCGTCAGTCAGCGGCTGGATCGGATCCTGGCGGAAACACTGCAGGACCCCAAGATCGATTTTCCCAATAAAGCCACGAAACTAGCCAAGACCATTGAACGCAGGAC  >  NC_002947/5709231‑5709521

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