Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,282	0	T	G	42.5%	124.2 / 14.2	75	V526G (GTT→GGT)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (12/30);  new base G (31/0);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 4.10e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTGACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCC  >  NC_000913/354132‑354431                                                                                                                                                       |                                                                                                                                                      ccTTTGACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCGGCCATCCGGGCGTTGCCGAAGTGGCggtggt                                                                                                                                                       <  2:1234741/151‑1 (MQ=255)     tGACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTgg                                                                                                                                                    <  2:503463/150‑1 (MQ=255)      gACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTggg                                                                                                                                                   <  2:104177/150‑1 (MQ=255)          ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTg                                                                                                                                                <  2:295812/149‑1 (MQ=255)           gggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGaaa                                                                                                                                             <  2:486353/151‑1 (MQ=255)              cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                                          <  2:575356/151‑1 (MQ=255)              cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTgggggg                                                                                                                                                >  1:443185/1‑143 (MQ=255)               aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATg                                                                                                                                         <  1:244715/151‑1 (MQ=255)                     gACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTg                                                                                                                                    <  2:606732/150‑1 (MQ=255)                     gACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAt                                                                                                                                          >  1:2836/1‑144 (MQ=255)                        gCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGaa                                                                                                                                  >  1:242147/1‑149 (MQ=255)                           gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAgg                                                                                                                               <  2:484826/149‑1 (MQ=255)                            aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCTGGAAggg                                                                                                                              >  2:1233733/1‑149 (MQ=255)                            aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGTAAAAGGCGCTGGAAggg                                                                                                                              >  2:506417/1‑149 (MQ=255)                             cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAg                                                                                                                           <  1:490601/151‑1 (MQ=255)                             cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAAGGGc                                                                                                                             >  1:168182/1‑149 (MQ=255)                              ggTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAgg                                                                                                                          <  2:1145529/151‑1 (MQ=255)                               gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGGGCTGGAAGGGGAgg                                                                                                                          >  1:1037140/1‑150 (MQ=255)                                ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAAGGCAgggg                                                                                                                        >  1:1082664/1‑149 (MQ=255)                                  aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtgg                                                                                                                       <  1:724193/150‑1 (MQ=255)                                  aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtgg                                                                                                                       <  2:967827/150‑1 (MQ=255)                                  aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGGGCTGGAAGGGGAggggg                                                                                                                       >  2:807015/1‑150 (MQ=255)                                   tCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGGAAGGGGAggtgggg                                                                                                                     >  2:1207715/1‑149 (MQ=255)                                     aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggt                                                                                                                   <  1:591545/151‑1 (MQ=255)                                               cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGttt                                                                                                            <  1:384335/148‑1 (MQ=255)                                               cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTg                                                                                                           <  1:1034505/149‑1 (MQ=255)                                               cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGt                                                                                                          >  2:769178/1‑150 (MQ=255)                                                gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCa                                                                                                        <  2:1233835/151‑1 (MQ=255)                                                  gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCAtt                                                                                                      <  1:1212359/151‑1 (MQ=255)                                                      aCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTTATTCCg                                                                                                   >  1:935022/1‑150 (MQ=255)                                                      aCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTcc                                                                                                    >  1:88290/1‑149 (MQ=255)                                                         gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGGGGGGGGAAAGACGGGCTCGAAGGGTAGGAGGTGGCGGCGGTTTTTATTcc                                                                                                    >  1:817867/1‑146 (MQ=255)                                                            gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagaga                                                                                            <  1:528460/151‑1 (MQ=255)                                                            gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagaga                                                                                            <  2:1024202/151‑1 (MQ=255)                                                            gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAAGGCAGGGGGCGGGGGGGTTTGTCATTCCGaaa                                                                                                >  1:354788/1‑147 (MQ=255)                                                             aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAAAAGCGCTGTAAGGGGAGGTGGCGGGGGGGGTTGTGAtt                                                                                                      >  1:683951/1‑140 (MQ=255)                                                               gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGGGAGTATCTCCAGCCATCCGGGGGTTGCCGAAGTGGCGGTGG‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTcc                                                                                                    <  1:842129/139‑1 (MQ=255)                                                                tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGCTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGGGGCGGTGGATTGGGGGAAAAGGGTGCCGTAAGGGCCGGTGGCGGTTGCGCTTTTCAATACGAAagggaatg                                                                                         >  1:582882/1‑147 (MQ=255)                                                                  aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCATTCCGAAAGAGAGCGAc                                                                                       >  1:111092/1‑150 (MQ=255)                                                                  aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCg                                                                                         >  1:571695/1‑148 (MQ=255)                                                                     aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGTCGGGGTGAAAGAAGAGCGTCAAGGGGTGGGGGGGGGGGCGTGTGTCCTTCCGaaaaaaagaggcag                                                                                     >  1:156720/1‑138 (MQ=255)                                                                        gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGGGCCGGAAGGGCAGGGGGGGGTGGGGGTTGTTATTTCGAAAGAGAGGGACAGTCTgg                                                                                >  1:1230669/1‑151 (MQ=255)                                                                        gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTTCGAAAGAGAGCGACAGTCTgg                                                                                >  1:200259/1‑151 (MQ=255)                                                                         ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCATTCCGCAAGAGAGGGACAGGCTggg                                                                               >  2:476107/1‑150 (MQ=255)                                                                             cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                           <  2:1240217/151‑1 (MQ=255)                                                                              gggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGGTATTCCGAAAGAGAGCGACAGTCTgg                                                                                >  2:904125/1‑145 (MQ=255)                                                                               ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGGGCTGAAAGGGGAGGGGGCGGGGGCGGTTGTCATTCCGGAAGAGAGCGGCAGGCTgg                                                                                >  2:176796/1‑144 (MQ=255)                                                                                 cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGCGGGGGCGGTTTTCATTCCGAAAGAGAGCGACAGTCTGGAAGAACGTg                                                                       >  1:1127661/1‑151 (MQ=255)                                                                                           ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGAGGGGCTGGAAGGGGAGGTGGCGGGGGGGGTTTTCATTCCGGAAGAGAGCGGCTGGCTGGGAGGCCGTGgtgttggg                                                               >  2:256521/1‑145 (MQ=255)                                                                                               gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAAGGCAGGGGGCGGTGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCCCTc                                                          >  1:219171/1‑150 (MQ=255)                                                                                                  cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCaa                                                      <  1:531051/151‑1 (MQ=255)                                                                                                    agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAaga                                                    <  1:1072781/151‑1 (MQ=255)                                                                                                         gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAgg                                               >  1:1111421/1‑151 (MQ=255)                                                                                                           agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCg                                             <  1:1104590/151‑1 (MQ=255)                                                                                                             agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAt                                           >  1:982045/1‑151 (MQ=255)                                                                                                              gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAtt                                          <  1:476109/151‑1 (MQ=255)                                                                                                               agTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTa                                         <  1:732338/151‑1 (MQ=255)                                                                                                                gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTAt                                        <  2:2836/151‑1 (MQ=255)                                                                                                                 tATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGAGGCGAAAGGGGAGGTGGCGGTGGGGTTTTTTATTCCGCAAGAGAGCGGGAGGGTTGGAAATCGTGATGGGGCGCACCCGCAAGAGAAGGGGAATa                                         >  2:813143/1‑149 (MQ=255)                                                                                                                       cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCCGGAAGGGGAGGGGGCGGGGGCGGTTGTTATTTCGAAAGAGAGCGACAGGCTGGGAGAACGTGGTGTGGGGCCCTCGCAAGAGAAGGCGGTTATGGCGGtg                                 >  2:518248/1‑151 (MQ=255)                                                                                                                         gCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggt                               <  1:318155/151‑1 (MQ=255)                                                                                                                          ccATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAGGGGAGGGGGCGGGGGGGGTTGTCATTCCGGAAGAGAGCGACAGTCTGGAAGATCGTGGTGTGGGGCACTCGCCAGAGAAGGCGGTTATGGCGGtgggg                              >  2:252694/1‑149 (MQ=255)                                                                                                                            aTCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAAATGCGCTGAAAGGGCAGGTGGGGGGGGCGTTTGTCCTTCCCAAAAAGAGCGACAGTCTGGAAGATCGTGGTGTGGGGCCCTCGCAAGAGAAGGTGGATATTGGGCcggtgg                             >  1:1023997/1‑150 (MQ=255)                                                                                                                            aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCCCTGAAAGGGCAGGGGGCGGGGGGGTTTGTCATTCCGAAAGAGGGCGGCAGTCTGGGAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGGTTATTGCGCtgg                                >  1:608520/1‑147 (MQ=255)                                                                                                                            aTCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGCGCTGAAAGGGCAGGTGGGGGGGGGGGTTTTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGGTGTGGGGCCCCCGCCAAAGAAGGGGGTTATGGCGCtggggg                             >  1:566212/1‑150 (MQ=255)                                                                                                                               cGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAg                         <  1:363621/151‑1 (MQ=255)                                                                                                                               cGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGGAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGGCAGTCTGGGAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAGTATGGCGCTGGGGGAc                           >  2:429429/1‑149 (MQ=255)                                                                                                                                gggCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGGTGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgggct                          >  1:1111940/1‑146 (MQ=255)                                                                                                                                gggCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCCTTTCGGAAGAGAGGGGCAGTCTGGGAGAACGGGGTGGGGGGCCCTCCCAAGAGAAGGGGGGTATTGGGCtgggggg                            >  2:15892/1‑146 (MQ=255)                                                                                                                                         cGAAGTGGCGGTCGGTGGGGTGGAAGATGCGCTCCAAGGGCAGGGGGGGGGGGCGTTTGGTATTCCGCAAGAGGGGGACAGTCTGGAAGATCGTGATGGGGCGCACCCGCAAGAGAAGGCGATTATTGCGCCTGTGGGCCGCCAGATTgg                >  2:1155141/1‑150 (MQ=255)                                                                                                                                          gAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGTGGGGGGGGGGGTTGTCATTCCGAAAGAGAGGGACAGGCTGGGAGAACGTGATGTGGCGCACTCGCAAGAGAAGGGGGTTATGGGGCTGGGGGACAGCCAGAtt                  >  2:477796/1‑147 (MQ=255)                                                                                                                                             gTGGCGGTGGTTGGGGGGGAAAAAGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGGGACAGTCTGGGAGATCGGGGTGTGGGGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGGCAGCCAGATTGGCaa             >  2:254758/1‑149 (MQ=255)                                                                                                                                               ggCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt         <  2:315491/151‑1 (MQ=255)                                                                                                                                               ggCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt         >  1:1098534/1‑151 (MQ=255)                                                                                                                                               ggCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt         <  1:755978/151‑1 (MQ=255)                                                                                                                                                 cGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTgg       >  2:416852/1‑151 (MQ=255)                                                                                                                                                   gtggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc     <  2:582425/151‑1 (MQ=255)                                                                                                                                                     ggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgc   >  1:790640/1‑151 (MQ=255)                                                                                                                                                      gtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgcc  <  2:443183/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CCTTTGACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCC  >  NC_000913/354132‑354431

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