Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,275 G→T 11.4% 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,2750GT11.4% 911.3 / 60.5 254H164N (CAT→AAT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (108/117);  new base T (17/12);  total (125/129)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 263 total aligned reads displayed.

CTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGAC  >  CP009273/3851126‑3851424
                                                                                                                                                     |                                                                                                                                                     
cTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATg                                                                                                                                                       <  2:4100808/150‑1 (MQ=255)
 tCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGACCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgc                                                                                                                                                      >  2:2286778/1‑150 (MQ=255)
   gtagtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATgcgc                                                                                                                                                    >  2:3967099/1‑150 (MQ=255)
    tagtaAAAATCAGAGTATTTCTCTAAATACTCCACCTGACAGGTGATCGGGTGTTGCGGCTCATCATTCGAAACCACCATAATGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTTCTGATACGATGCGcc                                                                                                                                                   <  2:4384517/150‑1 (MQ=255)
     agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGccc                                                                                                                                                  >  2:237378/1‑150 (MQ=255)
     agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGccc                                                                                                                                                  >  2:2551736/1‑150 (MQ=255)
     agtaAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGccc                                                                                                                                                  >  1:2729933/1‑150 (MQ=255)
       taAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAg                                                                                                                                                <  2:2039216/150‑1 (MQ=255)
       taAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAg                                                                                                                                                <  2:3857946/150‑1 (MQ=255)
        aaaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAGCGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAgg                                                                                                                                               >  1:4632077/1‑150 (MQ=255)
         aaaaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGt                                                                                                                                              <  2:688908/150‑1 (MQ=255)
           aaTCAGCGTATTTCTCTAAATACTCCTCCGGACATGTGAACGTGTGATGCGGCTCATCTTACGAAACCAACATCAGGAAAGGCTCGTCTTCGCGCGCAGGCTGTTGCTGAAAATCCACCGCCCGATTGCTTATACGATTCGCCCAGGTGa                                                                                                                                            <  1:3977682/150‑1 (MQ=255)
           aaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGa                                                                                                                                            >  2:3544513/1‑150 (MQ=255)
           aaTCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGa                                                                                                                                            <  1:2806709/150‑1 (MQ=255)
             tCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAg                                                                                                                                          >  1:3497493/1‑150 (MQ=255)
               aGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGt                                                                                                                                        >  1:4085655/1‑150 (MQ=255)
                gCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGtt                                                                                                                                       <  1:1476041/150‑1 (MQ=255)
                 cGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGttt                                                                                                                                      <  2:1091348/150‑1 (MQ=255)
                  gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                     <  1:1095697/150‑1 (MQ=255)
                  gTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtc                                                                                                                                     <  2:48460/150‑1 (MQ=255)
                   tATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcg                                                                                                                                    >  2:2720106/1‑150 (MQ=255)
                   tATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTtcg                                                                                                                                    <  1:1252813/150‑1 (MQ=255)
                       tctcTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGa                                                                                                                                >  2:1803790/1‑150 (MQ=255)
                           tAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATg                                                                                                                            <  2:1088386/150‑1 (MQ=255)
                           tAAATACTCCACCGGACAGGTGAAAGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTGTCGTCGATATg                                                                                                                            >  2:661512/1‑150 (MQ=255)
                            aaaTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATgg                                                                                                                           >  2:1010853/1‑150 (MQ=255)
                             aaTACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGt                                                                                                                          <  2:3414548/150‑1 (MQ=255)
                                 cTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGc                                                                                                                      >  2:1960884/1‑150 (MQ=255)
                                  tCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                     <  2:2517676/150‑1 (MQ=255)
                                  tCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                     >  2:863286/1‑150 (MQ=255)
                                  tCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                     >  2:246082/1‑150 (MQ=255)
                                  tCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGcc                                                                                                                     >  1:3728141/1‑150 (MQ=255)
                                    cACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTg                                                                                                                   <  2:1400783/150‑1 (MQ=255)
                                      ccGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTa                                                                                                                 <  1:2480992/150‑1 (MQ=255)
                                       cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                                <  2:2764899/150‑1 (MQ=255)
                                       cGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaa                                                                                                                >  2:4477976/1‑150 (MQ=255)
                                        ggACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATAGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTaaa                                                                                                               >  2:2977466/1‑150 (MQ=255)
                                           cAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCt                                                                                                            <  1:849827/150‑1 (MQ=255)
                                           cAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCt                                                                                                            <  1:3300747/150‑1 (MQ=255)
                                            aGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCtt                                                                                                           >  2:4194184/1‑150 (MQ=255)
                                             ggTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTc                                                                                                          >  1:1767327/1‑150 (MQ=255)
                                              gTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCg                                                                                                         <  1:1753991/150‑1 (MQ=255)
                                                gAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGAc                                                                                                       >  2:1429003/1‑150 (MQ=255)
                                                gAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGAc                                                                                                       <  1:3605875/150‑1 (MQ=255)
                                                   cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                    <  2:2812086/150‑1 (MQ=255)
                                                   cGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                    <  2:2737261/150‑1 (MQ=255)
                                                   cGGGTGATGCGGATCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                    <  1:166698/150‑1 (MQ=255)
                                                   cGGGTGATGCGGATCATAATACGAAACCACCATCAGGAAAGGCTCGACTGAGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCt                                                                                                    <  1:166564/150‑1 (MQ=255)
                                                    gggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTg                                                                                                   >  2:3275622/1‑150 (MQ=255)
                                                    gggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTg                                                                                                   <  2:2930703/150‑1 (MQ=255)
                                                    gggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTATCGTCGATATGGGTCGCCTGTAAATCTTCGACGCTg                                                                                                   >  2:3317711/1‑150 (MQ=255)
                                                     ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                                  <  1:3083081/150‑1 (MQ=255)
                                                     ggTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGt                                                                                                  <  2:187342/150‑1 (MQ=255)
                                                        gATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGGTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTa                                                                                               <  2:4331616/150‑1 (MQ=255)
                                                            cGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGcc                                                                                           >  2:3400400/1‑150 (MQ=255)
                                                               cTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCAtt                                                                                        >  1:388218/1‑150 (MQ=255)
                                                                  atcatACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTgcg                                                                                     <  1:2868786/150‑1 (MQ=255)
                                                                     atACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCca                                                                                  <  1:2308130/150‑1 (MQ=255)
                                                                      tACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcac                                                                                 >  1:243166/1‑150 (MQ=255)
                                                                       aCGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCcaca                                                                                <  2:3758580/150‑1 (MQ=255)
                                                                         gAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGTCATTGCGCCACAgg                                                                              >  1:4163753/1‑150 (MQ=255)
                                                                         gAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAgg                                                                              <  2:1177209/150‑1 (MQ=255)
                                                                             ccaccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTaa                                                                          <  2:220011/150‑1 (MQ=255)
                                                                             ccaccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTaa                                                                          <  1:1427816/150‑1 (MQ=255)
                                                                              caccaTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAt                                                                         <  1:2204429/150‑1 (MQ=255)
                                                                                 caTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAAtctc                                                                      >  2:4268521/1‑150 (MQ=255)
                                                                                  aTCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATctct                                                                     >  2:4352153/1‑150 (MQ=255)
                                                                                   tCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCtt                                                                    <  1:1960884/150‑1 (MQ=255)
                                                                                   acAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCtt                                                                    >  1:748024/2‑150 (MQ=255)
                                                                                    cAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttt                                                                   >  1:1154351/1‑150 (MQ=255)
                                                                                    cAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCttt                                                                   >  1:800362/1‑150 (MQ=255)
                                                                                       gAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTc                                                                <  2:3222751/150‑1 (MQ=255)
                                                                                        aaaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTcc                                                               >  1:2019400/1‑150 (MQ=255)
                                                                                         aaGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCg                                                              <  1:797273/150‑1 (MQ=255)
                                                                                          aGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGt                                                             >  1:855277/1‑150 (MQ=255)
                                                                                             cTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAg                                                          <  1:1884817/150‑1 (MQ=255)
                                                                                                gTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTc                                                       >  2:3240266/1‑150 (MQ=255)
                                                                                                  cTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGc                                                     >  2:4555194/1‑150 (MQ=255)
                                                                                                   tGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCt                                                    >  2:1234086/1‑150 (MQ=255)
                                                                                                   tGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCt                                                    <  1:1204995/150‑1 (MQ=255)
                                                                                                    gcgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTa                                                   <  2:1331308/150‑1 (MQ=255)
                                                                                                     cgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTaa                                                  <  1:4060109/150‑1 (MQ=255)
                                                                                                     cgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTaa                                                  >  2:495147/1‑150 (MQ=255)
                                                                                                     cgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTaa                                                  >  2:1100971/1‑150 (MQ=255)
                                                                                                     cgcgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCAAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCTTTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTaa                                                  <  1:358446/150‑1 (MQ=255)
                                                                                                       cgcgcAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGa                                                >  2:3337177/1‑150 (MQ=255)
                                                                                                           cAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGt                                            >  1:3847742/1‑150 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                           >  2:1558867/1‑150 (MQ=255)
                                                                                                            aGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGtt                                           <  1:4263634/150‑1 (MQ=255)
                                                                                                               cTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGc                                        >  1:443942/1‑150 (MQ=255)
                                                                                                                 gTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCTATATGGTTCGCCTGTAAATCTTCGACGATGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCAGTCAGTTCGCTAAGATAGTTCGccc                                      >  1:4179222/1‑150 (MQ=255)
                                                                                                                     cAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATc                                  <  1:3912204/150‑1 (MQ=255)
                                                                                                                      aGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCg                                 >  1:3161168/1‑150 (MQ=255)
                                                                                                                           aTCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCa                            >  1:2374029/1‑150 (MQ=255)
                                                                                                                              cACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTa                         >  2:406074/1‑150 (MQ=255)
                                                                                                                                     cGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCg                  >  1:956189/1‑150 (MQ=255)
                                                                                                                                        ttGCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTcc               >  2:3819143/1‑150 (MQ=255)
                                                                                                                                          gCTGATACGATTCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCa             <  1:4489748/150‑1 (MQ=255)
                                                                                                                                          gCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCa             <  2:375464/150‑1 (MQ=255)
                                                                                                                                              aTACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTcc         >  2:3709297/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCGTAGTAAAAATCAGCGTATTTCTCTAAATACTCCACCGGACAGGTGAACGGGTGATGCGGCTCATCATACGAAACCACCATCAGGAAAGGCTCGTCTGCGCGCGCAGGCTGTTGCAGAAAATCCACCGCCCGATTGCTGATACGATGCGCCCAGGTGAAGGTTTCGTCGATATGGTTCGCCTGTAAATCTTCGACGCTGTTTAGGCCATTGCGCCACAGGCTAATCTCTTTTTCCGTCAGTTCGCTAAGATAGTTCGCCCCATCGAACCAGTAATCAGCGTCCCACTCCGGCGGAC  >  CP009273/3851126‑3851424

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: ‑