Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,593,269	T→A	7.9%	I102N (ATC→AAC) ‡	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,593,269	0	T	A	7.9%	201.8 / 7.9	76	I102N (ATC→AAC) ‡	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base A (3/3);  total (38/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 = 7.55e-01

TCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCG  >  NC_002947/3593121‑3593415                                                                                                                                                     |                                                                                                                                                    tCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTAt                                                                                                                                                     >  2:150156/1‑149 (MQ=255) tCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTAt                                                                                                                                                     >  2:122293/1‑149 (MQ=255)     cTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATGATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAc                                                                                                                                                 <  3:21927/149‑1 (MQ=255)     cTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAc                                                                                                                                                 <  1:150156/149‑1 (MQ=255)         cACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcc‑tgg                                                                                                                                            <  3:78941/149‑1 (MQ=255)         cACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcc‑tgg                                                                                                                                            <  5:238732/149‑1 (MQ=255)               gCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcc‑tgggcctgg                                                                                                                                      <  1:167153/149‑1 (MQ=255)               gCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcc‑tgggcctg                                                                                                                                       <  4:266791/148‑1 (MQ=255)                    tGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGt                                                                                                                                 <  3:5922/149‑1 (MQ=255)                     gATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTg                                                                                                                                >  5:231558/1‑149 (MQ=255)                     gATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTg                                                                                                                                <  8:155675/149‑1 (MQ=255)                        cAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTgcgc                                                                                                                             >  7:251938/1‑149 (MQ=255)                         aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTgcgc                                                                                                                             >  7:73761/1‑148 (MQ=255)                         aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTgcgc                                                                                                                             <  3:137174/148‑1 (MQ=255)                         aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGcc                                                                                                                            >  2:379926/1‑149 (MQ=255)                         aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTACTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACc‑tgggcctgggc                                                                                                                                    >  7:265324/1‑141 (MQ=255)                         aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTACTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACc‑tgggcctgggc                                                                                                                                    <  8:265324/141‑1 (MQ=255)                          aCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTgcgc                                                                                                                             >  8:339753/1‑147 (MQ=255)                          aCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCa                                                                                                                           <  8:73761/149‑1 (MQ=255)                           cGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAg                                                                                                                          <  5:28703/149‑1 (MQ=255)                           cGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAg                                                                                                                          >  7:129287/1‑149 (MQ=255)                             aaaTCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGc                                                                                                                        <  5:310650/149‑1 (MQ=255)                              aaTCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCa                                                                                                                       <  4:342389/149‑1 (MQ=255)                                tCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCCCg                                                                                                                     >  6:297243/1‑149 (MQ=255)                                 cGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGt                                                                                                                    >  2:89548/1‑149 (MQ=255)                                     aaCGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCgg                                                                                                                >  1:154630/1‑149 (MQ=255)                                          cACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGa                                                                                                           <  8:129287/149‑1 (MQ=255)                                              cGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCCGGCACGTCCGGGCCGACCGt                                                                                                       >  8:71290/1‑149 (MQ=255)                                                acacATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTc                                                                                                     >  1:88879/1‑149 (MQ=255)                                                acacATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTc                                                                                                     <  1:122293/149‑1 (MQ=255)                                                   caTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCgg                                                                                                   <  1:274936/149‑1 (MQ=255)                                                   caTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCggg                                                                                                  <  8:251938/149‑1 (MQ=255)                                                    aTTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCggg                                                                                                  >  1:129847/1‑148 (MQ=255)                                                             ggTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGCGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGcgcg                                                                                         <  1:273520/148‑1 (MQ=255)                                                                    gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgtgccgc                                                                                  <  2:19438/148‑1 (MQ=255)                                                                    gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgc                                                                                  <  1:349805/148‑1 (MQ=255)                                                                    gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGACGTATCGACC‑TGGGCCTCGGCCGTGCGCCAGGCACGAACGGGACGACCGTTCGGGCATTGCGCGgcgccgcg                                                                                 >  2:320862/1‑149 (MQ=39)                                                                       gTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACACTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATAGCGCGgcgccgcgc                                                                                >  2:263273/1‑147 (MQ=255)                                                                        tATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgcgc                                                                                >  3:63657/1‑146 (MQ=255)                                                                          tCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCAt                                                                                               >  1:210567/1‑129 (MQ=255)                                                                          tCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCAt                                                                                               <  2:210567/129‑1 (MQ=255)                                                                           cATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGc                                                                          >  7:64776/1‑149 (MQ=255)                                                                              gctgcCCAACCATGCGCCGTTGGTGGCGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCgg                                                                        >  1:217725/1‑148 (MQ=255)                                                                                  cccAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGAttt                                                                   >  8:2400/1‑149 (MQ=255)                                                                                   ccAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                  >  1:364530/1‑149 (MQ=255)                                                                                    cAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCg                                                                 <  4:116332/149‑1 (MQ=255)                                                                                     aaCCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCgg                                                                >  1:317751/1‑149 (MQ=255)                                                                                         aTGGGCCGTTGGGGGTGGCGGAACAGTTTGGCACCCTGAATACATTCCACCCAGGTCGTATCGAGC‑TGGGCCTGGGGCGTGCGCCAGGCAGGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAg                                                            <  7:194697/149‑1 (MQ=255)                                                                                              ccGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACAt                                                       <  6:328931/149‑1 (MQ=255)                                                                                                   ggtggtggCCGAACTGTTTGGCACCCTCAATAATTTGCACCCTGGTCGTATCGCCC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGACCGGGATTTCGGCCAGGACATCGcc                                                   <  1:320862/148‑1 (MQ=255)                                                                                                    gtggtggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGa                                                 <  2:364530/149‑1 (MQ=255)                                                                                                          gCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCg                                           <  7:2400/149‑1 (MQ=255)                                                                                                          gCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGAGCCAGGCACGTACGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCg                                           >  8:30890/1‑149 (MQ=255)                                                                                                           ccGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                                          <  6:231558/149‑1 (MQ=255)                                                                                                           ccGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                                          >  7:68978/1‑149 (MQ=255)                                                                                                             gAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGa                                        >  7:166420/1‑149 (MQ=255)                                                                                                                   tttGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCt                                  >  4:128817/1‑149 (MQ=255)                                                                                                                     tGGCACCCTCAATACATTGCACCCAGGTCGTCTCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTg                                 <  7:339753/148‑1 (MQ=255)                                                                                                                         aCCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGa                            >  3:58578/1‑149 (MQ=255)                                                                                                                           ccTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGAc                           >  4:247739/1‑148 (MQ=255)                                                                                                                             tCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAAc                        >  7:348689/1‑149 (MQ=255)                                                                                                                              cAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACg                       >  7:172596/1‑149 (MQ=255)                                                                                                                               aaTACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCCTGATTACCTCGCCGACAACg                       >  5:384734/1‑148 (MQ=255)                                                                                                                                    cccaggcccAGGTCGATACGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACg                       >  8:48097/7‑143 (MQ=255)                                                                                                                                     ccaggcccAGGTCGATACGACC‑TGGGTCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTg                                         <  8:82138/119‑1 (MQ=255)                                                                                                                                     ccaggcccAGGTCGATACGACC‑TGGGTCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTg                                         >  7:82138/6‑124 (MQ=255)                                                                                                                                         aCCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCgg            <  1:89548/149‑1 (MQ=255)                                                                                                                                         aCCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCgg            <  1:379926/149‑1 (MQ=255)                                                                                                                                         aCCCAGGTCGTAGCGACC‑TGGGCCTGGCCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCgg            <  7:71290/149‑1 (MQ=255)                                                                                                                                           ccAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt          >  7:160806/1‑149 (MQ=255)                                                                                                                                           ccAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt          <  3:128817/149‑1 (MQ=255)                                                                                                                                           ccAGGTCGATACGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt          >  3:325294/1‑149 (MQ=255)                                                                                                                                               gTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTgcgc      <  5:174719/149‑1 (MQ=255)                                                                                                                                                tCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGcgcg     <  8:172596/149‑1 (MQ=255)                                                                                                                                                   tATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg  >  8:325330/1‑149 (MQ=255)                                                                                                                                                   tATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg  <  2:129847/149‑1 (MQ=255)                                                                                                                                                    aTCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg  >  8:121991/1‑148 (MQ=255)                                                                                                                                                    aTCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg  >  6:358160/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACC‑TGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCG  >  NC_002947/3593121‑3593415

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