Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 392,295 C→G 7.6% L857V (CTG→GTG) ‡ soxA → sarcosine oxidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947392,2950CG 7.6% 199.7 / 3.6 65L857V (CTG→GTG) ‡soxAsarcosine oxidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base G (3/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGC  >  NC_002947/392148‑392438
                                                                                                                                                   |                                                                                                                                               
gCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACt                                                                                                                                                <  5:100518/149‑1 (MQ=255)
  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                                                                                                                                    >  8:226548/1‑141 (MQ=255)
                    gTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGaa                                                                                                                                    <  7:226548/141‑1 (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                                                                                                         >  8:167509/1‑132 (MQ=255)
                                                        ggCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGacac                                                                                                         <  7:167509/132‑1 (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                                                                          <  4:259689/138‑1 (MQ=255)
                                                                                 gAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAAc                                                                          >  3:259689/1‑138 (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)
                                                                                                                                                   |                                                                                                                                               
GCCCGCGTGTTCCGTATCTCGTTCACCGGCGAGCTGTCGTACGAAGTGAACGTGCAGGCCAACTACGCCATGGGCGTGCTGGAACAGATCATCGAGGCCGGCAAGAAGTACAACCTGACCCCGTACGGCACCGAGACCATGCACGTACTGCGTGCAGAGAAGGGCTTCATCATCGTGGGCCAGGACACCGATGGTTCGATGAACCCGGACGACCTGAACATGAGCTGGTGCGTTGGCCGCAACAAGCCGTTCTCGTGGATCGGCCTGCGTGGCATGAACCGCGAAGACTGC  >  NC_002947/392148‑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: ‑