Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,735,448 G→T 5.3% L719M (CTG→ATG)  sucA ← 2‑oxoglutarate dehydrogenase subunit E1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,735,4480GT 5.3% 568.6 / 14.0 170L719M (CTG→ATG) sucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base G (82/79);  new base T (2/7);  total (84/86)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 180 total aligned reads displayed.

CAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCC  >  NC_002947/4735299‑4735596
                                                                                                                                                     |                                                                                                                                                    
cagGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAg                                                                                                                                                      <  1:2802524/150‑1 (MQ=255)
cagGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAg                                                                                                                                                      <  2:2737362/150‑1 (MQ=255)
    tggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcgg                                                                                                                                                  <  1:325479/150‑1 (MQ=255)
     ggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcggc                                                                                                                                                 <  1:3060002/150‑1 (MQ=255)
     ggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcggc                                                                                                                                                 <  1:2090555/150‑1 (MQ=255)
            cTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACtt                                                                                                                                          <  1:569814/150‑1 (MQ=255)
            cTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACtt                                                                                                                                          <  2:3616396/150‑1 (MQ=255)
             tgtgCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTtg                                                                                                                                         >  2:3887617/1‑150 (MQ=255)
               tgCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTtgtg                                                                                                                                       >  1:5103233/1‑150 (MQ=255)
                        cgtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACt                                                                                                                              <  2:5414065/150‑1 (MQ=255)
                        cgtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACt                                                                                                                              <  2:2104891/150‑1 (MQ=255)
                         gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                             <  2:5103233/150‑1 (MQ=255)
                         gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                             >  1:5154418/1‑150 (MQ=255)
                         gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                             <  2:2180283/150‑1 (MQ=255)
                         gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                             >  2:160632/1‑150 (MQ=255)
                         gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                             <  1:2969693/150‑1 (MQ=255)
                          tcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTgg                                                                                                                            >  1:2867422/1‑150 (MQ=255)
                          tcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTgg                                                                                                                            >  1:3002419/1‑150 (MQ=255)
                           cgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGt                                                                                                                           <  1:5193078/150‑1 (MQ=255)
                           cgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGt                                                                                                                           <  2:3737325/150‑1 (MQ=255)
                             gTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtga                                                                                                                         <  1:586926/150‑1 (MQ=255)
                                 gCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGaa                                                                                                                     <  2:4668851/150‑1 (MQ=255)
                                    gACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTg                                                                                                                  <  2:1774293/150‑1 (MQ=255)
                                       cTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTc                                                                                                               >  2:2248160/1‑150 (MQ=255)
                                          gATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAt                                                                                                            >  2:1888992/1‑150 (MQ=255)
                                            tgttgtGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATca                                                                                                          >  1:1621973/1‑150 (MQ=255)
                                            tgttgtGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATca                                                                                                          <  2:2529786/150‑1 (MQ=255)
                                                gtgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATcaccac                                                                                                      >  1:3819988/1‑150 (MQ=255)
                                                 tgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACt                                                                                                     >  1:3667384/1‑150 (MQ=255)
                                                  gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACtt                                                                                                    <  2:4600210/150‑1 (MQ=255)
                                                  gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACtt                                                                                                    <  1:4386292/150‑1 (MQ=255)
                                                   cTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTg                                                                                                   <  1:5282008/150‑1 (MQ=255)
                                                   cTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTg                                                                                                   >  1:283904/1‑150 (MQ=255)
                                                   cTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTg                                                                                                   >  2:1719661/1‑150 (MQ=255)
                                                    tCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgc                                                                                                  <  2:1267132/150‑1 (MQ=255)
                                                    tCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgc                                                                                                  >  2:5160922/1‑150 (MQ=255)
                                                     cGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgcg                                                                                                 >  2:2347157/1‑150 (MQ=255)
                                                      ggCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgcgc                                                                                                >  1:3257747/1‑150 (MQ=255)
                                                          cacaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCg                                                                                            >  1:4295423/1‑150 (MQ=255)
                                                           acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCTTAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                                                                           >  2:1183413/1‑150 (MQ=255)
                                                           acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCTTAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                                                                           >  2:3449671/1‑150 (MQ=255)
                                                           acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                                                                           >  2:2291638/1‑150 (MQ=255)
                                                               cTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGc                                                                                       <  2:4720075/150‑1 (MQ=255)
                                                                     gTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTc                                                                                 <  1:696062/150‑1 (MQ=255)
                                                                       aaCGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGc                                                                               >  1:3539402/1‑150 (MQ=255)
                                                                            tCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAAc                                                                          >  2:1994080/1‑150 (MQ=255)
                                                                            tCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAAc                                                                          >  2:393688/1‑150 (MQ=255)
                                                                              cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                                                                        <  1:4223139/150‑1 (MQ=255)
                                                                              cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                                                                        >  1:1552298/1‑150 (MQ=255)
                                                                               cgACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCAGAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTgg                                                                       <  1:3754091/149‑1 (MQ=255)
                                                                                gACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTggg                                                                      <  2:4295423/150‑1 (MQ=255)
                                                                                     gcgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTcc                                                                 <  2:3257747/150‑1 (MQ=255)
                                                                                      cgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTccc                                                                <  2:2079284/150‑1 (MQ=255)
                                                                                       ggaggaGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCa                                                               <  1:804717/150‑1 (MQ=255)
                                                                                       ggaggaGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCa                                                               <  2:4492968/150‑1 (MQ=255)
                                                                                       ggaggaGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCa                                                               <  2:3539402/150‑1 (MQ=255)
                                                                                       cgatgaatGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCa                                                               <  2:3671187/143‑1 (MQ=255)
                                                                                        gtggagTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCTGACCGCACAGGCGGCGCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAg                                                              <  2:1376279/148‑1 (MQ=255)
                                                                                           gagTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATc                                                           <  1:1719661/150‑1 (MQ=255)
                                                                                             gTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCAc                                                         >  2:4985031/1‑150 (MQ=255)
                                                                                             atGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGagaac                                                         <  2:3861289/149‑5 (MQ=255)
                                                                                             atGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGagaac                                                         <  2:1791039/149‑5 (MQ=255)
                                                                                              tGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCAcc                                                        >  2:845665/1‑150 (MQ=255)
                                                                                                cTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAg                                                      >  2:1037446/1‑150 (MQ=255)
                                                                                                cTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGa                                                             >  2:2914241/1‑143 (MQ=255)
                                                                                                cTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGa                                                             <  1:2914241/143‑1 (MQ=255)
                                                                                                cTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGAGAACCAg                                                      >  2:4913887/1‑150 (MQ=255)
                                                                                                 tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc                                                     >  1:1861684/1‑150 (MQ=255)
                                                                                                 tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc                                                     >  2:4010740/1‑150 (MQ=255)
                                                                                                  ccGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgcg                                                    >  1:1630088/1‑150 (MQ=255)
                                                                                                   cGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtgat                                                                                                                        >  1:2192938/1‑81 (MQ=255)
                                                                                                     gCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGt                                                 >  2:4749886/1‑150 (MQ=255)
                                                                                                       cccTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTc                                               <  1:1187111/150‑1 (MQ=255)
                                                                                                        ccTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCg                                              >  1:3282715/1‑150 (MQ=255)
                                                                                                        ccTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCg                                              >  1:952790/1‑150 (MQ=255)
                                                                                                               ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGtt                                                >  1:2442502/1‑141 (MQ=255)
                                                                                                               ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggt                                       >  2:2639732/1‑150 (MQ=255)
                                                                                                               ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggt                                       >  1:1177048/1‑150 (MQ=255)
                                                                                                               ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggt                                       >  2:5117930/1‑150 (MQ=255)
                                                                                                                tCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCa                                           >  2:3591867/1‑145 (MQ=255)
                                                                                                                tCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCa                                           <  1:3591867/145‑1 (MQ=255)
                                                                                                                    aGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggtggtgg                                  <  1:962826/150‑1 (MQ=255)
                                                                                                                    aGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggtggtgg                                  <  1:3705301/150‑1 (MQ=255)
                                                                                                                     gCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGa                                 <  2:1861684/150‑1 (MQ=255)
                                                                                                                          gtgGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGc                            >  1:730885/1‑150 (MQ=255)
                                                                                                                           tgGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                             <  2:5234912/52‑1 (MQ=255)
                                                                                                                           tgGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGcc                           >  2:1642754/1‑150 (MQ=255)
                                                                                                                                     aTGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAg                 <  1:4902704/150‑1 (MQ=255)
                                                                                                                                      tGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAgg                <  1:4929000/150‑1 (MQ=255)
                                                                                                                                        gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGcc                           <  2:2311995/137‑1 (MQ=255)
                                                                                                                                        gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGcc                           >  1:2311995/1‑137 (MQ=255)
                                                                                                                                        gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGcc              >  2:4731829/1‑150 (MQ=255)
                                                                                                                                        gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTGCGGCAGGGTGGTGGAGTAGCCGTATTCGAAGGcc              >  2:3984183/1‑150 (MQ=255)
                                                                                                                                         tCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCa             >  2:1366364/1‑150 (MQ=255)
                                                                                                                                           aGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGt           >  1:2595552/1‑150 (MQ=255)
                                                                                                                                             aCCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTAc         >  2:4486440/1‑150 (MQ=255)
                                                                                                                                                gCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGc      <  1:1758032/150‑1 (MQ=255)
                                                                                                                                                  acaGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCtt    <  1:3887617/150‑1 (MQ=255)
                                                                                                                                                    aGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCt     >  1:3757008/1‑147 (MQ=255)
                                                                                                                                                    aGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTcc  >  2:3041362/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCC  >  NC_002947/4735299‑4735596

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