Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,457,343	G→C	5.9%	R46G (CGC→GGC)	PP_1275 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,457,343	0	G	C	5.9%	210.7 / 3.0	68	R46G (CGC→GGC)	PP_1275	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base C (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGTAATCGCCT  >  NC_002947/1457196‑1457491                                                                                                                                                    |                                                                                                                                                      gTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCgg                                                                                                                                                      <  8:399739/149‑1 (MQ=255)      tACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCt                                                                                                                                                    <  7:300852/146‑1 (MQ=255)      tACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCt                                                                                                                                                    >  8:300852/1‑146 (MQ=255)      tACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTc                                                                                                                                                 <  1:274293/149‑1 (MQ=255)      tACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTc                                                                                                                                                 >  6:319739/1‑149 (MQ=255)          ccgccACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGc                                                                                                                                             <  5:319739/149‑1 (MQ=255)           cgccACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGc                                                                                                                                             >  3:428044/1‑148 (MQ=255)            gccACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCa                                                                                                                                           <  1:286271/149‑1 (MQ=255)              cACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGa                                                                                                                                         <  6:397932/149‑1 (MQ=255)              cACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGa                                                                                                                                         >  7:78900/1‑149 (MQ=255)               aCGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGAc                                                                                                                                        <  4:411208/149‑1 (MQ=255)               aCGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGAc                                                                                                                                        <  4:417757/149‑1 (MQ=255)               aCGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGAc                                                                                                                                        <  2:178304/149‑1 (MQ=255)                  cgtcAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACGTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCAc                                                                                                                                     <  6:33259/146‑1 (MQ=255)                    ccAGCGCAGTTCTGGCAGCGGCACGGCGTTGCCGTTCTTCATGCGTACCTGCAAGGCTCGGTCGGCGGGACAGTTTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGAGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGa                                                                                                                                   <  8:274755/149‑1 (MQ=255)                         gcAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAg                                                                                                                              <  3:109938/149‑1 (MQ=255)                          cAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAgg                                                                                                                             <  2:407253/149‑1 (MQ=255)                           aGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGc                                                                                                                            <  2:212009/149‑1 (MQ=255)                               cGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGccac                                                                                                                         <  7:93331/148‑1 (MQ=255)                                  tCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCAccag‑g                                                                                                                    >  8:112019/1‑149 (MQ=255)                                        gCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCAccag‑gcccagg                                                                                                              <  7:399493/149‑1 (MQ=255)                                           cgtcgtTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCca                                                                                                           <  1:74296/149‑1 (MQ=255)                                             tcgtTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCcaca                                                                                                         <  6:384322/149‑1 (MQ=255)                                                   ccGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAgc                                                                                                   >  5:385212/1‑149 (MQ=255)                                                    cGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAgcc                                                                                                  >  1:281168/1‑149 (MQ=255)                                                             tGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCAGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAgg                                                                                            >  1:298132/1‑146 (MQ=255)                                                               cGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGAt                                                                                       <  4:245745/149‑1 (MQ=255)                                                                 tACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATgg                                                                                     <  2:50571/149‑1 (MQ=255)                                                                          ggCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCt                                                                            >  3:163721/1‑149 (MQ=255)                                                                          ggCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCt                                                                            >  7:309346/1‑149 (MQ=255)                                                                               ggTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAGCTGGAAAGCCACACACAGCCGCCAGGGAGATGGGAAGATCCTTGGc                                                                        >  1:364183/1‑149 (MQ=255)                                                                                 tCGGGGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCGAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTGTTGCCAGACGACGAACAAGGCCACCAG‑GCGCATGCCACACACTGCCGCCAGGGAGATGGGAAGAGCCTTTGCTgg                                                                     <  4:83523/149‑1 (MQ=255)                                                                                 tCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTgg                                                                     >  5:216307/1‑149 (MQ=255)                                                                                    gcggGACAGTGTTCGGGGGCGTAGGTCACACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGt                                                                   >  1:247069/1‑148 (MQ=255)                                                                                        gACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACg                                                               >  5:405384/1‑148 (MQ=255)                                                                                           aGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGcaa                                                            >  3:321330/1‑148 (MQ=255)                                                                                               ttCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGcaacaac                                                        <  8:78900/148‑1 (MQ=255)                                                                                                  gggggCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCaacaacaac                                                     <  1:427512/148‑1 (MQ=255)                                                                                                    gggCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTaa                                                                 <  8:159353/134‑1 (MQ=255)                                                                                                    gggCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTaa                                                                 >  7:159353/1‑134 (MQ=255)                                                                                                    gggCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGAAACAACAACAcc                                                  >  5:11582/1‑149 (MQ=255)                                                                                                     ggCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACAcca                                                 <  6:216307/149‑1 (MQ=255)                                                                                                       cGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACAccagc                                               <  2:298132/149‑1 (MQ=255)                                                                                                          aGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACACCAccagccag                                            >  8:241953/1‑149 (MQ=255)                                                                                                          aGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCCGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGcaacaaca                                                       <  8:51997/138‑1 (MQ=255)                                                                                                          aGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCCGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGcaacaaca                                                       <  8:45379/138‑1 (MQ=255)                                                                                                          aGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCCGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGcaacaaca                                                       >  7:51997/1‑138 (MQ=255)                                                                                                          aGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCCGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGcaacaaca                                                       >  7:45379/1‑138 (MQ=255)                                                                                                            gTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGt                                          >  4:176787/1‑149 (MQ=255)                                                                                                               aGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGa                                        <  2:281168/148‑1 (MQ=255)                                                                                                               aGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAg                                       >  1:231313/1‑149 (MQ=255)                                                                                                                     aGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATcag                                 >  6:212503/1‑149 (MQ=255)                                                                                                                       gTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATcag                                 >  5:304755/1‑147 (MQ=255)                                                                                                                        tCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATcagcag                              >  2:64646/1‑149 (MQ=255)                                                                                                                          cAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAATGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAggg                            <  3:245595/149‑1 (MQ=255)                                                                                                                          cAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCGACACACAGCCGCCGGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAggg                            <  6:385212/149‑1 (MQ=255)                                                                                                                            ggCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCa                          <  4:163721/149‑1 (MQ=255)                                                                                                                             gcgcgcCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCAc                         <  4:86698/149‑1 (MQ=255)                                                                                                                                cgcCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACCGCCGCCAGGGAGATGGGAAGAGCCTTGGATGGGGAACGcaac                                                           >  3:184096/1‑112 (MQ=255)                                                                                                                                 gcCAGTCGAATGCTTTCGCGGCTGTCTTGCAAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCAccta                      >  4:6201/1‑146 (MQ=255)                                                                                                                                  ccAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATc                    <  2:247069/149‑1 (MQ=255)                                                                                                                                    aGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAg                  <  5:212503/149‑1 (MQ=255)                                                                                                                                    aGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCT‑GGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGTCACCGATCAGc                 <  8:309346/149‑1 (MQ=255)                                                                                                                                           tGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATcgcg           >  7:187432/1‑149 (MQ=255)                                                                                                                                             cTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGt          >  3:381842/1‑148 (MQ=255)                                                                                                                                             cTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGt          >  6:98395/1‑148 (MQ=255)                                                                                                                                              tttCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGa                                        >  3:97725/1‑117 (MQ=255)                                                                                                                                              tttCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGa                                        <  4:97725/117‑1 (MQ=255)                                                                                                                                              tttCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGTaa        >  6:419294/1‑149 (MQ=255)                                                                                                                                                 cgcgGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGTAATCg     >  2:235840/1‑149 (MQ=255)                                                                                                                                                  gcgGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGTAATCGc    >  7:3302/1‑149 (MQ=255)                                                                                                                                                    ggCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGTAATCGCCt  >  1:355887/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGTTCGGGGGCGTAGGTCAGACGCAGGTCCAGGCGCGCCAGTCGACTGCTTTCGCGGCTGTCTTGCCAGACCACGAACAAGGCCACCAG‑GCCCAGGCCACACACAGCCGCCAGGGAGATGGGAAGAGCCTTGGCTGGGTAACGCAACAACAACACCAGCCAGGTGAGGATCAGCAGGGCACCGATCAGCATCGCGTAATCGCCT  >  NC_002947/1457196‑1457491

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