Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	48.4%	82.3 / 17.5	62	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (5/27);  new base G (30/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 7.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354144‑354437                                                                                                                                                 |                                                                                                                                                       gCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                  <  2:272832/151‑1 (MQ=255)     ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcg                                                                                                                                              <  2:394688/151‑1 (MQ=255)           cGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGGGCT‑Gaa                                                                                                                                        >  1:168757/1‑150 (MQ=255)                aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAggg                                                                                                                                    >  1:538246/1‑149 (MQ=255)                 cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAggg                                                                                                                                    >  1:199766/1‑148 (MQ=255)                 cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GGAAggg                                                                                                                                    >  1:118232/1‑148 (MQ=255)                      aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGCTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGCAGTGGCGGTGGGTGGGGGGGGAGGAGGGCC‑GAAAAGGCAgggg                                                                                                                              >  2:321682/1‑147 (MQ=255)                           tttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCT‑GGAAGGGGAggggggggggg                                                                                                                       >  2:226312/1‑151 (MQ=255)                             tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAggtggcggtggc                                                                                                                      <  2:449030/150‑1 (MQ=255)                              aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGt                                                                                                                    <  2:549120/151‑1 (MQ=255)                                tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGtt                                                                                                                   <  2:180092/150‑1 (MQ=255)                                    gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCa                                                                                                              <  2:136571/151‑1 (MQ=255)                                     ggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGCGTTTGTCAt                                                                                                             >  1:767435/1‑151 (MQ=255)                                       cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAGAGAGGGGCG‑GAAAAGGGAGGGGGCGGTGGGGGTTGTCa                                                                                                              >  1:225933/1‑148 (MQ=255)                                          aCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGGGGGGGGGGAAAATCGCT‑GAAAGGGGAggtgggggggggg                                                                                                                     >  2:69065/1‑136 (MQ=255)                                          aCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAAAAGCGCT‑GGAAGGGGAGGGGGGGGGGGCGGTTGTCAtt                                                                                                            >  1:63288/1‑147 (MQ=255)                                           cTGAGGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATTTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCTGGTGGCGGTGGCGTTTGTCATTCCgt                                                                                                        <  2:214306/150‑2 (MQ=255)                                           cTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGGT‑GGAAGGGGAggggggggggggggtt                                                                                                                  >  2:26377/1‑135 (MQ=255)                                              aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGATGCGCT‑GTAAGGGCAGGGGGCGGGGGGGGTTGTCAtt                                                                                                            >  2:361057/1‑143 (MQ=255)                                                  tGTAATTAACGTTGGCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTTTCTCCAACCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCg                                                                                                         <  2:182422/142‑1 (MQ=255)                                                   gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGc                                                                                                <  1:256957/150‑1 (MQ=255)                                                   gTAATGAACGGTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCGCCAGCCATCCGGGGGTTGCCGATGTGGCGGTGGTTGGGGTGAAAGATGCGGT‑GAAAGGGCAGGTGGCGGTGGCGTTTGCCATTCGGAAag                                                                                                     <  2:189674/145‑1 (MQ=255)                                                    tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCg                                                                                               <  2:180850/150‑1 (MQ=255)                                                    tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAGTGCGCT‑GAAAGGGCAGGTGGGGGGGGGGGTTGTCATTTCCCAagaga                                                                                                  >  1:180869/1‑147 (MQ=255)                                                     aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGGAGGGGGGGGGGGGGGTTGTCATTTCGgaa                                                                                                      >  1:488186/1‑142 (MQ=255)                                                      aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACa                                                                                            <  2:34757/151‑1 (MQ=255)                                                      aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAAATGCGCT‑GAAAGGGGAGGTGGCGGGGGGGTTTGTCATTCCGAAAAAGAGCg                                                                                               >  1:708971/1‑148 (MQ=255)                                                       ttAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCAt                                                                                                             <  1:551758/133‑1 (MQ=255)                                                         aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAGGCGCTGGAAAGGGC‑GGGGGCGGTGGGGCTTGTCATTTCGCCAGGGAGCggccgt                                                                                          >  2:187258/1‑145 (MQ=255)                                                           cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAAAGCGCT‑GAAAGGGGAGGGGGCGGGGGCGTTTGTTATTTCGaaa                                                                                                      >  2:560278/1‑136 (MQ=255)                                                           cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGGAGGGGGCGGGGGGGGTTTTCATTCCGAAAAAGAGCGGCAGGCTg                                                                                       >  2:137996/1‑151 (MQ=255)                                                            gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                      <  1:555974/151‑1 (MQ=255)                                                            gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCT‑GCAAAGGGACGTGGGGGGGGGGGTTTGCATTTCGCAAGAGAGCGGAAGTCTgg                                                                                      >  2:619421/1‑151 (MQ=255)                                                               gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGaa                                                                                    <  2:29162/150‑1 (MQ=255)                                                                 cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGGT‑GAAAGGGCCGGGGGCGGGGGCGTTTGTCATTTCGCAAGAGAGAGACAGTCt                                                                                        >  2:559757/1‑144 (MQ=255)                                                                     cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAATGGGCC‑GGAAGGGGAGGTGGGGGGGGGGGTTGTTCTTCCGAAAGAAAGCGACAGTCttggaggacgtg                                                                             >  2:750532/1‑140 (MQ=255)                                                                       tCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGAGGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCCATTCGAAAGAGAGAGGCAGGCTgggaaagcg                                                                               >  1:421189/1‑140 (MQ=255)                                                                        cGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GGAAGGGGAGGGGGGGGGGGCGGTTGTCCTTCCGGAAGAGAGCGGCAGTCTGGGAGATCGTGAt                                                                           >  1:460403/1‑150 (MQ=255)                                                                           cTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCC‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTTCCAAAGAGAGCGACCGGCTGGGAGATCggggtgggg                                                                       >  2:111718/1‑143 (MQ=255)                                                                               ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa                                                                   <  1:760994/151‑1 (MQ=255)                                                                                   gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGGGGGGGGGGTTGTCATTCCGCAAGAGAGCGACAGTCTGGGAGATCGTGGTGTGGCGCACTc                                                                >  2:470249/1‑150 (MQ=255)                                                                                   gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCT‑GAAAGGGCAGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGgcgccc                                                                  >  2:757227/1‑146 (MQ=255)                                                                                     gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCa                                                             <  2:462415/151‑1 (MQ=255)                                                                                       gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                                           <  2:788094/151‑1 (MQ=255)                                                                                         gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GGAAGGGGAGGTGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGGCAGTCTGGAAGATCGTGATGTGGCGCACTCGCaa                                                            >  2:805304/1‑148 (MQ=255)                                                                                          aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAAAGCGCT‑GGAAGGGGCGGGGGCGGGGGCGGTTGTCATTCCGAAAAAGAGGGGCAGGCTGGAAGATCGGGATGTGGGGCACTCGCAag                                                           >  2:574643/1‑148 (MQ=255)                                                                                           ttGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCT‑GCAAGGGCAGGGGGCGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGGGCACTc                                                                >  2:540040/1‑142 (MQ=255)                                                                                            tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGaa                                                       <  2:216729/150‑1 (MQ=255)                                                                                                   agTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTa                                               <  1:513636/151‑1 (MQ=255)                                                                                                    gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTAt                                              <  1:551385/151‑1 (MQ=255)                                                                                                          ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcg                                          <  2:356100/149‑1 (MQ=255)                                                                                                           cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGGAAGAGAGGGGCAGGCTGGAAGATCGGGATGGGGGGCACCCGCAAGAGAAGGCGATTATgggggcg                                       >  2:201625/1‑146 (MQ=255)                                                                                                            aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAACGATGCGCC‑GAAAGGGCAGGGGGCGGGGGCGTTTGTTCTTCCGAAAGAGAGCGACAGGCTGGCAGATCGTGATGTGGCGCACTCGCCAGAGAAGGCGATTATGGCGCtgg                                      >  2:71937/1‑151 (MQ=255)                                                                                                            aGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GGAAGGGCAGGGGGCGGGGGCGTTTGTTCTTTCCAAAGAGAGCGACAGTCTGGAAGATAGGGGTGTGGGGCCCTCGCAAGAGAAGGCGGTTATGgcgc                                         >  2:444002/1‑148 (MQ=255)                                                                                                                aTCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtg                                    <  2:216588/149‑1 (MQ=255)                                                                                                                           gCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTg                       <  2:2521/151‑1 (MQ=255)                                                                                                                               aaGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGCCAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCa                    >  2:498022/1‑150 (MQ=255)                                                                                                                                 gTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCaa                   <  2:620583/149‑1 (MQ=255)                                                                                                                                    gCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg              >  2:615097/1‑151 (MQ=255)                                                                                                                                        tggtTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc           >  2:23534/1‑150 (MQ=255)                                                                                                                                              gggTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgc    <  2:279100/151‑1 (MQ=255)                                                                                                                                              gggTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgc    <  1:772224/151‑1 (MQ=255)                                                                                                                                                 tGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  >  2:659193/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                       GCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354144‑354437

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