Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,135,838 A→C 6.5% Y171D (TAC→GAC)  hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,135,8380AC 6.5% 638.5 / 36.2 230Y171D (TAC→GAC) hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base A (86/129);  new base C (5/10);  total (91/139)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 231 total aligned reads displayed.

CCAGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGC  >  CP009273/3135694‑3135985
                                                                                                                                                |                                                                                                                                                   
ccAGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGcc                                                                                                                                                  >  1:1414879/1‑148 (MQ=255)
  aGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCaa                                                                                                                                                >  2:359669/1‑147 (MQ=255)
  aGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCAt                                                                                                                                                <  1:1234623/148‑1 (MQ=255)
  aGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCAt                                                                                                                                                <  2:375975/148‑1 (MQ=255)
  aGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCAt                                                                                                                                                >  1:722594/1‑148 (MQ=255)
   ggTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCAtt                                                                                                                                               <  1:891028/148‑1 (MQ=255)
   ggTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCAtt                                                                                                                                               <  2:150245/148‑1 (MQ=255)
   ggTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCAtt                                                                                                                                               <  1:889705/148‑1 (MQ=255)
    gTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTa                                                                                                                                              <  1:2019319/148‑1 (MQ=255)
     ttCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAg                                                                                                                                             >  1:2683784/1‑148 (MQ=255)
     ttCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAg                                                                                                                                             >  1:2549644/1‑148 (MQ=255)
     ttCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAg                                                                                                                                             >  1:1341514/1‑148 (MQ=255)
     ttCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAg                                                                                                                                             >  2:1406816/1‑148 (MQ=255)
         gAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAGCGaa                                                                                                                                         >  2:1579606/1‑148 (MQ=255)
         gAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGaa                                                                                                                                         >  1:2914799/1‑148 (MQ=255)
         gAATGAGCGGCGTTATACCGCCAAGCAGCGCCCCGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGAAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGCATGACCCCAGTATCCATTAGCGaa                                                                                                                                         >  2:597435/1‑148 (MQ=255)
           aTGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGaaaa                                                                                                                                       <  1:2675120/148‑1 (MQ=255)
           aTGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGaaaa                                                                                                                                       >  1:414280/1‑148 (MQ=255)
            tgtgCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAAt                                                                                                                                      <  2:1049904/148‑1 (MQ=255)
              tgCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATAc                                                                                                                                    <  2:3437795/148‑1 (MQ=255)
                  gcgTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCaa                                                                                                                                <  1:2293705/148‑1 (MQ=255)
                     ttttACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAGCGAAAATACCCAACTg                                                                                                                             <  1:729302/148‑1 (MQ=255)
                     ttttACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAGCGAAAATACCCAACTg                                                                                                                             <  1:195665/148‑1 (MQ=255)
                      tttACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGa                                                                                                                            <  2:2265332/148‑1 (MQ=255)
                               cagcagCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCa                                                                                                                   >  2:3338492/1‑148 (MQ=255)
                                  cagcGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAAcc                                                                                                                <  1:3590793/148‑1 (MQ=255)
                                     cgcCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGa                                                                                                             <  1:367049/148‑1 (MQ=255)
                                     cgcCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGa                                                                                                             <  1:1566216/148‑1 (MQ=255)
                                       ccACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATc                                                                                                           >  1:485443/1‑148 (MQ=255)
                                       ccACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATc                                                                                                           >  1:2570618/1‑148 (MQ=255)
                                       ccACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATc                                                                                                           <  1:1086271/148‑1 (MQ=255)
                                        cACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCt                                                                                                          <  2:2838359/148‑1 (MQ=255)
                                        cACGACGCGGTTAGCGTAACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCt                                                                                                          <  1:3028481/148‑1 (MQ=255)
                                         acgacgCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCtt                                                                                                         <  2:2977788/148‑1 (MQ=255)
                                         acgacgCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCtt                                                                                                         <  2:421088/148‑1 (MQ=255)
                                         acgacgCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCtt                                                                                                         <  2:1852822/148‑1 (MQ=255)
                                          cgacgCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCttt                                                                                                        >  1:1756215/1‑148 (MQ=255)
                                            acgCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGa                                                                                                      <  2:3813430/148‑1 (MQ=255)
                                                 gTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTg                                                                                                 <  2:1800897/148‑1 (MQ=255)
                                                  ttAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGt                                                                                                >  2:954159/1‑148 (MQ=255)
                                                   tAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGtt                                                                                               <  1:638714/148‑1 (MQ=255)
                                                      cGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTg                                                                                            >  2:3302827/1‑148 (MQ=255)
                                                      cGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTg                                                                                            >  1:546630/1‑148 (MQ=255)
                                                       gTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGAAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGa                                                                                           <  1:1707740/148‑1 (MQ=255)
                                                       gTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGAAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGa                                                                                           <  1:1707160/148‑1 (MQ=255)
                                                        tCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGaa                                                                                          <  1:1379117/148‑1 (MQ=255)
                                                         cACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTACTGAAc                                                                                         >  1:2418890/1‑148 (MQ=255)
                                                          aCGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACGTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACt                                                                                        <  2:2335018/148‑1 (MQ=255)
                                                           cGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACtt                                                                                       >  2:2418232/1‑148 (MQ=255)
                                                           cGCTGGCACTCCAACGCCTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACtt                                                                                       >  2:3283536/1‑148 (MQ=255)
                                                               ggCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGt                                                                                   <  2:1555494/148‑1 (MQ=255)
                                                               ggCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGt                                                                                   <  2:1875769/148‑1 (MQ=255)
                                                                 cACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGa                                                                                 <  2:3182276/148‑1 (MQ=255)
                                                                 cACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGa                                                                                 <  2:1277779/148‑1 (MQ=255)
                                                                    tCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCACCTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAAct                                                                              <  2:2093280/148‑1 (MQ=255)
                                                                      cAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAActct                                                                            <  2:3764297/148‑1 (MQ=255)
                                                                       aaCGCGTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCtt                                                                           <  2:1003647/148‑1 (MQ=255)
                                                                            ttGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGa                                                                      <  1:1542970/148‑1 (MQ=255)
                                                                            ttGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGa                                                                      <  2:1736110/148‑1 (MQ=255)
                                                                               cAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTg                                                                   <  2:1777406/148‑1 (MQ=255)
                                                                               cAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTg                                                                   <  1:2979450/148‑1 (MQ=255)
                                                                               cAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTg                                                                   <  1:1796509/148‑1 (MQ=255)
                                                                                   tAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCa                                                               <  1:3435999/148‑1 (MQ=255)
                                                                                        cgcTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGc                                                          <  2:414291/148‑1 (MQ=255)
                                                                                          cTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCa                                                        <  1:1502979/148‑1 (MQ=255)
                                                                                          cTACCGCAATCAGGGTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTCGCCATTAGCGAAAATACCCACCTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCa                                                        <  2:2751274/148‑1 (MQ=255)
                                                                                           tACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAg                                                       >  1:1140020/1‑148 (MQ=255)
                                                                                            aCCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAgg                                                      >  1:443741/1‑148 (MQ=255)
                                                                                            aCCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAgg                                                      <  2:1068588/148‑1 (MQ=255)
                                                                                             ccGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGt                                                     >  1:3374108/1‑148 (MQ=255)
                                                                                                   tCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAAc                                               >  2:2098990/1‑148 (MQ=255)
                                                                                                   tCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAAc                                               >  1:1438099/1‑148 (MQ=255)
                                                                                                     aGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGc                                             <  1:400872/148‑1 (MQ=255)
                                                                                                      ggTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACTcc                                            <  1:182786/148‑1 (MQ=255)
                                                                                                       gTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCt                                           >  1:2130962/1‑148 (MQ=255)
                                                                                                       gTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCt                                           >  1:3260353/1‑148 (MQ=255)
                                                                                                        ttCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCtt                                          >  1:1097405/1‑148 (MQ=255)
                                                                                                         tCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCttt                                         >  2:2830452/1‑148 (MQ=255)
                                                                                                              tCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCa                                    <  1:1877696/148‑1 (MQ=255)
                                                                                                              tCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCa                                    >  1:646268/1‑148 (MQ=255)
                                                                                                              tCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCa                                    >  1:3537475/1‑148 (MQ=255)
                                                                                                               ccGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCAt                                   <  2:1185343/148‑1 (MQ=255)
                                                                                                               ccGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCAt                                   <  2:3034058/148‑1 (MQ=255)
                                                                                                                cggcggCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCAtt                                  <  2:2598125/148‑1 (MQ=255)
                                                                                                                   cggcAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCg                               >  1:2151007/1‑148 (MQ=255)
                                                                                                                     gcAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGa                             >  1:1847339/1‑148 (MQ=255)
                                                                                                                       aGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAgg                           <  2:3429134/148‑1 (MQ=255)
                                                                                                                        gTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGc                          <  1:1646269/148‑1 (MQ=255)
                                                                                                                         tttCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCt                         >  1:1099914/1‑148 (MQ=255)
                                                                                                                         tttCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCt                         >  1:3584856/1‑148 (MQ=255)
                                                                                                                         tttCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCt                         >  1:3779474/1‑148 (MQ=255)
                                                                                                                            cATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTttg                      <  2:3537475/148‑1 (MQ=255)
                                                                                                                            cATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTttg                      <  1:1686928/148‑1 (MQ=255)
                                                                                                                              tCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTttggt                    <  2:2001077/148‑1 (MQ=255)
                                                                                                                              tCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTttggt                    <  2:70044/148‑1 (MQ=255)
                                                                                                                               cGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTttggtt                   >  2:1056116/1‑148 (MQ=255)
                                                                                                                                gCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTttggttg                  <  2:2079016/148‑1 (MQ=255)
                                                                                                                                    gATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTGTGGTTGGGTc              >  1:2532898/1‑148 (MQ=255)
                                                                                                                                        aCCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCt          <  1:1330672/148‑1 (MQ=255)
                                                                                                                                                aGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTtgcagtgc  <  2:2940190/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCAGGTTCTGAATGTGCGGCGTTTTACCGCCCAGCAGCGCCACGACGCGGTTAGCGTCACGCTGGCACTCCAACGCTTGCAGGTAGTGCGCTACCGCAATCAGGTTCACTTCCGGCGGCAGTTTCATCGCCGGATGACCCCAGTAGCCATTAGCGAAAATACCCAACTGACCGCTGGCAACCAGATCTTTGATCTTGTTCTGAACTTTGGTGAACTCTTCCGGACTGTTCAGGTGCCAGGTCGAAACGCCTTTCAGCATTTCGGAGGCTTTGGTTGGGTCAGCTTGCAGTGC  >  CP009273/3135694‑3135985

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: