Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,735,471 C→G 5.8% S711T (AGT→ACT)  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,4710CG 5.8% 628.6 / 21.1 189S711T (AGT→ACT) sucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base C (91/87);  new base G (3/8);  total (94/95)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 210 total aligned reads displayed.

CGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTACAGA  >  NC_002947/4735323‑4735620
                                                                                                                                                    |                                                                                                                                                     
cgtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACt                                                                                                                                                      <  2:5414065/150‑1 (MQ=255)
 gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                                                     <  2:3884467/150‑1 (MQ=255)
 gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                                                     <  2:2180283/150‑1 (MQ=255)
 gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                                                     <  1:2969693/150‑1 (MQ=255)
 gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                                                     >  1:5154418/1‑150 (MQ=255)
 gtcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                                                     <  2:5103233/150‑1 (MQ=255)
  tcgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTgg                                                                                                                                                    >  1:3002419/1‑150 (MQ=255)
   cgGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGt                                                                                                                                                   <  2:3737325/150‑1 (MQ=255)
    ggTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGCAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGAGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtg                                                                                                                                                  >  2:2308879/1‑150 (MQ=255)
     gTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtga                                                                                                                                                 >  1:3556003/1‑150 (MQ=255)
     gTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtga                                                                                                                                                 <  1:586926/150‑1 (MQ=255)
             aCCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTgg                                                                                                                                         <  2:2012064/150‑1 (MQ=255)
               cTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTc                                                                                                                                       >  2:2248160/1‑150 (MQ=255)
                tGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCg                                                                                                                                      >  2:4169229/1‑150 (MQ=255)
                  gATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAt                                                                                                                                    >  2:1755289/1‑150 (MQ=255)
                    tgttgtGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATca                                                                                                                                  <  2:2529786/150‑1 (MQ=255)
                    tgttgtGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATca                                                                                                                                  >  1:1621973/1‑150 (MQ=255)
                        gtgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATcaccac                                                                                                                              >  1:3819988/1‑150 (MQ=255)
                         tgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACt                                                                                                                             >  1:3667384/1‑150 (MQ=255)
                          gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACtt                                                                                                                            <  1:4386292/150‑1 (MQ=255)
                            tCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgc                                                                                                                          >  2:5160922/1‑150 (MQ=255)
                             cGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTgcg                                                                                                                         >  2:2347157/1‑150 (MQ=255)
                                 gcACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCAcc                                                                                                                     >  2:1330063/1‑150 (MQ=255)
                                  cacaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCg                                                                                                                    >  1:4295423/1‑150 (MQ=255)
                                   acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCTTAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                                                                                                   >  2:3449671/1‑150 (MQ=255)
                                   acaGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGt                                                                                                                   <  2:4674281/150‑1 (MQ=255)
                                          cAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGaa                                                                                                            <  1:709671/150‑1 (MQ=255)
                                                  gCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCGGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGa                                                                                                    <  2:4786172/150‑1 (MQ=255)
                                                      cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                                                                                                >  1:1552298/1‑150 (MQ=255)
                                                      cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                                                                                                <  1:4223139/150‑1 (MQ=255)
                                                      cAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTg                                                                                                <  1:2128475/150‑1 (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)
                                                              cgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTccc                                                                                        <  1:2808173/150‑1 (MQ=255)
                                                              cgGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTccc                                                                                        <  1:2699407/150‑1 (MQ=255)
                                                               ggaggaGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCa                                                                                       <  2:3539402/150‑1 (MQ=255)
                                                               ggaggaGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCa                                                                                       <  1:804717/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)
                                                                   gaatGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGaga                                                                                   <  2:4497488/147‑3 (MQ=255)
                                                                     gTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCAc                                                                                 >  2:4985031/1‑150 (MQ=255)
                                                                     atGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGagaac                                                                                 <  2:1791039/149‑5 (MQ=255)
                                                                     atGCTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGagaac                                                                                 <  2:2657896/149‑5 (MQ=255)
                                                                      tGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCAcc                                                                                >  2:845665/1‑150 (MQ=255)
                                                                        cTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAg                                                                              >  2:1037446/1‑150 (MQ=255)
                                                                        cTCAGGACCCTGACCTTCATAGCCATGCGGCAGAAGGCAGACGAGACCAGACATTCTGAGCCACTTGCGTTCACCCGACGAAATGAACTGATCGAACACGACCTGCGCGCCGTTGGCGAAATCGCCGAACTGGGCTTCCCAGa                                                                                     <  1:2914241/143‑1 (MQ=255)
                                                                         tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc                                                                             >  1:1861684/1‑150 (MQ=255)
                                                                         tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc                                                                             >  2:4010740/1‑150 (MQ=255)
                                                                         tCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTAAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgc                                                                             >  2:4662022/1‑150 (MQ=255)
                                                                          ccGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAgcg                                                                            >  1:1630088/1‑150 (MQ=255)
                                                                                ccTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCg                                                                      >  1:952790/1‑150 (MQ=255)
                                                                                ccTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCg                                                                      >  1:3282681/1‑150 (MQ=255)
                                                                                 cTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCgg                                                                     >  2:226853/1‑150 (MQ=255)
                                                                                       ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggt                                                               >  1:2555285/1‑150 (MQ=255)
                                                                                       ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggt                                                               >  2:1609418/1‑150 (MQ=255)
                                                                                       ttCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggt                                                               >  2:5117930/1‑150 (MQ=255)
                                                                                        tCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCa                                                                   >  2:3591867/1‑145 (MQ=255)
                                                                                          aTAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggtggt                                                            <  2:1230017/150‑1 (MQ=255)
                                                                                             gCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGa                                                         <  2:1861684/150‑1 (MQ=255)
                                                                                               cGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGt                                                       <  1:2235241/150‑1 (MQ=255)
                                                                                                   tgGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTg                                                                                                                                                     >  1:5234912/1‑52 (MQ=255)
                                                                                                       agcagcaTGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTAt                                               >  1:3410033/1‑150 (MQ=255)
                                                                                                            caTGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGaa                                          >  2:1687355/1‑150 (MQ=255)
                                                                                                             aTGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAg                                         <  1:4902704/150‑1 (MQ=255)
                                                                                                                gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGcc                                                   >  1:2311995/1‑137 (MQ=255)
                                                                                                                gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGcc                                      >  2:4731829/1‑150 (MQ=255)
                                                                                                                gTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTGCGGCAGGGTGGTGGAGTAGCCGTATTCGAAGGcc                                      >  2:3984183/1‑150 (MQ=255)
                                                                                                                     aCCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTAc                                 >  2:4486440/1‑150 (MQ=255)
                                                                                                                          acaGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCtt                            <  1:3887617/150‑1 (MQ=255)
                                                                                                                            aGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCt                             >  1:3757008/1‑147 (MQ=255)
                                                                                                                            aGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTcc                          >  2:3041362/1‑150 (MQ=255)
                                                                                                                              gcggcCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTc                        <  1:108668/150‑1 (MQ=255)
                                                                                                                               cggcCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCg                       <  2:1450794/150‑1 (MQ=255)
                                                                                                                                ggcCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGGATGGTGGTGGGGTAGCCGGAGTCGGTGGCGAGTACCGCTGCCTCgg                      >  2:5257252/1‑150 (MQ=255)
                                                                                                                                  ccccACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGAc                    <  1:2864080/150‑1 (MQ=255)
                                                                                                                                   cccACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCAcc                                                                                >  1:4064199/1‑89 (MQ=255)
                                                                                                                                   cccACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACa                   >  1:2752023/1‑150 (MQ=255)
                                                                                                                                      aCTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGa                >  1:504753/1‑150 (MQ=255)
                                                                                                                                        ttGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAAc              >  1:308964/1‑150 (MQ=255)
                                                                                                                                        ttGGGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGTAAc              <  1:3487019/150‑1 (MQ=255)
                                                                                                                                         gttcatCACCAGTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGaa               <  1:667782/143‑1 (MQ=255)
                                                                                                                                           tgCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGaa           <  2:3801463/150‑1 (MQ=255)
                                                                                                                                             cTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGc                                       >  2:2108013/1‑120 (MQ=255)
                                                                                                                                             cTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATc         >  1:1127311/1‑150 (MQ=255)
                                                                                                                                             cTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATc         <  1:5163537/150‑1 (MQ=255)
                                                                                                                                             cTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATc         <  1:1183413/150‑1 (MQ=255)
                                                                                                                                              tCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATct        >  2:690649/1‑149 (MQ=255)
                                                                                                                                                accacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTa      <  2:3733152/150‑1 (MQ=255)
                                                                                                                                                accacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTa      <  1:1775338/150‑1 (MQ=255)
                                                                                                                                                 ccacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGaa           <  1:2250588/144‑1 (MQ=255)
                                                                                                                                                  cacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggtggtg                                                           <  2:4631635/95‑1 (MQ=255)
                                                                                                                                                  cacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggtggtg                                                           >  1:4631635/1‑95 (MQ=255)
                                                                                                                                                  cacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCAtggtggtg                                                           <  2:1891595/95‑1 (MQ=255)
                                                                                                                                                  cacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTAc                                 <  1:2043453/121‑1 (MQ=255)
                                                                                                                                                  cacTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTACa    >  2:4705201/1‑150 (MQ=255)
                                                                                                                                                   acTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTACAg   >  1:3231143/1‑150 (MQ=255)
                                                                                                                                                    cTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTACAGa  <  1:102926/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGATCACCACTTGCGCACCGTTGGCGAAGTCGCCGAACTGGGCTTCCCAGATCACCAGCGCGTTCGGCATGGTGGTGGAGTAGCCGTATTCGAAGGCCAGTACCGCTTCCTCGGACAGGAACGAATCGTACAGA  >  NC_002947/4735323‑4735620

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: