Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	813,440	A→G	5.4%	L315L (TTG→CTG)	PP_0700 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	813,440	0	A	G	5.4%	223.9 / 2.1	74	L315L (TTG→CTG)	PP_0700	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base G (2/2);  total (38/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 = 8.82e-01

ATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTC  >  NC_002947/813293‑813585                                                                                                                                                     |                                                                                                                                                  aTGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAgcaactgcaac                                                                                                                                                  >  2:185237/1‑149 (MQ=255)  tGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAgcaactgcaact                                                                                                                                                 <  6:222819/149‑1 (MQ=255)  tGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAgcaactgcaact                                                                                                                                                 >  7:300548/1‑149 (MQ=255)    tAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCaactgcaactgc                                                                                                                                               >  1:348790/1‑149 (MQ=255)             gcgcgcATTATTCTCTTGAACGCCCCCTCTTTCAACACTCATTTC‑AAAAAAAATTTGCACGCTATCAAAAACTTAGCGGCCAGCGCCGGTTTCAGCGCACCACAAGCCAGTGCTCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCa                                                                                                                                      <  1:62216/149‑1 (MQ=255)                  cATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                 <  1:79940/149‑1 (MQ=255)                  cATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                 <  2:100331/149‑1 (MQ=255)                  cATTATTCTCTCGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAg                                                                                                                                 <  7:273688/149‑1 (MQ=255)                   attattCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAgg                                                                                                                                <  7:60591/149‑1 (MQ=255)                      attCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGTGAt                                                                                                                             >  1:244619/1‑149 (MQ=255)                          tctTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTg                                                                                                                           >  5:433/1‑147 (MQ=255)                            ttGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAgg                                                                                                                       >  5:29856/1‑149 (MQ=255)                                aCGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCttt                                                                                                                   <  6:433/149‑1 (MQ=255)                                    ccccTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGc                                                                                                                >  2:66971/1‑148 (MQ=255)                                      ccTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAg                                                                                                             >  2:10211/1‑149 (MQ=255)                                          tGTCAACACTCATTTCAAAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAgg                                                                                                            <  8:320203/147‑1 (MQ=255)                                            tCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTg                                                                                                                  >  4:137078/1‑138 (MQ=255)                                            tCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTg                                                                                                                  <  3:137078/138‑1 (MQ=255)                                              aaCACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTc                                                                                                      >  3:345745/1‑148 (MQ=255)                                               acacTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTc                                                                                                      >  4:37255/1‑147 (MQ=255)                                                 acTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAgg                                                                                                  >  6:312750/1‑149 (MQ=255)                                                   tCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCGGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCt                                                                                                                                         <  2:380363/108‑1 (MQ=255)                                                   tCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCGGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCt                                                                                                                                         >  1:380363/1‑108 (MQ=255)                                                   tCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGCTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGtt                                                                                                >  5:180136/1‑149 (MQ=255)                                                     aTTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGa                                                                                              <  3:8931/149‑1 (MQ=255)                                                      tttcaaaaaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTg                                                                                               >  8:345352/6‑148 (MQ=255)                                                      tttC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGaa                                                                                             <  3:37255/149‑1 (MQ=255)                                                           aaaaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGc                                                                                                                          >  8:125494/1‑116 (MQ=255)                                                           aaaaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGc                                                                                                                          <  7:125494/116‑1 (MQ=255)                                                           aaaaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTcg                                                                                         <  8:352826/149‑1 (MQ=255)                                                           aaaaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCGTTGGCCAGGCTGGTCGAGGTTGAAGTcg                                                                                         >  2:274751/1‑149 (MQ=255)                                                             aaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTcgcc                                                                                       <  8:300548/149‑1 (MQ=255)                                                             aaaaaaTTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTcgcc                                                                                       >  7:223769/1‑149 (MQ=255)                                                                        aCGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGACGAGGTTGAAGTCGCCGCTTACACGGc                                                                            >  2:202960/1‑149 (MQ=255)                                                                           cTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGAt                                                                         <  3:20573/149‑1 (MQ=255)                                                                            tttCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGAt                                                                         >  7:18097/1‑148 (MQ=255)                                                                              tCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGa                                                                       <  7:345352/148‑1 (MQ=255)                                                                                aaaaaCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     <  1:185237/148‑1 (MQ=255)                                                                                aaaaaCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     >  7:19203/1‑148 (MQ=255)                                                                                aaaaaCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     >  2:237335/1‑148 (MQ=255)                                                                                  aaaCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCAc                                                                  >  3:206181/1‑149 (MQ=255)                                                                                  aaaCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCAc                                                                  <  1:10211/149‑1 (MQ=255)                                                                                    aCTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACtt                                                                <  8:18097/149‑1 (MQ=255)                                                                                     cTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTg                                                               <  8:223769/149‑1 (MQ=255)                                                                                      ttAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGt                                                              >  4:235698/1‑149 (MQ=255)                                                                                         gCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtc                                                            <  4:345745/148‑1 (MQ=255)                                                                                         gCGGCCAGCCCCGGCTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtc                                                            <  6:180136/148‑1 (MQ=255)                                                                                                 ccccTGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCATCTTCAACTTCTGTGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTAGATGTTGAAGTCGCCTCTTACACGGCGATGAGCCACTTGTtcatccatcat                                                   <  2:244619/149‑1 (MQ=255)                                                                                                 ccccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGAGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGa                                                                          >  8:272378/1‑126 (MQ=255)                                                                                                 ccccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGAGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGa                                                                          <  7:272378/126‑1 (MQ=255)                                                                                                  cccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCAc                                                                  <  1:57924/133‑1 (MQ=255)                                                                                                  cccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCAc                                                                  >  2:57924/1‑133 (MQ=255)                                                                                                  cccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtcatccatcatc                                                  <  6:29856/149‑1 (MQ=255)                                                                                                   ccGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtcatccatcatc                                                  >  2:90968/1‑148 (MQ=255)                                                                                                    cGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTtcatccatcatcca                                                <  2:348790/149‑1 (MQ=255)                                                                                                             cgcACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     <  4:234183/119‑1 (MQ=255)                                                                                                             cgcACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     >  3:234183/1‑119 (MQ=255)                                                                                                                  cacaAGCCAGTGCCCGAGCATGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAg                                                                                                             >  3:62280/1‑74 (MQ=255)                                                                                                                  cacaAGCCAGTGCCCGAGCATGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAg                                                                                                             <  4:62280/74‑1 (MQ=255)                                                                                                                   acaAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGt                                                              >  4:238156/1‑120 (MQ=255)                                                                                                                   acaAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGt                                                              <  3:238156/120‑1 (MQ=255)                                                                                                                   acaAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTg                                 <  3:235698/149‑1 (MQ=255)                                                                                                                     aaGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTgg                                >  8:341638/1‑148 (MQ=255)                                                                                                                       gCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGCCGCTGCCCTTGGTaa                             >  3:106789/1‑149 (MQ=255)                                                                                                                       gCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTaa                             >  1:355533/1‑149 (MQ=255)                                                                                                                             gCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGc                        >  5:42091/1‑148 (MQ=255)                                                                                                                               ccGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAg                     >  1:85614/1‑149 (MQ=255)                                                                                                                                 gAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCcaca                                              <  1:183214/122‑1 (MQ=255)                                                                                                                                 gAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCcaca                                              >  2:183214/1‑122 (MQ=255)                                                                                                                                   tgcagCCCAGCAGTTGCAGTTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAgg                    >  2:71643/5‑146 (MQ=255)                                                                                                                                   tgcagCCCAGCAGTTGCAGTTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAgg                    <  1:71643/142‑1 (MQ=255)                                                                                                                                    gcagCCCAGCAGTTGCAGTTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGGTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     >  6:343929/4‑96 (MQ=38)                                                                                                                                    gcagCCCAGCAGTTGCAGTTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGCTGAAGTCGCCGCTTACACGGCGATGAGc                                                                     <  5:343929/93‑1 (MQ=38)                                                                                                                                    cACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCt                <  7:341638/149‑1 (MQ=255)                                                                                                                                       gCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTc               >  8:276529/1‑147 (MQ=255)                                                                                                                                       gCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTc               >  2:51207/1‑147 (MQ=255)                                                                                                                                         ccAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAAc                            <  2:185671/132‑1 (MQ=255)                                                                                                                                         ccAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAAc                            >  1:185671/1‑132 (MQ=255)                                                                                                                                               actgcaactgcTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGc     <  4:206181/149‑1 (MQ=255)                                                                                                                                                   caactgcTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTc  >  2:241855/1‑148 (MQ=255)                                                                                                                                                    aactgcTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGgc                   >  2:258010/1‑130 (MQ=255)                                                                                                                                                    aactgcTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGgc                   <  1:258010/130‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTC‑AAAAAAAATTTGCACGCTTTCAAAAACTTAGCGGCCAGCCCCGGTTTCAGCGCACCACAAGCCAGTGCCCGAGCACGCCCAGCAACTGCAACTGCTGGGCTGCAGCCAGGGATTGCAGGCTTTGGCCAGGCTGGTCGAGGTTGAAGTCGCCGCTTACACGGCGATGAGCCACTTGTTCATCCATCATCCACAGACGCTGCCCTTGGTAACCGGCAAGGCGCTCCAGCACCTGCCTC  >  NC_002947/813293‑813585

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