Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,037,576 G→T 5.0% G697G (GGC→GGA)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,037,5760GT 5.0% 826.6 / 15.8 217G697G (GGC→GGA) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (113/93);  new base T (4/7);  total (117/100)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 227 total aligned reads displayed.

GAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCG  >  CP009273/1037427‑1037724
                                                                                                                                                     |                                                                                                                                                    
taaCTGTCTAATTCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAgcg                                                                                                                                                      <  2:4191303/149‑1 (MQ=255)
gaaCTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAgcg                                                                                                                                                      <  2:4191946/150‑1 (MQ=255)
 aaCTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAgcgc                                                                                                                                                     >  1:1079543/1‑150 (MQ=255)
   cTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCt                                                                                                                                                   <  1:1661468/150‑1 (MQ=255)
    tGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCtt                                                                                                                                                  <  1:1424182/150‑1 (MQ=255)
        tAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAat                                                                                                                                              <  2:2530452/150‑1 (MQ=255)
        tAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAat                                                                                                                                              <  2:4484954/150‑1 (MQ=255)
        tAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAat                                                                                                                                              >  1:2752378/1‑150 (MQ=255)
        tAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAat                                                                                                                                              <  2:2530501/150‑1 (MQ=255)
           tGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATAtt                                                                                                                                           >  1:70846/1‑150 (MQ=255)
           tGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCGTTAATAtt                                                                                                                                           >  1:3954659/1‑150 (MQ=255)
             cGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGa                                                                                                                                         <  2:3865062/150‑1 (MQ=255)
             cGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGa                                                                                                                                         <  2:2988355/150‑1 (MQ=255)
              gTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGAc                                                                                                                                        <  1:1855330/150‑1 (MQ=255)
               tAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATTGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACg                                                                                                                                       <  2:1514213/150‑1 (MQ=255)
                aGACCGGAAAAGGCGTTCACGCCGCCTCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgc                                                                                                                                      <  2:4012321/150‑1 (MQ=255)
                aGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgc                                                                                                                                      >  1:2944898/1‑150 (MQ=255)
                aGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgc                                                                                                                                      <  2:4530418/150‑1 (MQ=255)
                aGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgc                                                                                                                                      >  1:2708885/1‑150 (MQ=255)
                aGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgc                                                                                                                                      >  1:2707643/1‑150 (MQ=255)
                        aaaGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGt                                                                                                                              >  1:3554444/1‑150 (MQ=255)
                          aGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTc                                                                                                                            <  1:2990588/150‑1 (MQ=255)
                          aGGCGGTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTc                                                                                                                            <  2:3505091/150‑1 (MQ=255)
                              gTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                                                                                        >  1:2781912/1‑150 (MQ=255)
                              gTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                                                                                        >  1:2097601/1‑150 (MQ=255)
                               ttCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGAc                                                                                                                       <  1:2554153/150‑1 (MQ=255)
                               ttCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGAc                                                                                                                       <  2:1349452/150‑1 (MQ=255)
                               ttCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGAc                                                                                                                       <  1:4109323/150‑1 (MQ=255)
                                tCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACg                                                                                                                      >  2:2354218/1‑150 (MQ=255)
                                tCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACg                                                                                                                      >  2:2353918/1‑150 (MQ=255)
                                 cACGCCGCATCCGGCCACTTTCAGGTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGc                                                                                                                     >  2:1909936/1‑150 (MQ=255)
                                     ccgcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCa                                                                                                                 <  1:1041347/150‑1 (MQ=255)
                                       gcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                                                                               >  1:2721463/1‑150 (MQ=255)
                                       gcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                                                                               <  1:743820/150‑1 (MQ=255)
                                        cATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTg                                                                                                              >  2:921774/1‑150 (MQ=255)
                                        cATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTg                                                                                                              >  1:2888658/1‑150 (MQ=255)
                                         aTCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGa                                                                                                             >  1:659682/1‑150 (MQ=255)
                                           ccGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACa                                                                                                           >  2:2172422/1‑150 (MQ=255)
                                           ccGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACa                                                                                                           <  2:3304125/150‑1 (MQ=255)
                                              gCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAAc                                                                                                        >  2:2933226/1‑150 (MQ=255)
                                              gCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAAc                                                                                                        >  1:697299/1‑150 (MQ=255)
                                                cACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTc                                                                                                      >  1:578445/1‑150 (MQ=255)
                                                cACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTc                                                                                                      >  1:578717/1‑150 (MQ=255)
                                                  cTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCAc                                                                                                    <  2:4338304/150‑1 (MQ=255)
                                                    ttCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACct                                                                                                  <  2:1919786/150‑1 (MQ=255)
                                                        gTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCttt                                                                                              >  2:3411028/1‑150 (MQ=255)
                                                         atttACTCTATCACGGAGTAACTATAACAGTAATAGATATAGCCGAAACTGTAAGCTTAGCTGCCGAGATTAATCACACCATTGAGGATAGCGCCTAAAATATTAACTCATGCCTGTTCCAGACGCAGCATTGACAAAATCACCTCTATg                                                                                             >  2:388462/2‑150 (MQ=255)
                                                          tttACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTgg                                                                                            >  1:2360698/1‑150 (MQ=255)
                                                              ctctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgt                                                                                        >  2:3059052/1‑150 (MQ=255)
                                                               tctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtg                                                                                       >  2:4478813/1‑150 (MQ=255)
                                                               tctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATCGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtg                                                                                       <  1:1296248/150‑1 (MQ=255)
                                                                ctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtgt                                                                                      >  2:2463797/1‑150 (MQ=255)
                                                                ctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtgt                                                                                      >  2:483494/1‑150 (MQ=255)
                                                                 tttCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGtt                                                                                     <  1:2423612/150‑1 (MQ=255)
                                                                       ggAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGcc                                                                               >  1:1375056/1‑150 (MQ=255)
                                                                       ggAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGcc                                                                               <  2:4539809/150‑1 (MQ=255)
                                                                        gAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCa                                                                              <  1:1710966/150‑1 (MQ=255)
                                                                              cTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc                                                                        >  1:3602441/1‑150 (MQ=255)
                                                                              cTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc                                                                        >  1:438142/1‑150 (MQ=255)
                                                                              cTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCATTTAAAATTCACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgc                                                                        >  1:3696295/1‑150 (MQ=255)
                                                                               tataACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCCTTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcgcg                                                                       <  2:4320896/150‑1 (MQ=255)
                                                                                  aaCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCaa                                                                    >  1:4402369/1‑150 (MQ=255)
                                                                                    ccGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAAcc                                                                  >  2:3126402/1‑150 (MQ=255)
                                                                                      gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcag                                                                >  2:3678485/1‑150 (MQ=255)
                                                                                      gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcag                                                                <  2:2856301/150‑1 (MQ=255)
                                                                                       tAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagc                                                               >  1:2960287/1‑150 (MQ=255)
                                                                                        aaTAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagca                                                              >  2:2378682/1‑150 (MQ=255)
                                                                                             ttATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGGTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTg                                                         <  1:2237164/150‑1 (MQ=255)
                                                                                               ataGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTGGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGt                                                       <  2:1808421/150‑1 (MQ=255)
                                                                                                 aGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGc                                                     >  1:3131080/1‑150 (MQ=255)
                                                                                                 aGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGc                                                     <  1:779865/150‑1 (MQ=255)
                                                                                                  gCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGcc                                                    <  1:3328759/150‑1 (MQ=255)
                                                                                                      tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                >  2:1468288/1‑150 (MQ=255)
                                                                                                      tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                <  1:2676413/150‑1 (MQ=255)
                                                                                                       aaCTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAg                                               >  2:1434595/1‑150 (MQ=255)
                                                                                                          tGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAAc                                            <  2:2433733/150‑1 (MQ=255)
                                                                                                          tGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAAc                                            >  1:2414439/1‑150 (MQ=255)
                                                                                                              aGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGccc                                        >  2:1972998/1‑150 (MQ=255)
                                                                                                                  gTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACg                                    <  1:3226414/150‑1 (MQ=255)
                                                                                                                      tGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACcg                                >  1:199168/1‑150 (MQ=255)
                                                                                                                      tGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACcg                                <  2:653541/150‑1 (MQ=255)
                                                                                                                        gcgcGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACcgcg                              <  1:4479372/150‑1 (MQ=255)
                                                                                                                          gcgTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGc                            >  2:238469/1‑150 (MQ=255)
                                                                                                                          gcgTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGc                            <  2:4121812/150‑1 (MQ=255)
                                                                                                                           cgTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCa                           <  2:2390074/150‑1 (MQ=255)
                                                                                                                           cgTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCa                           <  2:2390054/150‑1 (MQ=255)
                                                                                                                             tttAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCAtc                         >  1:67255/1‑150 (MQ=255)
                                                                                                                                    acacCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCAcc                  >  2:4263590/1‑150 (MQ=255)
                                                                                                                                     cacCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCg                 <  1:3618788/150‑1 (MQ=255)
                                                                                                                                     cacCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCg                 >  1:795160/1‑150 (MQ=255)
                                                                                                                                      acCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGc                >  2:2070014/1‑150 (MQ=255)
                                                                                                                                       ccATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGcc               >  1:4454892/1‑150 (MQ=255)
                                                                                                                                        cATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCa              <  2:1416284/150‑1 (MQ=255)
                                                                                                                                        cATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCa              >  1:1567759/1‑150 (MQ=255)
                                                                                                                                          ttGAGGATAGCGCCTTTAATATTGATGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGc            <  1:3545803/150‑1 (MQ=255)
                                                                                                                                          ttGAGGATAGCGCCTTTAATATTGAGGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTGTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGc            <  1:3547116/150‑1 (MQ=255)
                                                                                                                                            gAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCAt          <  2:2149237/150‑1 (MQ=255)
                                                                                                                                             aGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATc         >  2:918320/1‑150 (MQ=255)
                                                                                                                                             aGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATc         >  1:3911001/1‑150 (MQ=255)
                                                                                                                                                    cgcCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGgcggcg  <  2:2988045/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCG  >  CP009273/1037427‑1037724

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