Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,275 G→T 9.7% 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,2750GT 9.7% 958.1 / 48.4 260H164N (CAT→AAT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (93/141);  new base T (9/17);  total (102/158)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

CTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGG  >  CP009273/3851126‑3851422
                                                                                                                                                     |                                                                                                                                                   
cTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATg                                                                                                                                                     <  2:3577178/150‑1 (MQ=255)
      gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCa                                                                                                                                               >  1:1700703/1‑150 (MQ=255)
      gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCa                                                                                                                                               <  2:3736473/150‑1 (MQ=255)
      gtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCa                                                                                                                                               >  1:529845/1‑150 (MQ=255)
            aTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGaa                                                                                                                                         <  2:278010/150‑1 (MQ=255)
                 cGTATTTCTCTAAATACTCCACCGGACAGGTGANCGGGTGATGNGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGttt                                                                                                                                    <  1:880956/150‑1 (MQ=255)
                 cGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGttt                                                                                                                                    >  2:3143319/1‑150 (MQ=255)
                  gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                   <  2:649003/150‑1 (MQ=255)
                  gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                   <  1:3653688/150‑1 (MQ=255)
                  gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                   <  2:2894333/150‑1 (MQ=255)
                     tttCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcgtc                                                                                                                                >  2:3027301/1‑150 (MQ=255)
                       tctcTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGa                                                                                                                              >  2:1255211/1‑150 (MQ=255)
                       tctcTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGa                                                                                                                              >  2:3348417/1‑150 (MQ=255)
                        ctctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGat                                                                                                                             <  1:239602/150‑1 (MQ=255)
                           tAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATg                                                                                                                          <  2:3106843/150‑1 (MQ=255)
                           tAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATg                                                                                                                          <  1:357682/150‑1 (MQ=255)
                              aTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGtt                                                                                                                       <  1:1014341/150‑1 (MQ=255)
                              aTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGtt                                                                                                                       <  2:796994/150‑1 (MQ=255)
                               tACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTc                                                                                                                      >  2:1232375/1‑150 (MQ=255)
                                  tCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                   >  2:1561892/1‑150 (MQ=255)
                                   ccACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCt                                                                                                                  >  1:4048415/1‑150 (MQ=255)
                                   ccACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCt                                                                                                                  >  2:3225951/1‑150 (MQ=255)
                                    cACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTg                                                                                                                 <  1:2251522/150‑1 (MQ=255)
                                     aCCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGt                                                                                                                <  2:3456298/150‑1 (MQ=255)
                                      ccGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTa                                                                                                               >  2:40858/1‑150 (MQ=255)
                                        ggACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaaa                                                                                                             >  1:43224/1‑150 (MQ=255)
                                         gACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAAt                                                                                                            >  2:563495/1‑150 (MQ=255)
                                           cAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCt                                                                                                          >  1:1269209/1‑150 (MQ=255)
                                           cAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCt                                                                                                          >  2:2638582/1‑150 (MQ=255)
                                             ggTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTc                                                                                                        <  2:1024286/150‑1 (MQ=255)
                                               tGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGa                                                                                                      >  2:286575/1‑150 (MQ=255)
                                                 aaCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACg                                                                                                    <  1:305314/150‑1 (MQ=255)
                                                 aaCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACg                                                                                                    <  1:3971795/150‑1 (MQ=255)
                                                   cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                  >  1:2944988/1‑150 (MQ=255)
                                                    gggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTg                                                                                                 >  1:245019/1‑150 (MQ=255)
                                                        gATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTa                                                                                             <  2:2284713/150‑1 (MQ=255)
                                                         aTGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAg                                                                                            >  1:2198345/1‑150 (MQ=255)
                                                          tGCGGCTCATCATACCAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAgg                                                                                           >  1:2271428/1‑150 (MQ=255)
                                                             ggCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCa                                                                                        >  2:2170256/1‑150 (MQ=255)
                                                             ggCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCa                                                                                        <  1:2046880/150‑1 (MQ=255)
                                                             ggCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCa                                                                                        <  1:2866106/150‑1 (MQ=255)
                                                                 catcatACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTgc                                                                                    >  2:1794337/1‑150 (MQ=255)
                                                                  atcatACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTgcg                                                                                   <  1:2324842/150‑1 (MQ=255)
                                                                    catACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGcc                                                                                 <  1:2121495/150‑1 (MQ=255)
                                                                     atACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCca                                                                                <  1:1436410/150‑1 (MQ=255)
                                                                       aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGGTTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                              <  2:3820416/150‑1 (MQ=255)
                                                                        cGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAg                                                                             >  1:594781/1‑150 (MQ=255)
                                                                        cGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAg                                                                             <  1:1139050/150‑1 (MQ=255)
                                                                          aaaCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGc                                                                           <  1:3558090/150‑1 (MQ=255)
                                                                          aaaCCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGc                                                                           >  2:1746006/1‑150 (MQ=255)
                                                                                ccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtct                                                                     <  1:822265/150‑1 (MQ=255)
                                                                                  aTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATctct                                                                   >  1:370607/1‑150 (MQ=255)
                                                                                    cAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttt                                                                 <  1:3607823/150‑1 (MQ=255)
                                                                                      ggAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttttt                                                               <  1:583487/150‑1 (MQ=255)
                                                                                       gAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGNAAATCCACCNCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTc                                                              <  1:2521760/150‑1 (MQ=255)
                                                                                       gAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTc                                                              <  2:1008342/150‑1 (MQ=255)
                                                                                       gAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTc                                                              <  1:1704667/150‑1 (MQ=255)
                                                                                           ggCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTc                                                          <  2:901767/150‑1 (MQ=255)
                                                                                            gCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCa                                                         <  1:1999865/150‑1 (MQ=255)
                                                                                             cTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAg                                                        <  2:477564/150‑1 (MQ=255)
                                                                                              tCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGt                                                       <  2:2904396/150‑1 (MQ=255)
                                                                                               cGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGtt                                                      <  1:1493877/150‑1 (MQ=255)
                                                                                                gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                     >  2:281173/1‑150 (MQ=255)
                                                                                                   tGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCTCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCt                                                  >  1:3680443/1‑150 (MQ=255)
                                                                                                     cgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTaa                                                <  1:2492750/150‑1 (MQ=255)
                                                                                                      gcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAg                                               <  1:3091025/150‑1 (MQ=255)
                                                                                                         cgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATa                                            <  2:3609168/150‑1 (MQ=255)
                                                                                                         cgcAGGCGGTTGCAGAACCTCCACCGCCCGATTGCTGATACGAGGCGCCCAGGGGAAGGTTTCGTCGCTAAGGTTCGCCTGTAACTCTTCGACGCCGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTGTCCGTCCGTTCGCTAAGATa                                            <  2:679174/150‑1 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                         >  2:972253/1‑150 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                         <  2:2078774/150‑1 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                         <  2:1731008/150‑1 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCATGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                         <  2:1427974/150‑1 (MQ=255)
                                                                                                             ggCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGTCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTc                                        <  1:3364641/150‑1 (MQ=255)
                                                                                                              gCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCg                                       <  1:616203/150‑1 (MQ=255)
                                                                                                               cTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGc                                      >  2:699178/1‑150 (MQ=255)
                                                                                                               cTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGc                                      >  1:3473078/1‑150 (MQ=255)
                                                                                                                tGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGcc                                     >  2:2592996/1‑150 (MQ=255)
                                                                                                                 gTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGccc                                    >  2:3281565/1‑150 (MQ=255)
                                                                                                                 gTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGccc                                    <  1:3422340/150‑1 (MQ=255)
                                                                                                                   tGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCa                                  >  2:745025/1‑150 (MQ=255)
                                                                                                                      aGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCg                               >  2:2442226/1‑150 (MQ=255)
                                                                                                                        aaaaTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTACGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCTCTAAGATAGTTCGCCCCATCGaa                             >  2:1068811/1‑150 (MQ=255)
                                                                                                                          aaTCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAAcc                           >  1:4166740/1‑150 (MQ=255)
                                                                                                                           aTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCa                          >  1:1033245/1‑150 (MQ=255)
                                                                                                                            tCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAg                         <  2:1358287/150‑1 (MQ=255)
                                                                                                                                ccGCCCGAGTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAAt                     <  2:4157225/150‑1 (MQ=255)
                                                                                                                                  gCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCa                   <  2:1771699/150‑1 (MQ=255)
                                                                                                                                    ccGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGc                 <  1:2586313/150‑1 (MQ=255)
                                                                                                                                     cGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCg                <  2:3035920/150‑1 (MQ=255)
                                                                                                                                       aTTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTc              >  1:374592/1‑150 (MQ=255)
                                                                                                                                        ttGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTcc             <  1:3187784/150‑1 (MQ=255)
                                                                                                                                        ttGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTcc             <  1:1817643/150‑1 (MQ=255)
                                                                                                                                           cTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCAc          >  1:3808584/1‑150 (MQ=255)
                                                                                                                                            tGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACt         >  2:1824652/1‑150 (MQ=255)
                                                                                                                                              aTACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTcc       <  1:1887586/150‑1 (MQ=255)
                                                                                                                                              aTACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTcc       <  2:3941800/150‑1 (MQ=255)
                                                                                                                                               tACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcg      <  1:194510/150‑1 (MQ=255)
                                                                                                                                                  gATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcggcg   <  1:4079373/150‑1 (MQ=255)
                                                                                                                                                  gATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcggcg   >  1:1357568/1‑150 (MQ=255)
                                                                                                                                                   aTGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcggcgg  >  1:1772279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGG  >  CP009273/3851126‑3851422

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: