Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	42.4%	97.0 / 14.5	59	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (6/28);  new base G (25/0);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCG  >  NC_000913/354142‑354429                                                                                                                                                   |                                                                                                                                               gggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGaa                                                                                                                                             >  1:430643/1‑150 (MQ=255)     aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGATg                                                                                                                                        >  1:427983/1‑151 (MQ=255)       ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                                                                       <  1:467197/150‑1 (MQ=255)       ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGATGAAAGATgc                                                                                                                                       <  1:33799/150‑1 (MQ=255)        cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGGGGGGGGGAAAATggg                                                                                                                                      >  2:306522/1‑148 (MQ=255)         gTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATgcg                                                                                                                                      >  2:189907/1‑149 (MQ=255)             cGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGaaa                                                                                                                                <  2:323842/151‑1 (MQ=255)              gCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑g                                                                                                                              <  1:126845/151‑1 (MQ=255)                 gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑gg                                                                                                                             <  2:604380/149‑1 (MQ=255)                  aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGGGCCGGAA‑gggg                                                                                                                           >  2:825140/1‑149 (MQ=255)                   cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAA‑ggg                                                                                                                            >  1:41387/1‑148 (MQ=255)                    ggTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAA‑GGGCAgg                                                                                                                        >  1:739048/1‑151 (MQ=255)                        aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAggtggc                                                                                                                    <  1:755896/151‑1 (MQ=255)                              ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAggtggcggtgg                                                                                                               <  1:425033/150‑1 (MQ=255)                                aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGAAAGAAGAGCTGGAA‑ggggcggggggggggg                                                                                                               >  2:780115/1‑135 (MQ=255)                                     cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAA‑GGGGAGGGGGGGGGGGGGGTTGt                                                                                                        >  2:266732/1‑150 (MQ=255)                                      gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCCCTTTAG‑GGGCGGGTGGCGGAGGCGGTTggtc                                                                                                      >  2:547165/1‑148 (MQ=255)                                        gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCAtt                                                                                                    <  2:1098478/151‑1 (MQ=255)                                         cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCGGAAA‑GGGGGGGGGGCGGGGGGGTTTTTTAtt                                                                                                    >  2:7810/1‑150 (MQ=255)                                          gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTcc                                                                                                  <  2:48316/151‑1 (MQ=255)                                           cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGGGGCGGTTGTCATTcc                                                                                                  >  1:1003944/1‑150 (MQ=255)                                               gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGaaa                                                                                              <  1:90119/150‑1 (MQ=255)                                                aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAaga                                                                                            <  2:300734/151‑1 (MQ=255)                                                   aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagagag                                                                                         <  1:74330/151‑1 (MQ=255)                                                    tGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGGGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagagag                                                                                         >  1:494714/1‑150 (MQ=255)                                                     gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGc                                                                                        <  1:372980/150‑1 (MQ=255)                                                     gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCg                                                                                       <  2:908327/151‑1 (MQ=255)                                                             cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTg                                                                               <  1:1068152/151‑1 (MQ=255)                                                                tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa                                                                             <  2:641301/150‑1 (MQ=255)                                                                tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                              >  1:810261/1‑149 (MQ=255)                                                                    gggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAAAGCGGTGAAAGGGGCAGG‑GGCGGGGGCGTTTGTTGTTTCCAAagagaggggaaggtcgg                                                                              >  2:618554/1‑134 (MQ=255)                                                                       cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAA‑GGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGTg                                                                     >  1:74440/1‑151 (MQ=255)                                                                          cGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTGTGTCATTCCGAAAGAGGGCGACAGTCTGGAAGATCGTGAtg                                                                  <  2:85756/151‑1 (MQ=255)                                                                              tGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGGTGTgg                                                               >  1:424442/1‑150 (MQ=255)                                                                               ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgc                                                              <  1:316141/150‑1 (MQ=255)                                                                                      cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                      <  2:559933/151‑1 (MQ=255)                                                                                      cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAGGGGCTGAAA‑AGGGAGGGGGGGGTGGGGGGTGTCATTTCGTAAGAGAGCGACAGTTCGGAAGGGTGGGATGTgg                                                               >  2:574816/1‑142 (MQ=255)                                                                                        cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAGATTGGCTGGAA‑GGGGAAGTGGCGGGGGCGGTTTTTATTTCGAAAGAGAGGGGCAGTCTGGGAAATCGTGGTGGGGGGCACTCg                                                       >  2:963287/1‑148 (MQ=255)                                                                                         gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAGATGGGCTGGAA‑GGGGGGGGGGCGGGGGGGTTTGTCATTTCCAAAGAGGGGGAAAGTCTGGGAGATTGTGATGGGGGGCa                                                           >  2:592298/1‑143 (MQ=255)                                                                                              tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGaa                                               <  2:416709/150‑1 (MQ=255)                                                                                              tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGTAAAATTCGCTGAAA‑GGGGGGGGGGCGGGGGCGTTTGTCATTTCGAAAGAGGGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                >  2:358214/1‑149 (MQ=255)                                                                                                aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGc                                            <  2:234699/151‑1 (MQ=255)                                                                                                aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGc                                            <  1:749209/151‑1 (MQ=255)                                                                                                 agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                                >  2:580298/1‑146 (MQ=255)                                                                                                      gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCCCTCGCAAGAGAAGGCGATTa                                       <  2:163198/150‑1 (MQ=255)                                                                                                           tcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGGCTGGGAGATCGAGATGTGGGGCACTCGCAAGAGAAGGCGATTTTGgcg                                  >  2:781083/1‑150 (MQ=255)                                                                                                             cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgc                                   <  2:606753/147‑1 (MQ=255)                                                                                                                ccATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCTGAAA‑GGGCAGGGGGCGGTGGCGTTTGTCATTCCGAAAGAGAGGGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGCTTATGGCGCtggtg                            >  2:213297/1‑151 (MQ=255)                                                                                                                ccATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtg                            >  1:520516/1‑151 (MQ=255)                                                                                                                 cATCCGGGCATTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCCGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGGGGGACAGTCTGGGAGAACGGGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                           >  1:602483/1‑151 (MQ=255)                                                                                                                    ccGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGGGGGGGTGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGGTTATGGCGGtggtgg                           >  2:867773/1‑148 (MQ=255)                                                                                                                     cGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGGGGTGGCGTTTGTCATTCCGAAAGAGAGCGGCCGTCGGGAAGATAGTGATGTGGCGCACTCGCCAGAGAAGGGGGTTATGGCGCtggggg                           >  1:458306/1‑147 (MQ=255)                                                                                                                        gCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCAAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCa                    <  1:271364/151‑1 (MQ=255)                                                                                                                        gCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGcc                     >  1:927949/1‑150 (MQ=255)                                                                                                                            tGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGGAA‑GGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGGGACAGTCTGGAAGATCGGGGTGTGGCGCACTCGGAAGAGAAGGGGGTTATGGCGCTGGTGGGCAGCCAGAtt                >  2:16654/1‑151 (MQ=255)                                                                                                                             gCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGAtt                <  2:633919/150‑1 (MQ=255)                                                                                                                                      gCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg      <  2:1053068/151‑1 (MQ=255)                                                                                                                                        ggtggtTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgc    >  1:569156/1‑151 (MQ=255)                                                                                                                                           ggtTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg  <  1:921640/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCG  >  NC_000913/354142‑354429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑