Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,529,153	T→C	5.8%	L59P (CTG→CCG)	ftsZ →	cell division protein FtsZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,529,153	0	T	C	5.8%	202.2 / 2.2	69	L59P (CTG→CCG)	ftsZ	cell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base T (36/29);  new base C (2/2);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGG  >  NC_002947/1529009‑1529300                                                                                                                                                 |                                                                                                                                                    gTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCaa                                                                                                                                                 <  6:198059/149‑1 (MQ=255)    cGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTCGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTg                                                                                                                                              >  8:94828/1‑149 (MQ=255)      gTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTggg                                                                                                                                            >  6:122252/1‑149 (MQ=255)      cgatgtAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTgg                                                                                                                                             <  7:71844/143‑1 (MQ=255)         aTTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCAc                                                                                                                                         >  7:113049/1‑149 (MQ=255)            aaaGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACgg                                                                                                                                       <  6:278736/148‑1 (MQ=255)             aaGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTggg                                                                                                                                            <  4:101113/142‑1 (MQ=255)             aaGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTggg                                                                                                                                            >  3:101113/1‑142 (MQ=255)             aaGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGAAAtttg                                                                                                                                             >  8:71844/1‑139 (MQ=255)              aGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGgtg                                                                                                                                    >  4:106024/1‑149 (MQ=255)               gTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGTGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGgtgt                                                                                                                                   <  3:71771/149‑1 (MQ=255)               gTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTTGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGgtgt                                                                                                                                   <  8:234394/149‑1 (MQ=255)                  aTCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGAc                                                                                                                                >  8:4811/1‑149 (MQ=255)                                  gtggCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgcc                                                                                                                >  2:247475/1‑149 (MQ=255)                                   tggCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccg                                                                                                               <  6:330028/149‑1 (MQ=255)                                   tggCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccg                                                                                                               <  7:4811/149‑1 (MQ=255)                                     gcggcAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccggt                                                                                                             <  7:249184/149‑1 (MQ=255)                                      cggcAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACTAAGGGCCTGggtgccggtg                                                                                                            <  6:189814/149‑1 (MQ=255)                                       ggcAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccggtgc                                                                                                           <  3:291871/149‑1 (MQ=255)                                        gcAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGggtgccggtgc                                                                                                           >  4:289351/1‑148 (MQ=255)                                          aaCGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCaa                                                                                                        >  1:209351/1‑149 (MQ=255)                                           aCGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAAt                                                                                                       >  1:31542/1‑149 (MQ=255)                                                gTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCgg                                                                                                   >  1:149463/1‑148 (MQ=255)                                                gTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGa                                                                                                  <  5:122252/149‑1 (MQ=255)                                                 tCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGaa                                                                                                 <  2:31542/149‑1 (MQ=255)                                                   aaTCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCACAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGTGCCTGGGTGCCGGTGCCAATCCGGAAGt                                                                                               >  1:255218/1‑149 (MQ=255)                                                        caTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTc                                                                                              <  3:233594/145‑1 (MQ=255)                                                        caTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTc                                                                                              >  4:233594/1‑145 (MQ=255)                                                        caTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGcc                                                                                          <  3:289351/149‑1 (MQ=255)                                                        caTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGcc                                                                                          <  6:325531/149‑1 (MQ=255)                                                        caTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGcc                                                                                          >  4:224171/1‑149 (MQ=255)                                                         aTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCCGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATcc                                                                                                     <  5:99585/137‑1 (MQ=255)                                                         aTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCCGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATcc                                                                                                     >  6:99585/1‑137 (MQ=255)                                                         aTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCCGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATcc                                                                                                     >  6:104170/1‑137 (MQ=255)                                                         aTGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCCGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATcc                                                                                                     <  5:104170/137‑1 (MQ=255)                                                             tCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAg                                                                                     >  7:152842/1‑149 (MQ=255)                                                               aaGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAgg                                                                                    <  2:209351/148‑1 (MQ=255)                                                                 gagCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCcg                                                                                 >  4:52645/1‑149 (MQ=255)                                                                 gagCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCAAGGCGCTGAAAAACATTGGCGCGCGCACCATACTGCAATTGGGCACGGGTGTGACAAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCcg                                                                                 >  8:106635/1‑149 (MQ=255)                                                                  agcagcATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCcgc                                                                                <  2:244782/149‑1 (MQ=255)                                                                    cagcaTCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCgcgc                                                                              >  6:250671/1‑149 (MQ=255)                                                                          cGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTgg                                                                           >  8:122707/1‑146 (MQ=255)                                                                           gatgtcGTGCAAGTCTTCTGCGCCATCACCGATGCCCAGGCGCTGAAATCCTTTGGCGCGCGAAGCATCCTGCAATTGGGCTCGGGTGTGAGCAAGGGCCTGGGTGGCTGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAg                                                                        <  7:283110/143‑1 (MQ=255)                                                                                 gTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGt                                                                   >  2:64100/1‑147 (MQ=255)                                                                                 gTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGa                                                                 <  6:35598/149‑1 (MQ=255)                                                                                      aTTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCa                                                            >  2:196521/1‑149 (MQ=255)                                                                                              gcgcCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCg                                                      >  2:293106/1‑147 (MQ=255)                                                                                                  cAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGtgc                                                >  7:73387/1‑149 (MQ=255)                                                                                                    acacCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGtgctg                                              <  7:122707/149‑1 (MQ=255)                                                                                                       ccGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAg                                           <  5:260726/149‑1 (MQ=255)                                                                                                        cGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAgg                                          >  3:242902/1‑149 (MQ=255)                                                                                                        cGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAgg                                          >  3:134271/1‑149 (MQ=255)                                                                                                          aTGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAggg                                         >  3:251687/1‑148 (MQ=255)                                                                                                          aTGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGc                                        <  3:52645/149‑1 (MQ=255)                                                                                                             cccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGGCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCAc                                      <  1:247475/148‑1 (MQ=255)                                                                                                             cccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCAcc                                     <  3:224171/149‑1 (MQ=255)                                                                                                              ccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGtgct                                               >  7:37754/1‑138 (MQ=255)                                                                                                              ccAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGtgct                                               <  8:37754/138‑1 (MQ=255)                                                                                                                aGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCaa                                   >  8:119325/1‑148 (MQ=255)                                                                                                                 ggCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACa                                 >  2:279/1‑149 (MQ=255)                                                                                                                   cgcTGAAAAAGATTGGCGCGCGCACGATCCTGCAATTGGGCACGGGTGTTACCAAGGGCCTGGTTGCCGGTTCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGATCGCATTGCAGAAGTTCTGCAGGGCACCAACATg                               <  4:134271/149‑1 (MQ=255)                                                                                                                      tGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGTCGAAGTGCTGCAGGGCACCAACATGgtg                            >  2:306564/1‑149 (MQ=255)                                                                                                                        aaaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGgtgt                           <  1:279/148‑1 (MQ=255)                                                                                                                         aaaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGtt                          >  6:110034/1‑148 (MQ=255)                                                                                                                          aaaCATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCAGAAGTGCTGCAGGGCACCAACATGGTGTtca                        >  8:81468/1‑149 (MQ=255)                                                                                                                            aCATTGGCGCGCGCACCATCCTGCAATTGCGCACGGGTGTGACAAAGGGCCAGGGTTCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCAGCACGGAACCAACATGGTGTtcatc                      >  1:132591/1‑149 (MQ=255)                                                                                                                                         cACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGTCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATgg         <  1:306564/149‑1 (MQ=255)                                                                                                                                            cATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTg      >  4:190227/1‑149 (MQ=255)                                                                                                                                            cATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCAGTCAGGCCGCGCTGGAAGAGCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTg      >  7:132426/1‑149 (MQ=255)                                                                                                                                                 tGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGcgg  >  6:81916/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GTCCGGTCATTAAAGTGATCGGCGTTGGTGGTGGTGGCGGCAACGCCGTCAATCACATGGTCAAGAGCAGCATCGAAGGCGTGGAATTCATCTGCGCCAACACCGATGCCCAGGCGCTGAAAAACATTGGCGCGCGCACCATCCTGCAATTGGGCACGGGTGTGACCAAGGGCCTGGGTGCCGGTGCCAATCCGGAAGTCGGCCGTCAGGCCGCGCTGGAAGACCGTGAGCGCATTGCCGAAGTGCTGCAGGGCACCAACATGGTGTTCATCACCACCGGCATGGGTGGCGG  >  NC_002947/1529009‑1529300

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