Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,049,751	G→T	7.7%	G768C (GGC→TGC)	quiA →	quinate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,049,751	0	G	T	7.7%	187.3 / 4.4	60	G768C (GGC→TGC)	quiA	quinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base T (3/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

CCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGC  >  NC_002947/4049609‑4049895                                                                                                                                               |                                                                                                                                                 ccTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGggcgg                                                                                                                                                 >  6:256194/1‑144 (MQ=255) ccTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGggcgg                                                                                                                                                 <  5:256194/144‑1 (MQ=255) ccTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGa                                                                                                                                            <  8:197385/149‑1 (MQ=255) ccTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGa                                                                                                                                            >  7:273400/1‑149 (MQ=255)   tGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCCGCCCTGGCTGccccc                                                                                                                                            >  2:234985/1‑145 (MQ=255)   tGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGAt                                                                                                                                           >  2:67716/1‑148 (MQ=255)   tGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCAACGACAGCAGCAACGGAAACCAAATCTGGAAGGCTCGCCTGCCCGTGGGCCGCCCTGGCTGcccac                                                                                                                                            <  1:234985/147‑3 (MQ=255)    gCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGAGGGCAGCCAGGGCGGCCCGATGa                                                                                                                                         >  2:413124/1‑149 (MQ=255)     ccAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGa                                                                                                                                         <  6:135662/148‑1 (MQ=255)     ccAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGAAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGAc                                                                                                                                        <  6:280759/149‑1 (MQ=255)              ggcggcACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGtctc                                                                                                                               >  8:69589/1‑149 (MQ=255)              ggcggcACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGtc                                                                                                                                 <  3:85832/147‑1 (MQ=255)               gcggcACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCg                                                                                                                              <  1:50514/149‑1 (MQ=255)               gcggcACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCg                                                                                                                              <  5:106703/149‑1 (MQ=255)                cggcACCCTGTCGACCCAAGGTGGCCTGGTGTTGATAGCAGGTACCGAGGACTTATACCTGCGCGCCTACGACAGCAGCTAGGGTAACGAAATCTGGTAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGAGCTACGTCTCGc                                                                                                                             <  8:38118/149‑1 (MQ=255)                cggcACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGc                                                                                                                             <  5:354821/149‑1 (MQ=255)                 ggcACCCTGTCGACCGAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGct                                                                                                                            <  8:184492/149‑2 (MQ=255)                    aCCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGGCCTACGTCTCGCCGaa                                                                                                                         <  6:112038/149‑1 (MQ=255)                    aCCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCg                                                                                                                           >  3:346767/1‑147 (MQ=255)                    aCCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCg                                                                                                                           >  6:46706/1‑147 (MQ=255)                    aCCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCg                                                                                                                           >  6:40413/1‑147 (MQ=255)                    aCCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGaa                                                                                                                         <  5:13114/149‑1 (MQ=255)                                         cTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCCGCCCTGGCTGcccac                                                                                                                                            >  2:45161/1‑107 (MQ=255)                                         cTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCCGCCCTGGCTGcccac                                                                                                                                            <  1:45161/109‑3 (MQ=255)                                            gtgtTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGa                                                                                                  >  6:45384/1‑148 (MQ=255)                                            gtgtTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGAt                                                                                                 <  7:69589/149‑1 (MQ=255)                                              gtTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCa                                                                                               <  8:273400/149‑1 (MQ=255)                                               ttCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGGGGTGATCAc                                                                                              >  7:342771/1‑149 (MQ=255)                                                   tCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGTAATCTGGAAGGCTCTCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCActgct                                                                                          <  2:177933/149‑1 (MQ=255)                                                      cAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTggc                                                                                       <  2:290326/149‑1 (MQ=255)                                                       aGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTggcg                                                                                      <  8:32303/149‑1 (MQ=255)                                                             ccAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGcgccc                                                                                >  2:289664/1‑149 (MQ=255)                                                                    ttCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTc                                                                         >  1:137612/1‑149 (MQ=255)                                                                      gcggatGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGa                                                                       <  4:199386/144‑1 (MQ=255)                                                                         ccTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGAc                                                                      <  5:45384/147‑1 (MQ=255)                                                                           tGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGAc                                                                  >  1:190576/1‑149 (MQ=255)                                                                              gcgcCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCGGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGt                                                               <  6:211944/149‑1 (MQ=255)                                                                                      gACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTa                                                       >  8:36926/1‑149 (MQ=255)                                                                                             gcaACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGAtt                                                >  6:284430/1‑149 (MQ=255)                                                                                                aCGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCt                                             >  2:149283/1‑149 (MQ=255)                                                                                                 cGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCtt                                            <  1:67716/149‑1 (MQ=255)                                                                                                    tAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACg                                         >  6:163898/1‑149 (MQ=255)                                                                                                     aaCGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGc                                        >  6:65861/1‑149 (MQ=255)                                                                                                      aCGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGct                                       >  2:122831/1‑149 (MQ=255)                                                                                                      aCGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGc                                        <  5:284430/148‑1 (MQ=255)                                                                                                       cGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGctc                                      >  6:417679/1‑149 (MQ=255)                                                                                                        gAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGctct                                     >  8:135214/1‑149 (MQ=255)                                                                                                         aaaTCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTg                                    <  5:65861/149‑1 (MQ=255)                                                                                                           aTCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGcc                                  <  1:254011/149‑1 (MQ=255)                                                                                                                gAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTaa                              <  5:417679/148‑1 (MQ=255)                                                                                                                gAAGGCTCGCCTGCCCGAGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCg                                 <  1:413124/145‑1 (MQ=255)                                                                                                                      tCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACcgcg                       >  5:21930/1‑149 (MQ=255)                                                                                                                          cTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGc                   <  5:163898/149‑1 (MQ=255)                                                                                                                            gCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTaaaa                            <  8:85957/138‑1 (MQ=255)                                                                                                                            gCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTaaaa                            >  7:85957/1‑138 (MQ=255)                                                                                                                             tcatcgggCCGCCCTGGCTGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTg                >  6:180266/6‑149 (MQ=255)                                                                                                                             cccGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTc                                              >  7:400085/1‑119 (MQ=255)                                                                                                                             cccGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTc                                              <  8:400085/119‑1 (MQ=255)                                                                                                                                     cAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCgg        <  1:149283/149‑1 (MQ=255)                                                                                                                                     cAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGGGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCgg        <  4:421680/149‑1 (MQ=255)                                                                                                                                        ccAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGc      >  5:318257/1‑148 (MQ=255)                                                                                                                                           gggCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGc  <  1:289664/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 CCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGC  >  NC_002947/4049609‑4049895

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