Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 293,197 A→T 6.8% D185V (GAC→GTC)  ssuD → alkanesulfonate monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947293,1970AT 6.8% 248.9 / 7.4 88D185V (GAC→GTC) ssuDalkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (42/40);  new base T (2/4);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CCCGTATCTGGCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGTGATCGTGC  >  NC_002947/293050‑293340
                                                                                                                                                   |                                                                                                                                                 
cccGTATCTGGCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGAc                                                                                                                                                  <  4:78074/149‑1 (MQ=255)
 ccGTATCTGGCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGAcc                                                                                                                                                 >  5:222244/1‑149 (MQ=255)
        tGGCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCAGGTCCTgggcggc                                                                                                                                        >  4:18752/1‑144 (MQ=255)
         ggCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGcc                                                                                                                                        >  2:37281/1‑149 (MQ=255)
         ggCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGAACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGACGCCGCCCAGGACCTGGCGGcc                                                                                                                                        >  3:327922/1‑149 (MQ=255)
          gcgcAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCAGGTCCTGGGCGGcct                                                                                                                                       <  1:204189/149‑2 (MQ=255)
          gcgcAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCAGGTCCTGGGCGGcct                                                                                                                                       <  3:18752/149‑2 (MQ=255)
           cgcAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCAGGTCCTGGGCGGcctc                                                                                                                                      <  7:29214/149‑3 (MQ=255)
            gcAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGcc                                                                                                                                        >  7:35061/1‑146 (MQ=255)
            gcAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGcc                                                                                                                                        <  8:35061/146‑1 (MQ=255)
             cAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCACGACCGGGCGGCCgca                                                                                                                                     >  6:228771/1‑146 (MQ=255)
                ggTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAgg                                                                                                                                 <  5:228771/149‑1 (MQ=255)
                     tGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAg                                                                                                                            >  2:245371/1‑149 (MQ=255)
                       gAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCTCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAg                                                                                                                            <  2:192209/147‑1 (MQ=255)
                           gCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCCCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTac                                                                                                                      >  6:280478/1‑149 (MQ=255)
                              aaaaCGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTacct                                                                                                                    >  5:243695/1‑148 (MQ=255)
                                      gATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTgggg                                                                                                           <  3:280676/149‑1 (MQ=255)
                                          acgacgGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGa                                                                                                        <  1:151651/148‑1 (MQ=255)
                                          acgacgGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGa                                                                                                        <  2:67569/148‑1 (MQ=255)
                                           cgacgGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGc                                                                                                      <  3:232356/149‑1 (MQ=255)
                                            gacgGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTacct                                                                                                                    <  7:88194/134‑1 (MQ=255)
                                            gacgGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTacct                                                                                                                    >  8:88194/1‑134 (MQ=255)
                                                gTAAGCACATCGAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGATCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCg                                                                                                 <  2:249149/149‑1 (MQ=255)
                                                 tAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCgg                                                                                                >  2:346245/1‑149 (MQ=255)
                                                             ggTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGcc                                                                                     >  2:140923/1‑148 (MQ=255)
                                                              gTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCga                                                                                   <  1:257130/149‑1 (MQ=255)
                                                                aaaGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGc                                                                                               <  7:194437/135‑1 (MQ=255)
                                                                aaaGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGc                                                                                               >  8:194437/1‑135 (MQ=255)
                                                                aaaGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCgag                                                                                  >  3:141578/1‑148 (MQ=255)
                                                                  aGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCgag                                                                                  >  2:126389/1‑146 (MQ=255)
                                                                     gcgcCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATc                                                                            <  2:288833/149‑1 (MQ=255)
                                                                      cgcCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTggg                                                                                                            <  1:217525/116‑1 (MQ=255)
                                                                      cgcCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTggg                                                                                                            >  2:217525/1‑116 (MQ=255)
                                                                        ccAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGca                                                                         >  7:185207/1‑149 (MQ=255)
                                                                        ccAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGca                                                                         <  5:280478/149‑1 (MQ=255)
                                                                         cAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCgg                                                                                                >  8:213226/1‑125 (MQ=255)
                                                                         cAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCgg                                                                                                <  7:213226/125‑1 (MQ=255)
                                                                                ctctACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGt                                                                    <  4:327922/146‑1 (MQ=255)
                                                                                ctctACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTgcg                                                                >  5:70484/1‑149 (MQ=255)
                                                                                  ctACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGc                                                                                      >  6:84720/1‑126 (MQ=255)
                                                                                   tACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTgcgc                                                               <  4:141578/147‑1 (MQ=255)
                                                                                   tACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTgcgc                                                               <  7:123219/147‑1 (MQ=255)
                                                                                     cccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGt                                                                    <  7:60394/141‑1 (MQ=255)
                                                                                     cccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGt                                                                    >  8:60394/1‑141 (MQ=255)
                                                                                      ccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCgag                                                                                  <  8:221351/126‑1 (MQ=255)
                                                                                      ccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGaaaa                                                          <  1:37281/149‑1 (MQ=255)
                                                                                      ccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGaaa                                                           >  4:192466/1‑148 (MQ=255)
                                                                                      ccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGaaa                                                           >  8:220238/1‑148 (MQ=255)
                                                                                      ccACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGaaa                                                           >  7:172622/1‑148 (MQ=255)
                                                                                      acACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTAGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCgag                                                                                  >  7:221351/2‑126 (MQ=255)
                                                                                        aCCGTTTCAGCAAGCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTTACCTGGGGCGAGCCACCGGCAGCTg                                                                                          <  2:249885/116‑1 (MQ=255)
                                                                                        aCCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCAAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTg                                                                                          >  1:249885/1‑116 (MQ=255)
                                                                                         ccGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGc                                                       <  1:346245/149‑1 (MQ=255)
                                                                                         ccGATTCAGCAACCACGCCCACCGCTGTATTTAGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGc                                                       <  6:243695/149‑1 (MQ=255)
                                                                                               cAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGACCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGgctgc                                                 >  5:43966/1‑149 (MQ=255)
                                                                                                aGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGgctgct                                                >  7:222259/1‑149 (MQ=255)
                                                                                                 gtaaCCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTc                                               <  1:141405/147‑1 (MQ=255)
                                                                                                   aaCCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTc                                               >  1:331867/1‑147 (MQ=255)
                                                                                                    aCCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAg                                            <  2:62584/149‑1 (MQ=255)
                                                                                                     ccACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAgg                                           >  6:295048/1‑149 (MQ=255)
                                                                                                     ccACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAgg                                           <  6:70484/149‑1 (MQ=255)
                                                                                                           ccACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGcgcg                                      >  1:329485/1‑148 (MQ=255)
                                                                                                            cACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGaa                                    <  8:172622/149‑1 (MQ=255)
                                                                                                               cGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTa                                 >  8:73419/1‑149 (MQ=255)
                                                                                                                  tGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCAc                                                                                                   <  4:302553/81‑1 (MQ=255)
                                                                                                                  tGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCAc                                                                                                   >  3:302553/1‑81 (MQ=255)
                                                                                                                     aTTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCga                                                                                   >  4:267877/1‑94 (MQ=255)
                                                                                                                     aTTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTc                           >  6:166704/1‑149 (MQ=255)
                                                                                                                      tttCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCg                          <  1:140923/149‑1 (MQ=255)
                                                                                                                      tttCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCg                          >  5:288540/1‑149 (MQ=255)
                                                                                                                      tttCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCg                          <  6:288540/149‑1 (MQ=255)
                                                                                                                            cggcTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATcc                    <  1:126389/149‑1 (MQ=255)
                                                                                                                            cggcTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATcc                    <  1:223813/149‑1 (MQ=255)
                                                                                                                             ggcTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCTGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATcc                    <  1:245371/149‑1 (MQ=255)
                                                                                                                               cTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGcc                 >  3:268147/1‑149 (MQ=255)
                                                                                                                                  gTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTg               >  7:10578/1‑148 (MQ=255)
                                                                                                                                  gTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGc              <  3:192466/149‑1 (MQ=255)
                                                                                                                                   tCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCa             >  5:326701/1‑149 (MQ=255)
                                                                                                                                     cGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACg           >  8:157347/1‑149 (MQ=255)
                                                                                                                                      gAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGt          >  3:326873/1‑149 (MQ=255)
                                                                                                                                        ggCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGTGa        <  8:10578/149‑1 (MQ=255)
                                                                                                                                         gccgccCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCg                          >  8:185894/1‑130 (MQ=255)
                                                                                                                                         gccgccCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCg                          <  7:185894/130‑1 (MQ=255)
                                                                                                                                          gccgccAGGTCCTGGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGaa                                    >  3:4488/5‑120 (MQ=255)
                                                                                                                                          gccgccAGGTCCTGGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGaa                                    <  4:4488/116‑1 (MQ=255)
                                                                                                                                          ccgccCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGTGATc      <  8:222259/149‑1 (MQ=255)
                                                                                                                                            gccCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGTGATCGt    >  4:46265/1‑149 (MQ=255)
                                                                                                                                              ccAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGTGATCGTGc  <  5:295048/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCGTATCTGGCGCAAGGTGCTGGAAGGCGAAAACGTCGATTACGACGGTAAGCACATCCAGGTAAAGGGCGCCAAGCTGCTCTACCCACCGATTCAGCAACCACGCCCACCGCTGTATTTCGGCGGCTCGTCCGAGGCCGCCCAGGACCTGGCGGCCGAACAGGTCGAGCTGTACCTGACCTGGGGCGAGCCACCGGCAGCTGTTGCCGAGAAGATCGCACAGGTGCGCGAAAAAGCGGCTGCTCAAGGGCGCGAAGTACGCTTCGGTATCCGCCTGCACGTGATCGTGC  >  NC_002947/293050‑293340

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: