Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,546,056 A→G 6.2% V302A (GTG→GCG)  PP_3131 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,546,0560AG 6.2% 188.6 / 2.5 64V302A (GTG→GCG) PP_3131membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

CGATGAACGGCAGCAGGTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCACGCTCATCACGTAAC  >  NC_002947/3545913‑3546197
                                                                                                                                                     |                                                                                                                                              
cGATGAACGGCAGCAGGTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGgcg                                                                                                                                           >  3:208460/1‑149 (MQ=255)
 gATGAACGGCAGCAGGTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGgcgc                                                                                                                                          <  6:364564/149‑1 (MQ=255)
       cGGCAGCAGGTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATcc                                                                                                                                    >  3:207550/1‑149 (MQ=255)
        ggCAGCAGGTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGCAGAAAGTAGCAGCGCCAGAGCCGATTGACGCCGCATTACCCCAGGGCGATGAAAACCACCACGAAGACCACgccgcc                                                                                                                                         >  1:115262/1‑148 (MQ=255)
               ggTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCa                                                                                                                            >  1:116246/1‑149 (MQ=255)
                   cAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCAcgc                                                                                                                         >  8:46973/1‑148 (MQ=255)
                   cAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCAcgc                                                                                                                         >  3:127496/1‑148 (MQ=255)
                   cAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCAcgc                                                                                                                         >  2:13487/1‑148 (MQ=255)
                    aGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCTGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGa                                                                                                                       >  1:238907/1‑149 (MQ=255)
                         cTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCTGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATg                                                                                                                  <  2:238907/149‑1 (MQ=255)
                         cTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACTGGc                                                                                                                              <  7:238832/137‑1 (MQ=255)
                         cTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGc                                                                                                                              >  8:238832/1‑137 (MQ=255)
                          ttttCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATgg                                                                                                                 >  8:107618/1‑149 (MQ=255)
                                 ccccAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCTACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGaa                                                                                                           <  1:268879/149‑1 (MQ=255)
                                  cccAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGaa                                                                                                           >  3:171603/1‑147 (MQ=255)
                                  cccAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACa                                                                                                         <  1:13487/149‑1 (MQ=255)
                                   ccAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAg                                                                                                        >  2:81044/1‑149 (MQ=255)
                                        cACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCcgg                                                                                                   <  2:116246/149‑1 (MQ=255)
                                               gTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCcag                                                                                            >  1:89851/1‑149 (MQ=255)
                                                      aGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAgg                                                                                     <  4:207550/149‑1 (MQ=255)
                                                              gCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCAt                                                                            >  5:179817/1‑149 (MQ=255)
                                                               ccAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATg                                                                           <  7:107618/149‑1 (MQ=255)
                                                                cAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATgg                                                                          >  8:263871/1‑149 (MQ=255)
                                                                    aTGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCtt                                                                      >  7:228565/1‑149 (MQ=255)
                                                                       cctgccGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTtgtg                                                                   <  5:290180/149‑1 (MQ=255)
                                                                       cctgccGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTtgt                                                                    <  2:89851/148‑1 (MQ=255)
                                                                        ctgccGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGc                                                                  >  6:363188/1‑149 (MQ=255)
                                                                           ccGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCtt                                                                      >  7:253389/1‑142 (MQ=255)
                                                                           ccGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCtt                                                                      <  8:253389/142‑1 (MQ=255)
                                                                            cGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCg                                                              <  1:81044/149‑1 (MQ=255)
                                                                             gC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGcc                                                                        >  4:57341/1‑138 (MQ=255)
                                                                             gC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGcc                                                                        <  3:57341/138‑1 (MQ=255)
                                                                              c‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCggg                                                            >  2:65683/1‑149 (MQ=255)
                                                                                    cTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGt                                                         >  6:261371/1‑149 (MQ=255)
                                                                                       gcagcGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGc                                                      <  2:115262/149‑1 (MQ=255)
                                                                                        cgctgctaGGTGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGca                                                     <  7:40985/142‑1 (MQ=255)
                                                                                        cagcGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCg                                                              <  8:345321/140‑1 (MQ=255)
                                                                                        cagcGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCg                                                              >  7:345321/1‑140 (MQ=255)
                                                                                        cagcGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGca                                                     <  4:127496/149‑1 (MQ=255)
                                                                                             ccAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCTACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAg                                                <  5:261371/149‑1 (MQ=255)
                                                                                             ccAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAACGCCCGCTGcacaagg                                                <  4:208460/149‑4 (MQ=255)
                                                                                                ggCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGccc                                             >  1:239467/1‑149 (MQ=255)
                                                                                                     gcTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTg                                           >  3:100582/1‑148 (MQ=255)
                                                                                                         gACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCGCGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACg                                       >  2:278079/1‑148 (MQ=255)
                                                                                                         gACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCGCGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACg                                       >  2:282857/1‑148 (MQ=255)
                                                                                                          aCGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGcaca                                                   >  3:166758/1‑135 (MQ=255)
                                                                                                          aCGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGcaca                                                   <  4:166758/135‑1 (MQ=255)
                                                                                                              cgcCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACgaacagg                                 <  6:179817/149‑1 (MQ=255)
                                                                                                                 cTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACgaacaggaac                              <  7:263871/149‑1 (MQ=255)
                                                                                                                    tGCCCAGGGCGATGAACACCACCACGATGACCGCGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACgaacaggaac                              <  1:282857/147‑1 (MQ=255)
                                                                                                                    tGCCCAGGGCGATGAACACCACCACGATGACCGCGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACgaacaggaac                              <  1:278079/147‑1 (MQ=255)
                                                                                                                      cccAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGca                                                     >  6:302385/1‑121 (MQ=255)
                                                                                                                      cccAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGca                                                     <  5:302385/121‑1 (MQ=255)
                                                                                                                          gggCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATg                     <  1:65683/149‑1 (MQ=255)
                                                                                                                              gATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCa                 <  8:228565/149‑1 (MQ=255)
                                                                                                                                tGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCACg               <  4:171603/149‑1 (MQ=255)
                                                                                                                                   acacCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGAt                      >  7:39262/1‑139 (MQ=255)
                                                                                                                                   acacCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGAt                      <  8:39262/139‑1 (MQ=255)
                                                                                                                                      ccaccacGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCAc                >  6:275059/1‑142 (MQ=255)
                                                                                                                                      ccaccacGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCAc                <  5:275059/142‑1 (MQ=255)
                                                                                                                                        accacGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCACGCTCATCAc       <  7:46973/149‑1 (MQ=255)
                                                                                                                                              gATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGc                          >  6:143979/1‑125 (MQ=255)
                                                                                                                                              gATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGc                          <  5:143979/125‑1 (MQ=255)
                                                                                                                                              gATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCACGCTCATCACGTAAc  >  7:189138/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CGATGAACGGCAGCAGGTCCAGCACCTTTTCGCCCCCAGGCACGACGGTCAGCAAGGGCAACGCCAGCATGCCTGCCGC‑‑‑ACCTAGCAGCGCCAGGCGCTTGACGGCGCCTTTGCCCAGGGCGATGAACACCACCACGATGACCACGGCGCCGATCCACGGGCCACGCGACAATGGCGCGAACAGCCCGGCGGCCAGCAGCAGGGCGCCCATGGCCTTGTGCAGCGGGCGTTGCACATAGCCCTGCACGAACAGGAACAGGCCGATGGCCACGCTCATCACGTAAC  >  NC_002947/3545913‑3546197

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: