Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,089,817 A→T 5.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,8170AT 5.6% 391.9 / 7.9 124T225S (ACC→TCC) fumC‑Iclass 2 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (63/54);  new base T (2/5);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 125 total aligned reads displayed.

TACGTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGG  >  NC_002947/1089670‑1089965
                                                                                                                                                   |                                                                                                                                                    
   gTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAAccgcc                                                                                                                                                 >  1:1145594/1‑150 (MQ=255)
   gTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAAccgcc                                                                                                                                                 >  1:5767/1‑150 (MQ=255)
   gTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCCCGCTGCAATGCTGCCAAGTCGCACATCGAAAGTGTGTATCCTGTCTTGTAGGCGCTGGCCCAGGGCGGAAccgcc                                                                                                                                                 >  1:1738357/1‑150 (MQ=255)
    tGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAAccgccg                                                                                                                                                >  1:2234800/1‑150 (MQ=255)
     gAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCCGGCGATGGCCCACGGCGGTTccgccct                                                                                                                                               <  2:1354262/150‑3 (MQ=255)
       aGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgg                                                                                                                                             <  1:984618/150‑1 (MQ=255)
            ggTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTAcc                                                                                                                                        <  1:113150/150‑1 (MQ=255)
                 tACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCAt                                                                                                                                   <  2:1438498/150‑1 (MQ=255)
                  aCTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATc                                                                                                                                  >  2:1201413/1‑150 (MQ=255)
                    tCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCaa                                                                                                                                >  1:1005256/1‑150 (MQ=255)
                     cACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAAc                                                                                                                               <  1:347170/150‑1 (MQ=255)
                      aCCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACg                                                                                                                              >  2:2284007/1‑150 (MQ=255)
                        cTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCAGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGcc                                                                                                                            >  1:1949236/1‑150 (MQ=255)
                         tGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGGTGCTGGACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGccc                                                                                                                           >  1:398113/1‑150 (MQ=255)
                         tGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCGAGGGCGGAACCGCCGTGGGTACCGGCATCAAc                                                                                                                               >  2:1943456/1‑146 (MQ=255)
                         tGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCGAGGGCGGAACCGCCGTGGGTACCGGCATCAAc                                                                                                                               <  1:1943456/146‑1 (MQ=255)
                         tGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgg                                                                                                                                             <  1:699631/132‑1 (MQ=255)
                             ggACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATcc                                                                                                                       >  1:1265516/1‑150 (MQ=255)
                              gACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCg                                                                                                                      >  2:766858/1‑150 (MQ=255)
                               aCGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGc                                                                                                                     >  1:1037497/1‑150 (MQ=255)
                                cgcgATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCa                                                                                                                    <  2:2234800/150‑1 (MQ=255)
                                   gATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGtt                                                                                                                 >  2:828801/1‑150 (MQ=255)
                                        ccGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGCTCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTg                                                                                                            >  2:924824/1‑150 (MQ=255)
                                        ccGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCa                                                                                                                    >  2:1307226/1‑142 (MQ=255)
                                        actggcgcATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCa                                                                                                                    <  1:1307226/138‑1 (MQ=255)
                                            gcgcATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCgg                                                                                                        >  1:368433/1‑150 (MQ=255)
                                              gcATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTg                                                                                                                                              <  1:610096/110‑1 (MQ=255)
                                              gcATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGtt                                                                                                      >  2:471471/1‑150 (MQ=255)
                                               cATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGttt                                                                                                     <  2:1265516/150‑1 (MQ=255)
                                                aTGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTtc                                                                                                    >  1:2276486/1‑150 (MQ=255)
                                                aTGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTtc                                                                                                    >  2:818432/1‑150 (MQ=255)
                                                 gcacccAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCACGGCGGTTCCGCCCTgg                                                                                                                                             >  1:2016065/5‑108 (MQ=255)
                                                    gCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTGGAAGAACGAAGGTGCGCTGCCGCACTTGCACGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTtcgctc                                                                                                >  2:1081422/1‑150 (MQ=255)
                                                          tGCTGGATGGCTGTGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGGTGCCGAACCTGGTGGCGCTGGCCCAGGGCGGAACCGCCGTGCGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGc                                                                                          <  1:549413/150‑1 (MQ=255)
                                                           gCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATc                                                                                                                        <  1:929365/119‑1 (MQ=255)
                                                           gCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATc                                                                                                                        >  2:929365/1‑119 (MQ=255)
                                                             tGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCa                                                                                       >  1:949126/1‑150 (MQ=255)
                                                               gATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGt                                                                                     <  1:718973/150‑1 (MQ=255)
                                                                aTGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTg                                                                                    >  1:391374/1‑150 (MQ=255)
                                                                aTGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTg                                                                                    <  2:2285428/150‑1 (MQ=255)
                                                                   gctgggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGGACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGcc                                                                                 <  1:2087748/150‑1 (MQ=255)
                                                                     tgggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTg                                                                               >  2:2015678/1‑150 (MQ=255)
                                                                     tgggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTg                                                                               >  1:653404/1‑150 (MQ=255)
                                                                       ggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGAc                                                                             <  1:766858/150‑1 (MQ=255)
                                                                       ggctgCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGAc                                                                             <  2:2276486/150‑1 (MQ=255)
                                                                           gcgcAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAg                                                                         <  2:391374/150‑1 (MQ=255)
                                                                           gcgcAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAg                                                                         <  2:1145594/150‑1 (MQ=255)
                                                                             gcAGATCAATGCTGCCAAGTCGCACATCGAAGGTATGCTTCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGt                                                                       >  2:1725304/1‑150 (MQ=255)
                                                                                gATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTc                                                                      <  1:2228117/148‑1 (MQ=255)
                                                                                gATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGa                                                                    <  1:828801/150‑1 (MQ=255)
                                                                                 aTCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCAAGGGCGACACCGCGGCAGACATCGGCATTACCGCTCTTTCGCAGTTCTCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAg                                                                   <  1:295456/150‑1 (MQ=255)
                                                                                  tCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGt                                                                  >  2:2152098/1‑150 (MQ=255)
                                                                                         tgcCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGcc                                                           >  2:981829/1‑150 (MQ=255)
                                                                                           ccAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTc                                                                                                          <  1:792901/101‑1 (MQ=255)
                                                                                           ccAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTc                                                                                                          >  2:792901/1‑101 (MQ=255)
                                                                                             aaGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCAt                                                                                                                                   >  2:334096/1‑74 (MQ=255)
                                                                                             aaGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCAt                                                                                                                                   <  1:334096/74‑1 (MQ=255)
                                                                                             aaGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCg                                                                                                               >  1:1247907/1‑94 (MQ=255)
                                                                                             aaGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCg                                                                                                               <  2:1247907/94‑1 (MQ=255)
                                                                                             aaGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGc                                                       >  2:1865848/1‑150 (MQ=255)
                                                                                                tCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAg                                                    >  2:1969378/1‑150 (MQ=255)
                                                                                                  gCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGaa                                                  <  1:2015678/150‑1 (MQ=255)
                                                                                                  gCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGaa                                                  >  2:1692342/1‑150 (MQ=255)
                                                                                                     caTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCt                                               <  1:1201413/150‑1 (MQ=255)
                                                                                                      aTCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTg                                              >  1:844287/1‑150 (MQ=255)
                                                                                                        cGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGCCGAGTTCACGCCCGGCCAGAACCTGtt                                            >  1:1288896/1‑150 (MQ=255)
                                                                                                         gAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTc                                           >  1:123842/1‑150 (MQ=255)
                                                                                                            ggTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCa                                        <  1:981829/150‑1 (MQ=255)
                                                                                                             gtgtGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCAc                                       <  2:653404/150‑1 (MQ=255)
                                                                                                                tgctgcCGCACTTGCAGGCGCTGGCCCAGGGAGGAACCGCCGTGGGTACCGGCATCAACGCACATCAGCAATGCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGa                                    >  2:1479174/1‑150 (MQ=255)
                                                                                                                  ctgcCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGAt                                   <  2:2107141/149‑1 (MQ=255)
                                                                                                                  ctgcCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGACCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGAt                                   >  1:2107141/1‑149 (MQ=255)
                                                                                                                        cACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTg                                                                                                            >  2:1733274/1‑70 (MQ=255)
                                                                                                                        cACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTg                                                                                                            <  1:1733274/70‑1 (MQ=255)
                                                                                                                          cTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCa                          <  2:5767/150‑1 (MQ=255)
                                                                                                                           ttGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCAc                                                              <  2:498646/113‑1 (MQ=255)
                                                                                                                               aGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGaca                     >  1:524397/1‑150 (MQ=255)
                                                                                                                               aGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGaca                     >  1:1477705/1‑150 (MQ=255)
                                                                                                                                ggCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGt                                                                                     >  2:2140290/1‑85 (MQ=255)
                                                                                                                                ggCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGt                                                                                     <  1:2140290/85‑1 (MQ=255)
                                                                                                                                    cTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGcc                <  1:1725304/150‑1 (MQ=255)
                                                                                                                                     tGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCg               >  2:779094/1‑150 (MQ=255)
                                                                                                                                      ggCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCAGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGt              <  1:422296/150‑1 (MQ=255)
                                                                                                                                       gCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTg             >  1:2089145/1‑150 (MQ=255)
                                                                                                                                        cccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTgg            <  2:1738357/150‑1 (MQ=255)
                                                                                                                                         ccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGc           >  2:664719/1‑150 (MQ=255)
                                                                                                                                         ccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGc           <  2:1037497/150‑1 (MQ=255)
                                                                                                                                         ccAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGc           >  1:1769588/1‑150 (MQ=255)
                                                                                                                                          cAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCTCAGTTCGCTGTCGGTGTCGCTCGCCAGCTCAGTGGCCTAACCAAGGTCGAGTTCACGCCCGCCCAGAACCTGTTCGCACTGATCGGCTCGCAGGCCACTGCCGTGGCt          >  2:413015/1‑150 (MQ=255)
                                                                                                                                          cAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGc                               <  2:2087627/129‑1 (MQ=255)
                                                                                                                                           aGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCtt         >  2:1193045/1‑150 (MQ=255)
                                                                                                                                              gCGGAACCGTCGTGGGTACCGGCATCAACGCCCATCCCAAACGTGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGt      <  1:2201976/150‑1 (MQ=255)
                                                                                                                                              gCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGt      <  2:368433/150‑1 (MQ=255)
                                                                                                                                               cGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGATTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTc     <  2:1183658/150‑1 (MQ=255)
                                                                                                                                                 gAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTg   <  1:664719/150‑1 (MQ=255)
                                                                                                                                                 gAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTg   <  2:949126/150‑1 (MQ=255)
                                                                                                                                                 gAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGAACGCTGCCGTGGCTTTGTCTg   >  2:624917/1‑150 (MQ=255)
                                                                                                                                                 gAACCGCCGTGGGTACCGACATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCt                                               <  2:472178/106‑1 (MQ=255)
                                                                                                                                                 gAACCGCCGTGGGTACCGACATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCt                                               >  1:472178/1‑106 (MQ=255)
                                                                                                                                                  aaCCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTTCCGTGGCTTTGTCTgg  >  1:1402359/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACGTGAAGACCGGTCGTACTCACCTGATGGACGCGATGCCCGTGCGCATGAGCCAGGTGCTGGATGGCTGGGCTGCGCAGATCAATGCTGCCAAGTCGCACATCGAAGGTGTGCTGCCGCACTTGCAGGCGCTGGCCCAGGGCGGAACCGCCGTGGGTACCGGCATCAACGCCCATCCGCAGTTCGCTGTCGGTTTCGCTCGCCAGCTCAGTGGCCTGACCAAGGTCGAGTTCACGCCCGGCCAGAACCTGTTCGCACTGATCGGCTCGCAGGACACTGCCGTGGCTTTGTCTGG  >  NC_002947/1089670‑1089965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: