Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,048,854	T→C	5.6%	F469L (TTC→CTC)	quiA →	quinate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,048,854	0	T	C	5.6%	207.5 / 2.1	72	F469L (TTC→CTC)	quiA	quinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (34/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 = 8.16e-01

TCCTACGACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACT  >  NC_002947/4048707‑4049000                                                                                                                                                    |                                                                                                                                                    tCCTACGACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGgtgt                                                                                                                                                     <  2:204591/148‑1 (MQ=255)  ccTACGACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTc                                                                                                                                                   >  2:11580/1‑149 (MQ=255)   cTACGACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTc                                                                                                                                                   >  7:215401/1‑148 (MQ=255)    tACGACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTc                                                                                                                                                   >  4:27852/1‑147 (MQ=255)       gACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGAc                                                                                                                                             >  8:234277/1‑149 (MQ=255)        aCCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGAcc                                                                                                                                            <  3:296157/149‑1 (MQ=255)           cGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTg                                                                                                                                         <  7:86463/149‑1 (MQ=255)              cAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGcac                                                                                                                                      <  2:120723/149‑1 (MQ=255)               aaTGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGcaca                                                                                                                                     <  5:48454/149‑1 (MQ=255)               aaTGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGcac                                                                                                                                      >  6:48454/1‑148 (MQ=255)                   aaCACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGa                                                                                                                                 >  6:132933/1‑149 (MQ=255)                     cacCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGtt                                                                                                                                                    >  7:296633/1‑129 (MQ=255)                     cacCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGtt                                                                                                                                                    <  8:296633/129‑1 (MQ=255)                      acCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACct                                                                                                                              >  4:231754/1‑149 (MQ=255)                       ccGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACctc                                                                                                                             >  5:70363/1‑149 (MQ=255)                        cGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACCACGACctct                                                                                                                            >  6:382521/1‑149 (MQ=255)                        cGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACctct                                                                                                                            >  1:195144/1‑149 (MQ=255)                          tgtTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTgg                                                                                                                          <  5:69165/149‑1 (MQ=255)                            ttCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTCACAGACCGTGCACAACGACCTCTgg                                                                                                                          <  7:234277/148‑1 (MQ=255)                            ttCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGa                                                                                                                        >  1:166877/1‑149 (MQ=255)                                tGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTc                                                                                                                    >  1:189960/1‑149 (MQ=255)                                   cGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGAc                                                                                                                 <  8:145552/149‑1 (MQ=255)                                     aTGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCt                                                                                                               <  3:231754/149‑1 (MQ=255)                                          atcctcctccACCGACAACTACGGCTTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGCCTTAGGACGCCACCACCGGCAACCAGAAAAGCGTGTTC‑CAGACCTTGCACAAAGACCTATGGGACTTCGACCTGCCCa                                                                                                          >  8:135190/5‑149 (MQ=255)                                           tCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCAt                                                                                                         >  5:379623/1‑149 (MQ=255)                                                 tccACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGc                                                                                                    >  3:117774/1‑148 (MQ=255)                                                     ccGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGc                                                                                               >  8:43045/1‑149 (MQ=255)                                                      cGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACg                                                                                                                                  <  1:107476/113‑1 (MQ=255)                                                      cGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACg                                                                                                                                  >  2:107476/1‑113 (MQ=255)                                                          aTCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCaa                                                                                                 <  6:34824/142‑1 (MQ=255)                                                          aTCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCaa                                                                                                 >  5:34824/1‑142 (MQ=255)                                                          aTCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGAt                                                                                          <  4:117774/149‑1 (MQ=255)                                                               cGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATc                                                                                         <  2:166877/145‑1 (MQ=255)                                                                ggCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACtt                                                                                    >  7:72722/1‑149 (MQ=255)                                                                  cGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCa                                                                                  <  6:313122/149‑1 (MQ=255)                                                                       aGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGAc                                                                                                                 >  7:89430/1‑113 (MQ=255)                                                                       aGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGAc                                                                                                                 <  8:89430/113‑1 (MQ=255)                                                                       aGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAg                                                                             <  6:70363/149‑1 (MQ=255)                                                                         cgcAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGga                                                                           >  2:345389/1‑149 (MQ=255)                                                                                aaCTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGc                                                                    >  3:335830/1‑149 (MQ=255)                                                                                 aCTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAgg                                                                            <  4:207305/140‑1 (MQ=255)                                                                                 aCTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAgg                                                                            >  3:207305/1‑140 (MQ=255)                                                                                     gACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTgg                                                                                                                          <  6:121023/90‑1 (MQ=255)                                                                                     gACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTgg                                                                                                                          >  5:121023/1‑90 (MQ=255)                                                                                     gACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGcc                                                                   >  3:12622/1‑145 (MQ=255)                                                                                     gACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGcc                                                                   <  4:12622/145‑1 (MQ=255)                                                                                       ccACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGACCGACTTCACCAAGGACGACGGCCAGTCgg                                                             >  8:282655/1‑149 (MQ=255)                                                                                        cacaCTTTCGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTTGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGt                                                            <  6:379623/149‑1 (MQ=255)                                                                                        cacaCTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGt                                                            <  5:132933/149‑1 (MQ=255)                                                                                           acTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGcc                                                         >  6:118709/1‑149 (MQ=255)                                                                                           acTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGCTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGcc                                                         >  4:127017/1‑149 (MQ=255)                                                                                           acTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGCTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGcc                                                         >  6:219902/1‑149 (MQ=255)                                                                                                cGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCAAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCg                                                     <  2:189960/148‑1 (MQ=255)                                                                                                      ttCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGa                                              >  3:214176/1‑149 (MQ=255)                                                                                                        cGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGAt                                             <  6:82614/148‑1 (MQ=255)                                                                                                          gTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCgg                                          >  3:159529/1‑149 (MQ=255)                                                                                                             cTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCa                                        <  3:27852/148‑1 (MQ=255)                                                                                                              tGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTTGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCggtg                                                  <  6:277940/137‑1 (MQ=255)                                                                                                              tGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCggtg                                                  >  5:277940/1‑137 (MQ=255)                                                                                                               ggCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGATC‑CAGACCGTGCACAACGACCTCTCGGACTACGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGAGCCTCCGGTGCTGATCGGCAcc                                      >  3:4970/1‑148 (MQ=255)                                                                                                                 cTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAg                                   <  5:382521/149‑1 (MQ=255)                                                                                                                      gACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACGTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCg                               <  1:11580/148‑1 (MQ=255)                                                                                                                            accaccGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGAt                        <  8:72722/149‑1 (MQ=255)                                                                                                                                ccGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTAc                    <  4:335830/149‑1 (MQ=255)                                                                                                                                  ggCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGt                  <  4:159529/149‑1 (MQ=255)                                                                                                                                        cAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGa            >  5:169328/1‑149 (MQ=255)                                                                                                                                        cAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGCCGATCaccaccgcagg                                 <  6:200246/128‑6 (MQ=255)                                                                                                                                        cAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGCCGATCaccaccgcagg                                 >  5:200246/1‑123 (MQ=255)                                                                                                                                          gAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCaa                                    <  2:256651/123‑1 (MQ=255)                                                                                                                                          gAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGAACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCaa                                    >  1:256651/1‑123 (MQ=255)                                                                                                                                          gAAATGGGTGCTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGAcc          <  5:219902/149‑1 (MQ=255)                                                                                                                                          gAAATGGGTGCTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGAcc          <  3:127017/149‑1 (MQ=255)                                                                                                                                           aaaTGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTAc                    >  7:254725/1‑138 (MQ=255)                                                                                                                                           aaaTGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTAc                    <  8:254725/138‑1 (MQ=255)                                                                                                                                            aaTGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGACCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACcgc        <  7:282655/149‑1 (MQ=255)                                                                                                                                                  tgtTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGTAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACt  >  2:363099/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCCTACGACCCGGCAATGAACACCGTGTTCCTGCCGATGGGCTCGTCCTCCACCGACATCTACGGCGTTGAGCGCAGCAAACTGGACCACACTTACGGCGCTTCGGTGCTGGCTTTGGACGCCACCACCGGCAACCAGAAATGGGTGTTC‑CAGACCGTGCACAACGACCTCTGGGACTTCGACCTGCCCATGCAGCCAAGCCTGATCGACTTCACCAAGGACGACGGCCAGTCGGTGCCTGCGGTGGTGATCGGCACCAAGGCCGGGCAGATCTACGTGCTCGACCGCGCCACT  >  NC_002947/4048707‑4049000

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