Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,275 G→T 19.2% H164N (CAT→AAT)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,851,2750GT19.2% 895.5 / 139.6 305H164N (CAT→AAT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (103/143);  new base T (28/31);  total (131/174)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

CGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGA  >  CP009273/3851128‑3851423
                                                                                                                                                   |                                                                                                                                                    
  tagtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGcc                                                                                                                                                  <  2:2008461/150‑1 (MQ=255)
   agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGccc                                                                                                                                                 <  1:1083653/150‑1 (MQ=255)
   agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGccc                                                                                                                                                 <  1:3535915/150‑1 (MQ=255)
   agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGccc                                                                                                                                                 <  1:3231164/150‑1 (MQ=255)
   agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGccc                                                                                                                                                 <  1:1676057/150‑1 (MQ=255)
    gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCa                                                                                                                                                >  1:1931814/1‑150 (MQ=255)
    gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCa                                                                                                                                                <  1:63576/150‑1 (MQ=255)
    gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCa                                                                                                                                                <  1:3685092/150‑1 (MQ=255)
      aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAgg                                                                                                                                              <  1:4464288/150‑1 (MQ=255)
      aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAgg                                                                                                                                              >  1:1147712/1‑150 (MQ=255)
      aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAgg                                                                                                                                              >  2:2412572/1‑150 (MQ=255)
       aaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGt                                                                                                                                             >  2:3149364/1‑150 (MQ=255)
       aaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGt                                                                                                                                             >  2:2785864/1‑150 (MQ=255)
       aaaaTCAGCGTATTTCTCTAAAGACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGt                                                                                                                                             <  2:431209/150‑1 (MQ=255)
        aaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTg                                                                                                                                            >  1:2039278/1‑150 (MQ=255)
        aaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTg                                                                                                                                            <  2:4145496/150‑1 (MQ=255)
        aaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGAGTCGCCCAGGTg                                                                                                                                            <  2:2043776/150‑1 (MQ=255)
            cAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAgg                                                                                                                                        <  2:1309865/150‑1 (MQ=255)
              gCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGtt                                                                                                                                      <  2:56892/150‑1 (MQ=255)
                gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                    <  2:3030642/150‑1 (MQ=255)
                  aTTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgt                                                                                                                                  <  1:287694/150‑1 (MQ=255)
                    ttCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgtcg                                                                                                                                >  1:3999921/1‑150 (MQ=255)
                     tctcTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGa                                                                                                                               <  2:2681183/150‑1 (MQ=255)
                      ctctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGat                                                                                                                              >  2:4604848/1‑150 (MQ=255)
                        ctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGatat                                                                                                                            <  1:1588370/150‑1 (MQ=255)
                        ctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGatat                                                                                                                            >  1:999063/1‑150 (MQ=255)
                          aaaTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATgg                                                                                                                          <  1:273492/150‑1 (MQ=255)
                          aaaTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATgg                                                                                                                          >  2:3055658/1‑150 (MQ=255)
                            aTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGtt                                                                                                                        >  2:2624920/1‑150 (MQ=255)
                             tACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTc                                                                                                                       >  1:3442348/1‑150 (MQ=255)
                             tACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTc                                                                                                                       <  2:4367320/150‑1 (MQ=255)
                                 ccACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCt                                                                                                                   <  2:3021524/150‑1 (MQ=255)
                                 ccACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCt                                                                                                                   <  2:2589086/150‑1 (MQ=255)
                                  cACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTg                                                                                                                  <  1:3492375/150‑1 (MQ=255)
                                   aCCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGt                                                                                                                 <  2:3169535/150‑1 (MQ=255)
                                   aCCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGt                                                                                                                 <  2:4076261/150‑1 (MQ=255)
                                     cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                               <  2:901079/150‑1 (MQ=255)
                                      ggACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaaa                                                                                                              <  1:2352994/150‑1 (MQ=255)
                                      ggACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaaa                                                                                                              >  2:4226344/1‑150 (MQ=255)
                                          aGGTGAACGGGTGATGCGGCTCATCATACGNAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCtt                                                                                                          >  1:4284304/1‑150 (MQ=255)
                                          aGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCtt                                                                                                          <  1:3948933/150‑1 (MQ=255)
                                           ggTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTc                                                                                                         <  1:3020610/150‑1 (MQ=255)
                                              gAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGAc                                                                                                      <  2:2477179/150‑1 (MQ=255)
                                                 cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                   <  2:817518/150‑1 (MQ=255)
                                                   ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                                 >  2:1430663/1‑150 (MQ=255)
                                                        tGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAgg                                                                                            >  1:3769557/1‑150 (MQ=255)
                                                         gCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGc                                                                                           >  2:394799/1‑150 (MQ=255)
                                                          cGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGcc                                                                                          <  2:730507/150‑1 (MQ=255)
                                                             cTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCAtt                                                                                       <  1:3652721/150‑1 (MQ=255)
                                                             cTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCAtt                                                                                       >  1:3394524/1‑150 (MQ=255)
                                                              tcatcaTACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTg                                                                                      >  2:3874481/1‑150 (MQ=255)
                                                              tcatcaTACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCCCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTg                                                                                      >  2:2460125/1‑150 (MQ=255)
                                                                 tcatACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTgcgc                                                                                   <  2:1091135/150‑1 (MQ=255)
                                                                   atACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCca                                                                                 <  2:2222098/150‑1 (MQ=255)
                                                                   atACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCca                                                                                 <  1:812756/150‑1 (MQ=255)
                                                                     aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCAACGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                               <  1:4543797/150‑1 (MQ=255)
                                                                      cGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATCTGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAg                                                                              >  2:1608960/1‑150 (MQ=255)
                                                                      cGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAg                                                                              <  1:2568727/150‑1 (MQ=255)
                                                                       gAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAgg                                                                             >  2:3202008/1‑150 (MQ=255)
                                                                        aaaCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGc                                                                            <  1:3293252/150‑1 (MQ=255)
                                                                             accaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtc                                                                       >  2:761628/1‑150 (MQ=255)
                                                                             accaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtc                                                                       >  1:2161639/1‑150 (MQ=255)
                                                                                aTCAGGAAAGGCTCGTCTGCGCGCACAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATctct                                                                    <  1:171133/150‑1 (MQ=255)
                                                                                 tCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCtt                                                                   <  2:3418912/150‑1 (MQ=255)
                                                                                    ggAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttttt                                                                >  1:3097312/1‑150 (MQ=255)
                                                                                      aaaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTcc                                                              >  2:4542042/1‑150 (MQ=255)
                                                                                       aaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCg                                                             <  2:3582798/150‑1 (MQ=255)
                                                                                         ggCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTc                                                           >  1:1378844/1‑150 (MQ=255)
                                                                                          gCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCa                                                          <  2:3296181/150‑1 (MQ=255)
                                                                                              gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                      >  1:3165458/1‑150 (MQ=255)
                                                                                              gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                      <  2:497214/150‑1 (MQ=255)
                                                                                               tCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCg                                                     <  1:2854694/150‑1 (MQ=255)
                                                                                                 tGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCt                                                   <  1:539684/150‑1 (MQ=255)
                                                                                                   cgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTaa                                                 <  1:3849101/150‑1 (MQ=255)
                                                                                                      gcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGAt                                              <  1:1190601/150‑1 (MQ=255)
                                                                                                       cgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATa                                             >  1:596733/1‑150 (MQ=255)
                                                                                                          aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                          <  1:3004596/150‑1 (MQ=255)
                                                                                                          aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                          >  2:872147/1‑150 (MQ=255)
                                                                                                          aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                          <  1:2406673/150‑1 (MQ=255)
                                                                                                          aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                          >  2:2349657/1‑150 (MQ=255)
                                                                                                           ggCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTc                                         >  2:2512554/1‑150 (MQ=255)
                                                                                                              tGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGcc                                      >  1:2460579/1‑150 (MQ=255)
                                                                                                              tGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGcc                                      <  1:1844784/150‑1 (MQ=255)
                                                                                                               gTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGccc                                     >  2:3493275/1‑150 (MQ=255)
                                                                                                                 tGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCa                                   <  2:2912410/150‑1 (MQ=255)
                                                                                                                     gAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGa                               <  1:229564/150‑1 (MQ=255)
                                                                                                                          tCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAg                          <  2:3841271/150‑1 (MQ=255)
                                                                                                                            cACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCCTCGAACCAGTa                        >  1:2952536/1‑150 (MQ=255)
                                                                                                                             aCCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGtca                       >  1:4251947/1‑148 (MQ=255)
                                                                                                                             aCCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTGTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTaa                       >  2:374916/1‑150 (MQ=255)
                                                                                                                               cGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATc                     >  1:1405124/1‑150 (MQ=255)
                                                                                                                                gCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTCATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                    >  1:1220111/1‑150 (MQ=255)
                                                                                                                                gCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                    >  2:3812594/1‑150 (MQ=255)
                                                                                                                                gCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                    <  1:734342/150‑1 (MQ=255)
                                                                                                                                gCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                    <  1:1134418/150‑1 (MQ=255)
                                                                                                                                 cccGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAg                   <  2:2826397/150‑1 (MQ=255)
                                                                                                                                     aTTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTc               >  2:4452555/1‑150 (MQ=255)
                                                                                                                                     aTTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTc               >  1:1103090/1‑150 (MQ=255)
                                                                                                                                      ttGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTcc              >  2:62172/1‑150 (MQ=255)
                                                                                                                                             tACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcg       >  2:1345107/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGA  >  CP009273/3851128‑3851423

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: