Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,089,812 A→T 6.6% T225S (ACC→TCC)  fumC‑I → class 2 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,089,8120AT 6.6% 217.1 / 3.4 74T225S (ACC→TCC) fumC‑Iclass 2 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base T (2/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGT  >  NC_002947/1089665‑1089956
                                                                                                                                                   |                                                                                                                                                
tACGTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAAc                                                                                                                                                 <  5:14212/149‑1 (MQ=255)
     gAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTccgccc                                                                                                                                            >  7:209603/1‑148 (MQ=255)
     gAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCAAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGAACTTGCAGCCGCTGGCCCACTGAGGAAccgcct                                                                                                                                            >  3:202882/1‑148 (MQ=255)
      aaGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGt                                                                                                                                           <  2:201715/149‑1 (MQ=255)
      aaGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGt                                                                                                                                           <  8:120977/149‑1 (MQ=255)
      aaGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGt                                                                                                                                           <  6:331332/149‑1 (MQ=255)
        gACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGGGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgg                                                                                                                                         <  7:151213/149‑1 (MQ=255)
          ccGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgggc                                                                                                                                       <  8:209603/149‑2 (MQ=255)
                   cTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATc                                                                                                                              <  3:47897/149‑1 (MQ=255)
                    tCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATc                                                                                                                              >  8:166943/1‑148 (MQ=255)
                     cACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgg                                                                                                                                         >  8:112042/1‑136 (MQ=255)
                     cACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgg                                                                                                                                         <  7:112042/136‑1 (MQ=255)
                      aCCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAAc                                                                                                                           >  7:232977/1‑149 (MQ=255)
                                   gATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGt                                                                                                              >  5:65744/1‑149 (MQ=255)
                                    aTGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGtt                                                                                                             <  1:11670/149‑1 (MQ=255)
                                        ccGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCg                                                                                                           >  5:301718/1‑147 (MQ=255)
                                        ccGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCg                                                                                                           <  6:301718/147‑1 (MQ=255)
                                        ccGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCt                                                                                                         >  3:20757/1‑149 (MQ=255)
                                        ccGTGAGCATCAGCCAGGTGCTGGATGGCTGGTCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGATGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCt                                                                                                         >  6:125524/1‑149 (MQ=255)
                                          gTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGt                                                                                                       <  3:200672/149‑1 (MQ=255)
                                              gcATGTGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGt                                                                                                   >  2:2394/1‑149 (MQ=255)
                                                aTGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGttt                                                                                                 >  6:10250/1‑149 (MQ=255)
                                                     ccAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTtcgctc                                                                                            <  3:200655/149‑1 (MQ=255)
                                                     ccAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTtcgctc                                                                                            <  5:10250/149‑1 (MQ=255)
                                                       aGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTtcgctcgc                                                                                          <  7:166943/149‑1 (MQ=255)
                                                            cTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGc                                                                                      >  2:333622/1‑148 (MQ=255)
                                                             tGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCt                                                                                     >  5:78789/1‑148 (MQ=255)
                                                               gATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCa                                                                                   <  3:74104/148‑1 (MQ=255)
                                                               gATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCa                                                                                   <  8:316575/148‑1 (MQ=255)
                                                               gATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAg                                                                                  >  6:67433/1‑149 (MQ=255)
                                                                   gctgggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGc                                                                              <  1:2394/149‑1 (MQ=255)
                                                                    ctgggctgCGCAGATAAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGAGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGcc                                                                             >  1:219814/1‑149 (MQ=255)
                                                                     tgggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCt                                                                            >  4:55073/1‑149 (MQ=255)
                                                                                  tCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGa                                                                <  6:65744/148‑1 (MQ=255)
                                                                                  tCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAg                                                               <  1:333622/149‑1 (MQ=255)
                                                                                        ctgcCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCAc                                                          <  8:232977/148‑1 (MQ=255)
                                                                                        ctgcCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCAc                                                          <  5:67433/148‑1 (MQ=255)
                                                                                        ctgcCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTACACg                                                         >  8:50178/1‑149 (MQ=255)
                                                                                          gcCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGAAGTTCGCTGTCGCTTTCGCTCGCCCGCGCAGTGGCCTGACCAAGGTCGACTTAACGcc                                                       >  7:268929/1‑149 (MQ=255)
                                                                                           ccAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGccc                                                      >  5:98978/1‑149 (MQ=255)
                                                                                              aGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGc                                                   <  6:78789/149‑1 (MQ=255)
                                                                                               gTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCt                                                                            >  1:231729/1‑123 (MQ=255)
                                                                                               gTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCt                                                                            <  2:231729/123‑1 (MQ=255)
                                                                                                  gCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAAc                                                                                                                           <  2:162337/73‑1 (MQ=255)
                                                                                                  gCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAg                                                >  1:319756/1‑148 (MQ=255)
                                                                                                   cacaTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCgg                                                    <  5:69786/143‑1 (MQ=255)
                                                                                                   cacaTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCgg                                                    >  6:69786/1‑143 (MQ=255)
                                                                                                    acaTCGAAGGTGGGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAAc                                             <  6:169008/149‑1 (MQ=255)
                                                                                                      aTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGt                                                                                 <  8:98726/111‑1 (MQ=255)
                                                                                                      aTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGt                                                                                 >  7:98726/1‑111 (MQ=255)
                                                                                                       tCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTg                                          >  5:154573/1‑149 (MQ=255)
                                                                                                       tCGAAGGTGTGCTGCCGCACTTGCAGGCGCAGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTg                                          >  1:32000/1‑149 (MQ=255)
                                                                                                         gAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGtt                                        <  6:98978/149‑1 (MQ=255)
                                                                                                          aaGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTc                                       >  8:185398/1‑149 (MQ=255)
                                                                                                           aGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCg                                                                                                                                   <  4:137884/56‑1 (MQ=255)
                                                                                                           aGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCg                                                                                                                                   >  3:137884/1‑56 (MQ=255)
                                                                                                           aGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAATGTCGAGTACACGCCCGGCCAGAACCTGTTCg                                      <  7:50178/149‑1 (MQ=255)
                                                                                                                tgctgcCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGGGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGCACCTGTTCGCACTg                                 >  7:266144/1‑149 (MQ=255)
                                                                                                                    gcCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGc                                                   >  7:110793/1‑127 (MQ=255)
                                                                                                                    gcCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGc                                                   <  8:110793/127‑1 (MQ=255)
                                                                                                                     ccGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATc                              <  2:32000/147‑1 (MQ=255)
                                                                                                                     ccGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCgg                            >  5:291019/1‑149 (MQ=255)
                                                                                                                          cTTGCTGGCGCTGGGCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCGGGAGTTCGCTGTCGGTTTGGCTCGCCATCGCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGc                       <  2:219814/149‑1 (MQ=255)
                                                                                                                             gCAGGAGCTGGCCCAGGGCGGAACCGACGTGGGTACCGGCATCAACGACCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTAGCACTGATCGGCTCGCAg                     <  8:268929/148‑1 (MQ=255)
                                                                                                                                 gcgcTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGt                                         >  8:12377/1‑124 (MQ=255)
                                                                                                                                 gcgcTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGt                                         <  7:12377/124‑1 (MQ=255)
                                                                                                                                   gcTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGacac                <  2:299897/147‑1 (MQ=255)
                                                                                                                                   gcTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGacac                >  1:299897/1‑147 (MQ=255)
                                                                                                                                   gcTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTg              >  8:163242/1‑149 (MQ=255)
                                                                                                                                         ccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTgg        >  6:332438/1‑149 (MQ=255)
                                                                                                                                         ccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTgg        >  5:54989/1‑149 (MQ=255)
                                                                                                                                          cAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGc       <  4:20757/149‑1 (MQ=255)
                                                                                                                                            gggCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCtt     >  2:274988/1‑149 (MQ=255)
                                                                                                                                               cGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGt  <  6:54989/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TACGTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGT  >  NC_002947/1089665‑1089956

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: