Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,710,595 C→A 6.2% intergenic (+29/‑75) yiaD → / → ghrB multicopy suppressor of bamB; outer membrane lipoprotein/glyoxylate/hydroxypyruvate reductase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,710,5950CA 6.2% 815.2 / 19.8 210intergenic (+29/‑75)yiaD/ghrBmulticopy suppressor of bamB; outer membrane lipoprotein/glyoxylate/hydroxypyruvate reductase B
Reads supporting (aligned to +/- strand):  ref base C (84/112);  new base A (5/9);  total (89/121)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

CGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCAC  >  CP009273/3710446‑3710740
                                                                                                                                                     |                                                                                                                                                 
cGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATc                                                                                                                                                   <  2:2791325/150‑1 (MQ=255)
cGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATc                                                                                                                                                   <  2:703108/150‑1 (MQ=255)
  tGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCAc                                                                                                                                                 >  2:3993780/1‑150 (MQ=255)
    gACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCGGGTGATGCAGAAATCACCt                                                                                                                                               <  2:3252861/150‑1 (MQ=255)
      cGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTgg                                                                                                                                             >  1:431884/1‑150 (MQ=255)
            ccGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCCTGCCAGGTGATGCAGAAATCACCTGGCAttct                                                                                                                                       <  2:1014225/150‑1 (MQ=255)
             cGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCAttctt                                                                                                                                      <  1:2070025/150‑1 (MQ=255)
               cATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCa                                                                                                                                    >  1:685744/1‑150 (MQ=255)
                 tCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGTTTATGCAGAAATCACCTGGCATTCTTCAGt                                                                                                                                  <  1:2659781/150‑1 (MQ=255)
                  ccGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGtt                                                                                                                                 <  2:2708473/150‑1 (MQ=255)
                  ccGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGtt                                                                                                                                 <  1:3115907/150‑1 (MQ=255)
                    gTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCCTGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGtttt                                                                                                                               <  2:1562844/150‑1 (MQ=255)
                    gTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGtttt                                                                                                                               >  2:1552475/1‑150 (MQ=255)
                    gTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGtttt                                                                                                                               >  1:424655/1‑150 (MQ=255)
                      aCTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTcc                                                                                                                             >  2:3927096/1‑150 (MQ=255)
                        tCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCtt                                                                                                                           <  2:2885894/150‑1 (MQ=255)
                         cAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTtc                                                                                                                          >  1:285047/1‑150 (MQ=255)
                           gggCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTtcat                                                                                                                        <  1:1924583/150‑1 (MQ=255)
                             gCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTtcatca                                                                                                                      >  1:3102656/1‑150 (MQ=255)
                              ccTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTcatcat                                                                                                                     >  1:953855/1‑150 (MQ=255)
                              ccTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTcatcat                                                                                                                     >  1:1910078/1‑150 (MQ=255)
                               cTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCata                                                                                                                    >  2:1755982/1‑150 (MQ=255)
                                ttGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCatat                                                                                                                   >  2:1863262/1‑150 (MQ=255)
                                ttGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCatat                                                                                                                   <  2:343385/150‑1 (MQ=255)
                                 tGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATAtt                                                                                                                  <  1:482682/150‑1 (MQ=255)
                                  ggCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATAttt                                                                                                                 <  1:3665651/150‑1 (MQ=255)
                                  ggCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATAACCTGGCATTCTTCAGTTTTCCTTCATCATAttt                                                                                                                 >  1:417682/1‑150 (MQ=255)
                                     ccGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAg                                                                                                              <  2:2369809/150‑1 (MQ=255)
                                      cGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAgg                                                                                                             <  1:2099443/150‑1 (MQ=255)
                                           aaCCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAg                                                                                                        <  2:1377572/150‑1 (MQ=255)
                                            aCCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAgg                                                                                                       <  1:3310557/150‑1 (MQ=255)
                                               cAATCGCCAGCAACAGCACCGCAGAAGTTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGa                                                                                                    <  1:3348752/150‑1 (MQ=255)
                                               cAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGa                                                                                                    >  2:913266/1‑150 (MQ=255)
                                                 aTCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATc                                                                                                  >  1:3801584/1‑150 (MQ=255)
                                                   cGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGc                                                                                                >  2:1956664/1‑150 (MQ=255)
                                                   cGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGc                                                                                                <  2:4088416/150‑1 (MQ=255)
                                                    gCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGct                                                                                               <  2:836494/150‑2 (MQ=255)
                                                      cAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCtt                                                                                             <  1:710267/150‑1 (MQ=255)
                                                      cAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCtt                                                                                             <  2:1612297/150‑1 (MQ=255)
                                                        gCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTAt                                                                                           >  1:2512094/1‑150 (MQ=255)
                                                         cAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATAACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATc                                                                                          <  2:1978862/150‑1 (MQ=255)
                                                             aGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTg                                                                                      >  1:934548/1‑150 (MQ=255)
                                                              gCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATAACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGa                                                                                     <  1:2344767/150‑1 (MQ=255)
                                                                aCCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAAt                                                                                   <  1:930121/150‑1 (MQ=255)
                                                                aCCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAAt                                                                                   >  2:240887/1‑150 (MQ=255)
                                                                aCCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAAt                                                                                   <  1:1841630/150‑1 (MQ=255)
                                                                 ccGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATg                                                                                  <  2:1357747/150‑1 (MQ=255)
                                                                   gCAGAAGGTTAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGga                                                                                <  2:2521462/150‑1 (MQ=255)
                                                                     aGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGgaga                                                                              <  2:3497832/150‑1 (MQ=255)
                                                                         ggTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAg                                                                          <  1:1813054/150‑1 (MQ=255)
                                                                           tAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCa                                                                        >  1:1327000/1‑150 (MQ=255)
                                                                            aaGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCAt                                                                       <  1:3652942/150‑1 (MQ=255)
                                                                             aGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATg                                                                      >  2:2095444/1‑150 (MQ=255)
                                                                             aGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATg                                                                      >  1:2266328/1‑150 (MQ=255)
                                                                                 gcAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGc                                                                  <  2:330253/150‑1 (MQ=255)
                                                                                   aaaaCCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTTAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGccg                                                                <  2:4069245/150‑1 (MQ=255)
                                                                                    aaaCCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCGTATCAGTGAATGGAGAGAAGCATGAAGccgt                                                               >  2:1988894/1‑150 (MQ=255)
                                                                                        cgccgTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGccgtccgt                                                           <  1:1468258/150‑1 (MQ=255)
                                                                                        cgccgTGTAGAAAGTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGccgtccgt                                                           <  1:1398158/150‑1 (MQ=255)
                                                                                          ccgTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTa                                                         >  2:1214468/1‑150 (MQ=255)
                                                                                          ccgTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTa                                                         >  2:2073303/1‑150 (MQ=255)
                                                                                             tgtAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATcc                                                      <  2:3528581/150‑1 (MQ=255)
                                                                                             tgtAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATcc                                                      >  2:3724129/1‑150 (MQ=255)
                                                                                               tAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCctc                                                    <  1:1107643/150‑1 (MQ=255)
                                                                                               tAGAAATTACCTTAAGCCCGCTGTAAGCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCctc                                                    <  2:1312195/150‑1 (MQ=255)
                                                                                                agtaaTTACCTTAAGCGCGCTGTAATCGCTTTCATGCCAGGTGATGCAGATATAACCTGGCATTCTTCAGTTTTCCTTCATCTTATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCctct                                                   <  2:3891922/147‑1 (MQ=255)
                                                                                                aGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCctct                                                   >  1:337984/1‑150 (MQ=255)
                                                                                                 gAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTa                                                  <  1:3513954/150‑1 (MQ=255)
                                                                                                 gAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATAACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGTATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTa                                                  <  1:45842/150‑1 (MQ=255)
                                                                                                  aaaTTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTAc                                                 >  1:1252111/1‑150 (MQ=255)
                                                                                                  aaaTTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTAc                                                 >  2:756107/1‑150 (MQ=255)
                                                                                                  aaaTTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTAc                                                 >  1:26765/1‑150 (MQ=255)
                                                                                                    aTTACCTTAAGCCCGCTGTAATCCCTGTCATGCCAGGTGATGCAGAAATCACCTGGCAGTCTTCAGGTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGCATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACaa                                               >  2:3365334/1‑150 (MQ=255)
                                                                                                      tACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAg                                             >  1:4038289/1‑150 (MQ=255)
                                                                                                         cTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCt                                          <  2:3072844/150‑1 (MQ=255)
                                                                                                          ttAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCtt                                         <  1:2161391/150‑1 (MQ=255)
                                                                                                                ccGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCtga                                   >  1:3827032/1‑150 (MQ=255)
                                                                                                                 cGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCtgat                                  >  2:4003208/1‑150 (MQ=255)
                                                                                                                      tAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGAttt                             >  1:1704647/1‑150 (MQ=255)
                                                                                                                        aTCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTAc                           >  2:4170730/1‑150 (MQ=255)
                                                                                                                        aTCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTAc                           <  2:4133064/150‑1 (MQ=255)
                                                                                                                        aTCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTAc                           >  2:3279093/1‑150 (MQ=255)
                                                                                                                         tCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACt                          >  2:1853863/1‑150 (MQ=255)
                                                                                                                          cccTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATAATTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTg                         >  2:356783/1‑150 (MQ=255)
                                                                                                                                 aTGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGcc                  <  1:3913772/150‑1 (MQ=255)
                                                                                                                                 aTGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGcc                  >  1:157586/1‑150 (MQ=255)
                                                                                                                                   gCCAGGTGCTGCAGAAATCACCTGGCATTCTTCCGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCACCGCCTg                <  1:3180834/150‑1 (MQ=255)
                                                                                                                                   gCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTg                <  1:1661476/150‑1 (MQ=255)
                                                                                                                                    ccAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGc               >  2:2938145/1‑150 (MQ=255)
                                                                                                                                     cAGGTGATGCAGAAATCACCTGGCATTCTTAAGTTTTCCTTCATCATCTTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGACGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCa              <  2:218390/150‑1 (MQ=255)
                                                                                                                                      aGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCaa             >  2:489068/1‑150 (MQ=255)
                                                                                                                                      aGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCaa             >  1:2120762/1‑150 (MQ=255)
                                                                                                                                       ggTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAATCCTTACCTGATGATTTACTGCAACGCCTGCAag            <  2:1671580/150‑1 (MQ=255)
                                                                                                                                       ggTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAag            <  2:1671514/150‑1 (MQ=255)
                                                                                                                                        gTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAaga           >  2:67942/1‑150 (MQ=255)
                                                                                                                                           aTGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCa        <  1:3533233/150‑1 (MQ=255)
                                                                                                                                           aTGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCa        >  2:1351161/1‑150 (MQ=255)
                                                                                                                                             gCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCAtt      <  1:1064123/150‑1 (MQ=255)
                                                                                                                                               aGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTc    >  2:960412/1‑150 (MQ=255)
                                                                                                                                                 aaaTCACCTGGCATTCTTCAGTTTTCCTTCACCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCAc  <  1:4099823/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTTCATGCCAGGTGATGCAGAAATCACCTGGCATTCTTCAGTTTTCCTTCATCATATTTCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCAC  >  CP009273/3710446‑3710740

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: