Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,773,938 C→T 6.7% G1777S (GGC→AGC)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,773,9380CT 6.7% 176.6 / 2.2 62G1777S (GGC→AGC) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/2);  total (31/31)
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.55e-01

CCACTGCACACCACGCGGCGAACGCTGCGGCACTGCTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTGCAGGCGGTTGA  >  NC_002947/4773790‑4774084
                                                                                                                                                    |                                                                                                                                                  
ccACTGCACACCACGCGGCGAACGCTGCGGCACTGCTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGcc                                                                                                                                                    >  6:232134/1‑149 (MQ=255)
                         tGCGGCACTGCTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGgcgcgcg                                                                                                                           >  5:27588/1‑149 (MQ=255)
                             gCACTGCTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCCCGGCCAGGCGCGCGCCgg                                                                                                                       >  8:267183/1‑149 (MQ=255)
                               aCTGCTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGt                                                                                                                      >  7:146432/1‑148 (MQ=255)
                                   cTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGtcggca                                                                                                                 >  8:247659/1‑145 (MQ=255)
                                   cTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGtcat                                                                                                                   >  4:58629/1‑147 (MQ=255)
                                   cTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGtcat                                                                                                                   >  8:189619/1‑147 (MQ=255)
                                      gACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCTGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCTCGCCGGTCAGCAGCg                                                                                                              <  2:190293/149‑1 (MQ=255)
                                      gACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAg                                                                                                                <  5:132075/147‑1 (MQ=255)
                                        ccTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGc                                                                                                            >  2:58774/1‑149 (MQ=255)
                                          tCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGc                                                                                                            >  4:133755/1‑147 (MQ=255)
                                           cGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCTGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGc                                                                                                            >  6:13803/1‑146 (MQ=255)
                                             gcgcTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGGCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCagaa                                                                                                       <  3:221037/149‑1 (MQ=255)
                                              cgcTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCTGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCagaag                                                                                                      <  5:13803/149‑1 (MQ=255)
                                                cTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAgaagaa                                                                                                    >  3:159498/1‑149 (MQ=255)
                                                 tGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAAc                                                                                                   <  6:101883/149‑1 (MQ=255)
                                                   gCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAAc                                                                                                   <  1:58774/147‑1 (MQ=255)
                                                    ccATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTcc                                                                                                <  3:58629/149‑1 (MQ=255)
                                                     cATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTcc                                                                                                >  8:195046/1‑148 (MQ=255)
                                                       tGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGgcgcgcgc                                                                                                                          <  5:294356/120‑1 (MQ=255)
                                                       tGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGgcgcgcgc                                                                                                                          >  6:294356/1‑120 (MQ=255)
                                                         aaGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCacac                                                                                           <  4:159498/149‑1 (MQ=255)
                                                            gCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGa                                                                                        <  3:133755/149‑1 (MQ=255)
                                                              cTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGaca                                                                                      >  8:71161/1‑149 (MQ=255)
                                                              cTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGaca                                                                                      >  2:302857/1‑149 (MQ=255)
                                                                gcagcaTCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCTGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGAcaca                                                                                    >  5:316058/1‑149 (MQ=255)
                                                                  agcaTCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATc                                                                                  >  5:148009/1‑149 (MQ=255)
                                                                     aTCGAAGGCATGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGaa                                                                               <  2:12409/149‑1 (MQ=255)
                                                                      tCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAg                                                                              <  3:66961/149‑1 (MQ=255)
                                                                      tCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAg                                                                              >  4:55754/1‑149 (MQ=255)
                                                                      tCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAg                                                                              >  2:19812/1‑149 (MQ=255)
                                                                       cGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGc                                                                             >  4:133236/1‑149 (MQ=255)
                                                                              aCGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCgg                                                                                                             >  1:76903/1‑110 (MQ=255)
                                                                              aCGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCgg                                                                                                             <  2:76903/110‑1 (MQ=255)
                                                                              aCGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACgg                                                                      >  1:227127/1‑149 (MQ=255)
                                                                                    tGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGtttt                                                                <  8:146432/149‑1 (MQ=255)
                                                                                      cAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCt                                                              >  1:101068/1‑149 (MQ=255)
                                                                                      cAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCt                                                              <  6:27588/149‑1 (MQ=255)
                                                                                       aGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTg                                                             >  7:103031/1‑149 (MQ=255)
                                                                                         ggTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTgca                                                           <  3:55754/149‑1 (MQ=255)
                                                                                               gACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCg                                                     >  5:234374/1‑149 (MQ=255)
                                                                                                    ggtggtCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCg                                                >  7:79942/1‑149 (MQ=255)
                                                                                                      tggtCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGTAGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCt                                              <  1:119539/149‑1 (MQ=255)
                                                                                                        gtCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTg                                             <  1:302857/148‑1 (MQ=255)
                                                                                                           gatgatGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCAt                                         >  8:261248/1‑149 (MQ=255)
                                                                                                            atgatgGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATc                                        <  2:227127/149‑1 (MQ=255)
                                                                                                                tgGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAg                                     <  5:232134/148‑1 (MQ=255)
                                                                                                                 ggTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTc                                                                                                 <  3:353153/87‑1 (MQ=255)
                                                                                                                 ggTCTGGCACAGGAGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTc                                                                                                 >  4:353153/1‑87 (MQ=255)
                                                                                                                  ggcTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGcc                                  <  7:336401/147‑1 (MQ=255)
                                                                                                                  gTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGc                                   <  8:103031/148‑1 (MQ=255)
                                                                                                                  gTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGc                                   <  2:101068/148‑1 (MQ=255)
                                                                                                                   tCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCa                                 <  7:195046/149‑1 (MQ=255)
                                                                                                                         acaGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGgcg                           >  7:326329/1‑149 (MQ=255)
                                                                                                                             gcgcTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCtt                       <  8:326329/149‑1 (MQ=255)
                                                                                                                                   cGGGTCGCGGTGATCGCTGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATcc                 <  6:316058/149‑1 (MQ=255)
                                                                                                                                    gggTCGCTGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCa                <  7:261248/149‑1 (MQ=255)
                                                                                                                                       tCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTg             <  1:19812/149‑1 (MQ=255)
                                                                                                                                       tCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTg             >  6:66927/1‑149 (MQ=255)
                                                                                                                                          cgGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTGCAg          >  6:116001/1‑149 (MQ=255)
                                                                                                                                           ggTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTGCAgg         >  3:139054/1‑149 (MQ=255)
                                                                                                                                                  gCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTGCAGGCGGTTGa  <  7:71161/149‑1 (MQ=255)
                                                                                                                                                  gCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTGCAGGCGGTTGa  <  7:189619/149‑1 (MQ=255)
                                                                                                                                                   ccGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCAt                   >  1:19124/1‑131 (MQ=255)
                                                                                                                                                   ccGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCAt                   <  2:19124/131‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCACTGCACACCACGCGGCGAACGCTGCGGCACTGCTCGACCTCGGCGCTGGCCATGAAGGCCTGCAGCATCGAAGGCACGAAGTGCAGGGTGCTGACCTGGTGGTCGATGATGGTCTGGCACAGGCGCTGCGGGTCGCGGTGATCGCCGGGCAAGGCCACGGCCAGGCGCGCGCCGGTCATCAGCGGCCAGAAGAACTCCCACACCGACACATCGAAGCTGAACGGGGTTTTCTGCAGCACCGTGTCGCTGGCATCCAGGCCATAGGCGCCTTGCATCCACTGCAGGCGGTTGA  >  NC_002947/4773790‑4774084

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