Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,275 G→T 30.1% 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,2750GT30.1% 671.1 / 284.3 342H164N (CAT→AAT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (119/119);  new base T (45/59);  total (164/178)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

CTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGA  >  CP009273/3851126‑3851423
                                                                                                                                                     |                                                                                                                                                    
cTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATg                                                                                                                                                      <  1:2319185/150‑1 (MQ=255)
 tCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgc                                                                                                                                                     <  1:2216330/150‑1 (MQ=255)
 tCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgc                                                                                                                                                     <  2:935650/150‑1 (MQ=255)
  cGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgcg                                                                                                                                                    >  1:703032/1‑150 (MQ=255)
  cGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgcg                                                                                                                                                    >  2:4432712/1‑150 (MQ=255)
  cGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgcg                                                                                                                                                    >  1:2916124/1‑150 (MQ=255)
       taAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAg                                                                                                                                               >  1:233730/1‑150 (MQ=255)
        aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAgg                                                                                                                                              <  2:2299223/150‑1 (MQ=255)
        aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAgg                                                                                                                                              <  2:2621181/150‑1 (MQ=255)
        aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAgg                                                                                                                                              >  1:993553/1‑150 (MQ=255)
        aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAgg                                                                                                                                              >  2:3134306/1‑150 (MQ=255)
          aaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTg                                                                                                                                            <  2:2184713/150‑1 (MQ=255)
           aaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGa                                                                                                                                           <  1:1625453/150‑1 (MQ=255)
           aaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGa                                                                                                                                           <  2:1619234/150‑1 (MQ=255)
            aTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGaa                                                                                                                                          >  2:2125328/1‑150 (MQ=255)
                 cGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGttt                                                                                                                                     >  1:483060/1‑150 (MQ=255)
                         tctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGata                                                                                                                             <  2:1865730/150‑1 (MQ=255)
                          ctAACTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGatat                                                                                                                            <  1:3742316/150‑1 (MQ=255)
                          ctAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGatat                                                                                                                            >  1:1839093/1‑150 (MQ=255)
                           tAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATg                                                                                                                           >  1:4649716/1‑150 (MQ=255)
                             aagaCTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGt                                                                                                                         <  1:3689839/147‑1 (MQ=255)
                             aaTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGt                                                                                                                         >  1:2803845/1‑150 (MQ=255)
                               tACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTc                                                                                                                       >  1:4523903/1‑150 (MQ=255)
                                 cTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGc                                                                                                                     >  1:573897/1‑150 (MQ=255)
                                 cTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGc                                                                                                                     >  1:818394/1‑150 (MQ=255)
                                   ccACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCGGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCTTCGATATGGTTCGCCt                                                                                                                   >  2:2814699/1‑150 (MQ=255)
                                     aCCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGt                                                                                                                 >  2:3291706/1‑150 (MQ=255)
                                      ccGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTa                                                                                                                <  1:3363627/150‑1 (MQ=255)
                                      ccGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTa                                                                                                                >  1:160164/1‑150 (MQ=255)
                                       cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAGTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                               >  2:1053259/1‑150 (MQ=255)
                                       cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                               >  2:1475000/1‑150 (MQ=255)
                                           cAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCt                                                                                                           <  1:2603333/150‑1 (MQ=255)
                                            aGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCtt                                                                                                          >  2:3317314/1‑150 (MQ=255)
                                            aGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCtt                                                                                                          >  2:920705/1‑150 (MQ=255)
                                               tGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGa                                                                                                       >  1:3765103/1‑150 (MQ=255)
                                                 aaCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACg                                                                                                     <  2:4600672/150‑1 (MQ=255)
                                                 aaCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACg                                                                                                     <  1:3035/150‑1 (MQ=255)
                                                   cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                   >  1:4687787/1‑150 (MQ=255)
                                                    gggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTg                                                                                                  >  1:3878972/1‑150 (MQ=255)
                                                     ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                                 <  1:1554582/150‑1 (MQ=255)
                                                     ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                                 <  1:743726/150‑1 (MQ=255)
                                                     ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                                 <  1:3132208/150‑1 (MQ=255)
                                                      gTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGtt                                                                                                <  1:1039192/150‑1 (MQ=255)
                                                       tGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGttt                                                                                               <  1:3416364/150‑1 (MQ=255)
                                                          tGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAgg                                                                                            >  1:3964177/1‑150 (MQ=255)
                                                            cGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGcc                                                                                          >  1:3295775/1‑150 (MQ=255)
                                                                tcatcaTACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTg                                                                                      <  1:60806/150‑1 (MQ=255)
                                                                tcatcaTACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTg                                                                                      <  1:4418165/150‑1 (MQ=255)
                                                                    catACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGcc                                                                                  <  1:2417367/150‑1 (MQ=255)
                                                                    catACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGcc                                                                                  <  2:3765103/150‑1 (MQ=255)
                                                                     atACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCca                                                                                 <  2:533752/150‑1 (MQ=255)
                                                                      tACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcac                                                                                >  1:2114376/1‑150 (MQ=255)
                                                                       aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                               <  2:3681125/150‑1 (MQ=255)
                                                                       aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                               <  2:1273052/150‑1 (MQ=255)
                                                                        cGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAg                                                                              <  2:3036861/150‑1 (MQ=255)
                                                                         gAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAgg                                                                             >  2:3866614/1‑150 (MQ=255)
                                                                                ccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtct                                                                      <  2:1308318/150‑1 (MQ=255)
                                                                                 caTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtctc                                                                     <  2:171091/150‑1 (MQ=255)
                                                                                    cAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttt                                                                  <  1:4424688/150‑1 (MQ=255)
                                                                                     aGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCtttt                                                                 <  1:829688/150‑1 (MQ=255)
                                                                                     aGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCtttt                                                                 <  2:289520/150‑1 (MQ=255)
                                                                                     aGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCtttt                                                                 <  1:207628/150‑1 (MQ=255)
                                                                                      ggAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttttt                                                                <  2:3789689/150‑1 (MQ=255)
                                                                                        aaaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTcc                                                              <  2:2630938/150‑1 (MQ=255)
                                                                                        aaaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTcc                                                              >  2:2499611/1‑150 (MQ=255)
                                                                                         aaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCg                                                             >  2:4262202/1‑150 (MQ=255)
                                                                                             cTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAg                                                         <  2:3076946/150‑1 (MQ=255)
                                                                                              tCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGt                                                        <  1:2000146/150‑1 (MQ=255)
                                                                                              tCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGt                                                        <  1:4264439/150‑1 (MQ=255)
                                                                                                   tGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCt                                                   >  1:2719653/1‑150 (MQ=255)
                                                                                                       cgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGa                                               <  2:993514/150‑1 (MQ=255)
                                                                                                         cgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATa                                             <  2:1082885/150‑1 (MQ=255)
                                                                                                         cgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATa                                             >  2:4099590/1‑150 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                          <  1:2926965/150‑1 (MQ=255)
                                                                                                              gCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCg                                        <  2:230464/150‑1 (MQ=255)
                                                                                                                      aGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCg                                >  2:1032899/1‑150 (MQ=255)
                                                                                                                      aGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCg                                <  2:3539725/150‑1 (MQ=255)
                                                                                                                        aaaaTCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGaa                              >  2:2042678/1‑150 (MQ=255)
                                                                                                                        aaaaTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGaa                              >  1:1122794/1‑150 (MQ=255)
                                                                                                                         aaaTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAAc                             >  1:4721890/1‑150 (MQ=255)
                                                                                                                         aaaTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAAc                             >  2:797507/1‑150 (MQ=255)
                                                                                                                         aaaTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAAc                             <  2:463102/150‑1 (MQ=255)
                                                                                                                          aaTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAAcc                            >  2:4555597/1‑150 (MQ=255)
                                                                                                                           aTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCa                           <  1:1755092/150‑1 (MQ=255)
                                                                                                                           aTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCa                           >  2:4763969/1‑150 (MQ=255)
                                                                                                                             ccACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTGAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGt                         <  2:3452789/150‑1 (MQ=255)
                                                                                                                                ccGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAAt                      >  2:1807996/1‑150 (MQ=255)
                                                                                                                                ccGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAAt                      <  2:612778/150‑1 (MQ=255)
                                                                                                                                       aTTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTc               <  1:2151991/150‑1 (MQ=255)
                                                                                                                                       aTTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTc               >  1:4085310/1‑150 (MQ=255)
                                                                                                                                       aTTGCTGATACGATGCGCCCAGGAGAAGGTTTCGTCGATATGGTTCGCCTTTAAATCTTCGACTCTGTTTAGGCCATTGCGCCACAGGCTAATATCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGAGTc               <  2:4657206/150‑1 (MQ=255)
                                                                                                                                        ttGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTcc              >  2:1895977/1‑150 (MQ=255)
                                                                                                                                            tGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACt          >  1:3224822/1‑150 (MQ=255)
                                                                                                                                            tGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACt          >  1:1921354/1‑150 (MQ=255)
                                                                                                                                            tGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACt          >  1:1382606/1‑150 (MQ=255)
                                                                                                                                               tACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcg       >  1:2988642/1‑150 (MQ=255)
                                                                                                                                               tACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcg       <  2:4270298/150‑1 (MQ=255)
                                                                                                                                                aCGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcgg      >  2:3908328/1‑150 (MQ=255)
                                                                                                                                                 cGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcggc     >  1:3923530/1‑150 (MQ=255)
                                                                                                                                                 cGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCcggc     <  1:952020/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGA  >  CP009273/3851126‑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: