Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,049,742	A→C	8.1%	S765R (AGC→CGC)	quiA →	quinate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,049,742	0	A	C	8.1%	181.4 / 6.9	62	S765R (AGC→CGC)	quiA	quinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base C (3/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGG  >  NC_002947/4049596‑4049889                                                                                                                                                   |                                                                                                                                                    cGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCC‑ACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGcc                                                                                                                                                  <  5:294362/149‑1 (MQ=255)      aaGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAggg                                                                                                                                              <  3:204655/149‑1 (MQ=255)              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)                                                                                                                                                   |                                                                                                                                                    CGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGGCGGCCCGATGACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGG  >  NC_002947/4049596‑4049889

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