Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	392,297	G→A	7.8%	L857L (CTG→CTA) ‡	soxA →	sarcosine oxidase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	392,297	0	G	A	7.8%	184.5 / 4.4	64	L857L (CTG→CTA) ‡	soxA	sarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (34/25);  new base A (3/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

CCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGC  >  NC_002947/392150‑392438                                                                                                                                                    |                                                                                                                                              ccGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGc                                                                                                                                              >  6:35372/1‑149 (MQ=255) ccGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACACCGTACGGCACCGAGACCATGCACGTACTGc                                                                                                                                              >  2:18596/1‑149 (MQ=255)   gcgTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGt                                                                                                                                            >  5:307545/1‑149 (MQ=255)           cGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCagag                                                                                                                                      >  2:117037/1‑147 (MQ=255)              aTCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAggg                                                                                                                                 >  5:55468/1‑149 (MQ=255)                ctcGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCt                                                                                                                               >  5:150563/1‑149 (MQ=255)                 tcGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCtt                                                                                                                              >  1:61079/1‑149 (MQ=255)                   gTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGaa                                                                                                                                    <  7:226548/141‑1 (MQ=255)                   gTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGaa                                                                                                                                    >  8:226548/1‑141 (MQ=255)                          ggCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGt                                                                                                                     <  3:149010/149‑1 (MQ=255)                               gCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGcc                                                                                                                <  6:307545/149‑1 (MQ=255)                                  gTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAg                                                                                                              >  7:105016/1‑148 (MQ=255)                                      tACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                         >  1:317248/1‑149 (MQ=255)                                      tACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                         >  1:312324/1‑149 (MQ=255)                                       aCGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGCGAAGGTCGTCATCATCGTGGGCCATGaca                                                                                                          >  4:230654/1‑147 (MQ=255)                                          aaGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGAt                                                                                                     <  2:179904/149‑1 (MQ=255)                                           aGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTgg                                                                                                                   >  7:63365/1‑134 (MQ=255)                                           aGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTgg                                                                                                                   <  8:63365/134‑1 (MQ=255)                                               aaCGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTc                                                                                                >  6:50021/1‑149 (MQ=255)                                               aaCGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTc                                                                                                >  1:41546/1‑149 (MQ=255)                                                aCGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTc                                                                                                >  4:81100/1‑148 (MQ=255)                                                aCGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCg                                                                                               <  5:35372/149‑1 (MQ=255)                                                  gTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCTTCGAGGCGGGCAAGATGTACATCGTGACCCCGTACGGCAGCGAGACCATGCACGTTCTGCGTGCTGAGAAGGGCTTCATCATCGTGGGCCAGGACACGGATGGTTCGa                                                                                              <  6:50522/148‑1 (MQ=255)                                                  gTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGAt                                                                                             <  1:48685/149‑1 (MQ=255)                                                    gCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGCAGTACGTg                                                                                                                                           <  6:205040/101‑1 (MQ=255)                                                    gCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGCAGTACGTg                                                                                                                                           >  5:205040/1‑101 (MQ=255)                                                      aGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAAc                                                                                         >  3:90057/1‑149 (MQ=255)                                                      aGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAAc                                                                                         <  8:105016/149‑1 (MQ=255)                                                       ggCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                         <  7:167509/132‑1 (MQ=255)                                                       ggCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                         >  8:167509/1‑132 (MQ=255)                                                                cGCCATGGGCGTGCTGGAACATATCATCGAGGCCTGCAAGAAGAACAACCTGACCCCGTACCGCAGCGAGACCATGCCCGTACTGCGTTCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCAATGGTTCGATGAACCCGgacgac                                                                                <  2:312324/148‑1 (MQ=255)                                                                cGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGgacgac                                                                                <  6:150563/148‑1 (MQ=255)                                                                cGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCAATGGTTCGATGAACCCGgacgac                                                                                <  2:317248/148‑1 (MQ=255)                                                                 gCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCt                                                                              >  4:6704/1‑149 (MQ=255)                                                                   cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTtcatca                                                                                                                         >  4:263796/1‑104 (MQ=255)                                                                   cATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTtcatca                                                                                                                         <  3:263796/104‑1 (MQ=255)                                                                         cGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGa                                                                      <  1:117037/149‑1 (MQ=255)                                                                           tGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGc                                                                    <  2:61079/149‑1 (MQ=255)                                                                              tGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTgg                                                                 >  8:64209/1‑149 (MQ=255)                                                                                gAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAAc                                                                          >  3:259689/1‑138 (MQ=255)                                                                                gAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAAc                                                                          <  4:259689/138‑1 (MQ=255)                                                                                         atcGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCg                                                      >  1:296594/1‑149 (MQ=255)                                                                                                  gtcAAGAATTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTTCTGCGTGCAGTGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCTATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGcc                                             <  4:90057/147‑1 (MQ=255)                                                                                                  ggCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGcc                                             <  7:64209/149‑1 (MQ=255)                                                                                                          gTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGt                                     >  4:44966/1‑149 (MQ=255)                                                                                                              aaCCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTtct                                        >  5:315151/1‑142 (MQ=255)                                                                                                              aaCCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTtct                                        <  6:315151/142‑1 (MQ=255)                                                                                                              aaCCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGAt                                 >  5:248348/1‑149 (MQ=255)                                                                                                                 cTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCgg                              <  3:44966/149‑1 (MQ=255)                                                                                                                      cccGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGc                         >  2:296977/1‑149 (MQ=255)                                                                                                                         gTACGGCCCCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGATGGCCGCAACAAGCCGTTCTCGTGTATCGGCCTGCGTg                      >  4:175562/1‑149 (MQ=255)                                                                                                                             ggCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCAt                  >  8:2725/1‑149 (MQ=255)                                                                                                                               cACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACAt                                                                        <  4:197712/93‑1 (MQ=255)                                                                                                                                 ccGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAAc              >  3:82630/1‑149 (MQ=255)                                                                                                                                     gaCCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACcgc           >  1:132394/1‑148 (MQ=255)                                                                                                                                       ccATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGaa        <  7:2725/149‑1 (MQ=255)                                                                                                                                       ccATGCACGTACTGCGTGCAGAGAAGGCCTTCATCATCGTTGGCCAGGTCACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGaa        <  3:230654/149‑1 (MQ=255)                                                                                                                                         ctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTGGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCa                   <  7:67706/135‑1 (MQ=255)                                                                                                                                         ctGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTAGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCa                   >  8:67706/2‑136 (MQ=255)                                                                                                                                         aTGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGa      >  6:56348/1‑149 (MQ=255)                                                                                                                                         aTGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGa      >  1:9469/1‑149 (MQ=255)                                                                                                                                         aTGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCA‑GAACCGCGAAGAc     >  5:58629/1‑149 (MQ=255)                                                                                                                                          tGCACGCAGTACGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGAc     >  1:85674/1‑149 (MQ=255)                                                                                                                                             aCGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACtgc  <  1:296977/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              CCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGC  >  NC_002947/392150‑392438

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