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

GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑354437                                                                                                                                                      |                                                                                                                                                       ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTg                                                                                                                                                       <  1:492694/149‑1 (MQ=255)            gACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGATGC‑GCt                                                                                                                                           >  1:617954/1‑149 (MQ=255)              cGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGaaa                                                                                                                                       <  1:45820/151‑1 (MQ=255)                      gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAggt                                                                                                                               <  2:32081/151‑1 (MQ=255)                        tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAggtg                                                                                                                              <  2:452354/150‑1 (MQ=255)                          tCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGGGGAAGATGC‑GCCGGAAGGGGAgggggcg                                                                                                                           >  2:97319/1‑151 (MQ=255)                             cTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGGGGAAGAAGC‑GCTGGAAGGGCGgggggggggg                                                                                                                        >  2:717611/1‑149 (MQ=255)                              tttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGGAAAATGCGGCGGAAAgggaggggggggcg                                                                                                                         >  2:553394/1‑139 (MQ=255)                               ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGGAAGATGC‑GGTGGAAGGGGAggtggggggggg                                                                                                                      >  2:269315/1‑150 (MQ=255)                                   tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGGAAGGTGG‑GGTGGAAGGGGAggtgggggggggggt                                                                                                                   >  2:142469/1‑146 (MQ=255)                                      cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc                                                                                                               <  2:961769/151‑1 (MQ=255)                                      cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGATGC‑GCTGGAAGGGCAGGTGGCGGGGGCGGTTGt                                                                                                                >  1:534779/1‑150 (MQ=255)                                         gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCAtt                                                                                                            <  2:696296/151‑1 (MQ=255)                                          cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGGAAGATGC‑GCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCCtt                                                                                                            >  1:1007911/1‑150 (MQ=255)                                          cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑G‑GGGTGGGGGGGAAGATGC‑GCGGGAAGGGGAGGGGGCGGGGGGGGTTggtattt                                                                                                           >  2:26732/1‑145 (MQ=255)                                               tGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑A‑GGGTGGGGGGGAAAAGGC‑GCGGGAAAGGGAGGGGGGGGGGGGGTTTGTCCTTcccga                                                                                                       >  2:799433/1‑147 (MQ=255)                                                  cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                                                                                   <  1:386176/151‑1 (MQ=255)                                                   gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGATGC‑GCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCATTCCGaaa                                                                                                      >  2:331354/1‑147 (MQ=255)                                                         aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGGAAAGAAG‑‑GCCGGAAGGGGCGGGGGGGGCGGCGGTTGGCATTTCGCAAGAGGGCg                                                                                               >  2:95566/1‑148 (MQ=255)                                                           tAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                                                                                           <  1:927656/150‑1 (MQ=255)                                                             aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCt                                                                                        <  2:604394/151‑1 (MQ=255)                                                               gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑G‑GGGTGGGGGGGAAAATGC‑GCTGGAAGGGGAGGGGGCGGGGGGGGTTGTCTTTTCGCAAAAGAGGGACAGGCTgg                                                                                      >  2:691876/1‑151 (MQ=255)                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGaa                                                                                    <  1:793232/150‑1 (MQ=255)                                                                    cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGAAGG‑GCTGGAAGGGGCGGGGGGGGGGGGGGTTGTTATTTCGGAAGAGAGCGACAGGCTggga                                                                                    >  2:727484/1‑146 (MQ=255)                                                                    cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑C‑GGTTGGGGTGAAAGATGG‑GCTGCAAGGGGCGGGGGCGGGGGGGGTTGCCGTTTCGAAAAAGCGCGACAGtggggga                                                                                    >  1:520164/1‑142 (MQ=255)                                                                     gggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCG‑GGTGGGGGGGAAAGA‑GC‑GCTGAAAGGGCCGGGGGGGGGGGCGGTTGTCATTCCGCAAGAAAGAGACAGACTGGaa                                                                                    >  2:944108/1‑147 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGa                                                                            <  1:163464/151‑1 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGGAAAATGC‑GCTGAAAAGGGAGGGGGCGGTGGGGTTTGTCCTTCCCAAAGAGAGCGACAGTCTGGAAGATCGTg                                                                             >  1:207497/1‑150 (MQ=255)                                                                          tCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGGGAAAGATGC‑GCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTg                                                                             >  1:288997/1‑149 (MQ=255)                                                                            ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAtgt                                                                         <  1:725381/151‑1 (MQ=255)                                                                            ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGGGGAAGATGC‑GCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTTCGGAAGAGAGCGACAGGCTGGAAGATCGTg                                                                             >  1:429635/1‑147 (MQ=255)                                                                                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTgg                                                                       <  1:507043/149‑1 (MQ=255)                                                                                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgc                                                                      <  2:917575/150‑1 (MQ=255)                                                                                    aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACt                                                                 <  1:788833/151‑1 (MQ=255)                                                                                    aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGGAAGATGC‑GCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCATTTCGGAAGAGAGGGACAGGCTTGAAGATCGTGGTGGGGGGCAc                                                                  >  1:592545/1‑150 (MQ=255)                                                                                     cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                                <  1:507398/151‑1 (MQ=255)                                                                                       cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                              <  2:686551/151‑1 (MQ=255)                                                                                       cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                              <  1:133689/151‑1 (MQ=255)                                                                                       cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGAAAGAGAGGGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                                >  1:961103/1‑149 (MQ=255)                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCaa                                                            <  2:127929/151‑1 (MQ=255)                                                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAaga                                                          <  1:667773/151‑1 (MQ=255)                                                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGGTGC‑CCCGAAAAGGGCGGGGGCGGGGGGGGTTGTCATTCCGGAAGGGAGCGACACTCTGGAAGATGGTGGAAGGGCGCGCTCGc                                                              >  1:782037/1‑147 (MQ=255)                                                                                            gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAagag                                                         <  1:739119/151‑1 (MQ=255)                                                                                               tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGaa                                                       <  2:42205/150‑1 (MQ=255)                                                                                                gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAgg                                                     <  1:535676/151‑1 (MQ=255)                                                                                                   gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑A‑GGGTGGGGGGAGAAATGC‑GCTGGAAGGGGCGGGGGGGGGGGCGGTTGGGAGTTCGGAAGAGAGCGAGAGGCTTGGAGATCGGGATGTTGCGGACTCGCCAGGGAAgg                                                     >  2:314758/1‑148 (MQ=255)                                                                                                    agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGATGC‑GCTGAAAGGGGAGGGGGCGGGGGCGGTTGTCATTTCGAAAGAGAGCGACAGTCTGGGAGATCGTGGTGGGGCGCACTCGCAAGAGAAgg                                                     >  2:740539/1‑147 (MQ=255)                                                                                                          tctcCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGGGAAAGATGC‑GCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCCTTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGGGCACTCGCAAGAGAAGGCggatttgg                                            >  1:834454/1‑143 (MQ=255)                                                                                                                   aTCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTTCGAAAGAGAGCGACAGTCTTGAAGATCGTGATGTGGGGCACTCGCAAGAGAAGGCGGTTATTGCGCtggtgg                                   >  1:665646/1‑150 (MQ=255)                                                                                                                   aTCCGGGCGTTGCCGAAGTGGCGG‑T‑GGGTGGGGGGAAAAATGC‑GCTGCAAGGGGAGGGGGCGGGGGGGGTTTTCGTTTCGAAAGAGAGGGGCAGGCTTGAAGATCGTGGTGTGGGGCGCTCGCAAAAGAAGGCGGTTATGGCGCcggtgg                                   >  2:499355/1‑150 (MQ=255)                                                                                                                     ccGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACa                                <  1:172665/151‑1 (MQ=255)                                                                                                                      cGGGCGTTGCCGAAGTGGCGG‑A‑GGGTGGGGGGAAAGATGC‑GCTGAAAGGGGAGGGGGGGGGGGCGTTTGTCATTCCCAAAGAGAGCGACAGTCTGGAAGATCGTGATGGGGGGCACCCGCAAGAGAAGGGGGTTATGGCGCtggggg                                   >  2:639554/1‑147 (MQ=255)                                                                                                                              gCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTg                       <  1:441513/151‑1 (MQ=255)                                                                                                                                       gCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg              >  1:896324/1‑151 (MQ=255)                                                                                                                                        cGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTgg             >  2:630539/1‑151 (MQ=255)                                                                                                                                          g‑t‑ggtTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc           <  2:643139/151‑1 (MQ=255)                                                                                                                                               gtTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgcc        >  2:53117/1‑151 (MQ=255)                                                                                                                                                     gTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  >  2:406239/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑T‑GGTTGGGGTGAAAGATGC‑GCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑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: ‑