Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,109,140	T→A	5.4%	D38V (GAC→GTC)	PP_0968 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,109,140	0	T	A	5.4%	220.2 / 2.6	74	D38V (GAC→GTC)	PP_0968	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

TTTACCATTGATCCCCCCGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTG  >  NC_002947/1108992‑1109288                                                                                                                                                     |                                                                                                                                                     tttACCATTGATCCCCCCGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGt                                                                                                                                                      >  5:336099/1‑149 (MQ=255)          gATCCCCCCGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCACGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCg                                                                                                                                             >  5:223098/1‑149 (MQ=255)            tCCCCCCGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCTGTGACccggcccggc                                                                                                                                           <  7:210790/149‑1 (MQ=255)              cccccGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTg                                                                                                                                         <  5:153098/149‑1 (MQ=255)               ccccGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTg                                                                                                                                         >  2:242598/1‑148 (MQ=255)                 ccGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATGGCCCAATTGCAGTTACTGCTGCAATTTCACGCCTTCAAGTCCGATGAGCAAGTTTTTATGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGctg                                                                                                                                      <  4:50345/149‑3 (MQ=255)                    gCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGc                                                                                                                                   <  2:26268/149‑1 (MQ=255)                       tAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAg                                                                                                                                >  1:162631/1‑149 (MQ=255)                       tAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAg                                                                                                                                <  2:140822/149‑1 (MQ=255)                       tAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAg                                                                                                                                <  2:162631/149‑1 (MQ=255)                        aaGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAg                                                                                                                                >  2:90659/1‑148 (MQ=255)                          gCCTTGCCGTAGTGCAGCCGCTGCTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGaa                                                                                                                             <  5:23935/149‑1 (MQ=255)                           ccTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGaa                                                                                                                             >  7:21560/1‑148 (MQ=255)                             gtGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAg                                                                                                                          <  8:74176/148‑1 (MQ=255)                               gCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGc                                                                                                                         <  3:33272/148‑1 (MQ=255)                               gCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGc                                                                                                                         <  3:37854/148‑1 (MQ=255)                                     gTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACa                                                                                                                  <  1:90659/149‑1 (MQ=255)                                     gTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCGCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACa                                                                                                                  <  2:341123/149‑1 (MQ=255)                                     gTGCAGCCACTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATaa                                                                                                                    >  6:94804/1‑147 (MQ=255)                                               tGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGGAAGGTTTTAGTGCGGCTCGGGTGACGGGCCGGGTCACAGCCCGGCTGCAGGGCGAGGAACAGCAATATAAGCAAACATga                                                                                                        >  8:146040/1‑148 (MQ=255)                                                 tttCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCaaa                                                                                                             >  7:222564/1‑142 (MQ=255)                                                 tttCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCaaa                                                                                                             <  8:222564/142‑1 (MQ=255)                                                 tttCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGc                                                                                                        >  3:263574/1‑147 (MQ=255)                                                   tCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAg                                                                                                     <  6:336099/148‑1 (MQ=255)                                                     aGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCg                                                                                                  <  4:68389/149‑1 (MQ=255)                                                       gTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCa                                                                                                >  6:311140/1‑149 (MQ=255)                                                        ttGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCAt                                                                                               <  5:311140/149‑1 (MQ=255)                                                           gTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACACCAAACTTGCCAGTCGCATcgg                                                                                            >  1:239617/1‑149 (MQ=255)                                                              gCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTc                                                                                         <  6:335535/149‑1 (MQ=255)                                                               cATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCg                                                                                        >  7:276977/1‑149 (MQ=255)                                                                 tCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCa                                                                                      <  7:123386/149‑1 (MQ=255)                                                                    cccAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATcc                                                                                   >  1:327445/1‑149 (MQ=255)                                                                       aaTTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCtt                                                                                 >  2:90626/1‑148 (MQ=255)                                                                           gCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCTGTGACccggcccggcac                                                                                                                                         >  7:352008/1‑86 (MQ=38)                                                                           gCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCTGTGACccggcccggcac                                                                                                                                         <  8:352008/88‑3 (MQ=38)                                                                            cAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAg                                                                            >  5:109398/1‑148 (MQ=255)                                                                            cAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAg                                                                            >  5:166707/1‑148 (MQ=255)                                                                            cAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGc                                                                           <  7:339150/149‑1 (MQ=255)                                                                                    gctgcAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAAc                                                                    >  3:239790/1‑148 (MQ=255)                                                                                       gcAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGc                                                                >  8:16402/1‑149 (MQ=255)                                                                                          aTTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGtt                                                             <  3:246225/149‑1 (MQ=255)                                                                                           ttGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTg                                                            <  4:263574/149‑1 (MQ=255)                                                                                             gCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGc                                                           <  2:327445/148‑1 (MQ=255)                                                                                              cACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGcaa                                                         >  1:340302/1‑149 (MQ=255)                                                                                                 gCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGca                                                          <  6:138840/145‑1 (MQ=255)                                                                                                 gCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGca                                                          >  5:138840/1‑145 (MQ=255)                                                                                                 gCCGTCAAGGCCGATGAGCAAGGGTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAAc                                                                    >  1:211527/1‑135 (MQ=255)                                                                                                 gCCGTCAAGGCCGATGAGCAAGGGTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAAc                                                                    <  2:211527/135‑1 (MQ=255)                                                                                                    gTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGt                                                   >  5:18216/1‑149 (MQ=255)                                                                                                     tCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTg                                                  <  6:192103/149‑1 (MQ=255)                                                                                                          gCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATg                                             <  2:340302/149‑1 (MQ=255)                                                                                                              aTGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCca                                         <  2:239617/149‑1 (MQ=255)                                                                                                                 aGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACaa                                      <  1:90626/149‑1 (MQ=255)                                                                                                                   cAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATgg                                                                                             <  1:324775/92‑1 (MQ=255)                                                                                                                   cAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATgg                                                                                             >  2:324775/1‑92 (MQ=255)                                                                                                                   cAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACt                                    >  1:237957/1‑149 (MQ=255)                                                                                                                   cAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACt                                    >  4:318482/1‑149 (MQ=255)                                                                                                                    aaGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTg                                   >  4:356834/1‑149 (MQ=255)                                                                                                                     aGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTgg                                  >  8:315443/1‑149 (MQ=255)                                                                                                                              ggcggcTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATAc                         <  7:315443/149‑1 (MQ=255)                                                                                                                              ggcggcTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGCAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACACAACGAGCGTTGCAACAAAGTGGCATTCCCACAACTGGCAGAAATAc                         >  4:272373/1‑149 (MQ=255)                                                                                                                               ccggcACGGGAGCCGGGCCGGGTCACACCCCGCCTGCAGGGCCAGGAACAGCAATAACAGCAAACTACCCAGACGCAAGGCGTCGCATCCTTGAAAGCACGCAACGACCGATCAAACAAAGAGGCATCACCCCAACAGGCAGATATAc                         >  8:343619/2‑148 (MQ=255)                                                                                                                                cggcTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCa                       <  3:243938/149‑1 (MQ=255)                                                                                                                                    tCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAggtg                   >  5:35487/1‑149 (MQ=255)                                                                                                                                      gggTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAggtggt                 >  8:274321/1‑149 (MQ=255)                                                                                                                                       ggTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGtt                <  5:94804/149‑1 (MQ=255)                                                                                                                                          gccgggctgtGACCCGGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGACCCATCCTTGAAAGAACGCAACGAGCGTTGCAACaaa                                                     >  4:138631/1‑109 (MQ=255)                                                                                                                                          gccgggccggGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGc                                 <  6:66102/129‑1 (MQ=255)                                                                                                                                          gccgggccggGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGc                                 >  5:66102/1‑129 (MQ=255)                                                                                                                                            cgggccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACct           <  1:242598/149‑1 (MQ=255)                                                                                                                                            cgggccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACAAGTTCGTATAcc            >  8:293374/1‑148 (MQ=255)                                                                                                                                               gccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCtt        <  4:239790/149‑1 (MQ=255)                                                                                                                                                ccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCttt       <  2:249426/149‑1 (MQ=255)                                                                                                                                                  gggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGc     <  6:109398/149‑1 (MQ=255)                                                                                                                                                   ggTCACAGCCCGGCTGCAGGGCCAGGAACCGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGcc    >  6:199826/1‑149 (MQ=255)                                                                                                                                                     tCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTg  <  6:166707/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TTTACCATTGATCCCCCCGGCGTAAGCCTTGCCGTAGTGCAGCCGCTGTTTCAGGTTGGTCGCATCGCCCAATTGCAGGTACTGCTGCAATTGCACGCCGTCAAGGCCGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTG  >  NC_002947/1108992‑1109288

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