Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,037,576 G→T 7.4% 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 7.4% 888.5 / 29.0 235G697G (GGC→GGA) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (122/95);  new base T (12/6);  total (134/101)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

AACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCGT  >  CP009273/1037428‑1037725
                                                                                                                                                    |                                                                                                                                                     
     tCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTa                                                                                                                                                 <  2:2750006/150‑1 (MQ=255)
     tCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTa                                                                                                                                                 <  1:230423/150‑1 (MQ=255)
        aaTGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAata                                                                                                                                              >  2:733531/1‑150 (MQ=255)
             gTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGAc                                                                                                                                         <  1:3183925/150‑1 (MQ=255)
              tAGACCGGAAAAGTCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTATTATTGTCg                                                                                                                                        >  2:4258336/1‑150 (MQ=255)
              tAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTTGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACg                                                                                                                                        >  2:1576909/1‑150 (MQ=255)
              tAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTTGGAGTAACTATAACCGTAATAGGTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCAGTGAGGATAGCGCCTTTAATATTGACg                                                                                                                                        >  2:1531437/1‑150 (MQ=255)
              tAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACg                                                                                                                                        >  1:2391711/1‑150 (MQ=255)
                 aCCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcct                                                                                                                                     >  2:1632016/1‑150 (MQ=255)
                 aCCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcct                                                                                                                                     >  1:3500860/1‑150 (MQ=255)
                    ggAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcctgcc                                                                                                                                  >  1:4141994/1‑150 (MQ=255)
                     gAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACgcctgcct                                                                                                                                 <  2:4023168/150‑1 (MQ=255)
                      aaaaGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTATAATATTGACGcctgcctg                                                                                                                                >  2:657662/1‑150 (MQ=255)
                       aaaGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGt                                                                                                                               <  2:2539018/150‑1 (MQ=255)
                        aaGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGtt                                                                                                                              >  1:1658129/1‑150 (MQ=255)
                        aaGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGtt                                                                                                                              >  2:749353/1‑150 (MQ=255)
                         aGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTc                                                                                                                             >  2:580944/1‑150 (MQ=255)
                          ggCGTTCTCGCCGCATAGGGCCACTTTCAGTTTTACTCTTTCTCGGAGTATATATTACCGTATTAGTTATAGACGTAACTGTAAGCGGTGCTGGCGCGTTTAATCAATCCATTGAGGTTTGTGACTTTAATATTGACGCCTGCATGTTcc                                                                                                                            <  2:2249087/150‑1 (MQ=255)
                           gCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCa                                                                                                                           <  1:4015759/150‑1 (MQ=255)
                           gCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCa                                                                                                                           <  1:1329942/150‑1 (MQ=255)
                            cGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAg                                                                                                                          <  2:89064/150‑1 (MQ=255)
                             gTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                                                                                         >  2:2284816/1‑150 (MQ=255)
                             gTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGa                                                                                                                         <  2:2454607/150‑1 (MQ=255)
                               tCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACg                                                                                                                       >  2:391762/1‑150 (MQ=255)
                                  cgccgcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTg                                                                                                                    >  1:1722402/1‑150 (MQ=255)
                                      gcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                                                                                <  2:3678361/150‑1 (MQ=255)
                                      gcATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCAtt                                                                                                                >  1:2550815/1‑150 (MQ=255)
                                        aTCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGa                                                                                                              >  1:519061/1‑150 (MQ=255)
                                         tCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGAc                                                                                                             >  1:4258025/1‑150 (MQ=255)
                                          ccGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACa                                                                                                            <  1:1522160/150‑1 (MQ=255)
                                          ccGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACa                                                                                                            <  1:1522384/150‑1 (MQ=255)
                                           cGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCGGTTCCAGACGCTGCATTGACaa                                                                                                           >  1:3539221/1‑150 (MQ=255)
                                            ggCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACaaa                                                                                                          >  1:2449306/1‑150 (MQ=255)
                                             gCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAAc                                                                                                         <  2:4605379/150‑1 (MQ=255)
                                              ccACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACt                                                                                                        <  1:2452293/150‑1 (MQ=255)
                                              ccACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACt                                                                                                        <  1:778647/150‑1 (MQ=255)
                                                aCTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCa                                                                                                      >  2:2342010/1‑150 (MQ=255)
                                                   ttCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACct                                                                                                   <  2:3107346/150‑1 (MQ=255)
                                                      aGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCtt                                                                                                >  1:681114/1‑150 (MQ=255)
                                                       gTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCttt                                                                                               <  1:1265610/150‑1 (MQ=255)
                                                          ttACTCTTTCTCGGTGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGc                                                                                            <  2:2078491/150‑1 (MQ=255)
                                                            aCTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCAGTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCgg                                                                                          >  2:2841462/1‑150 (MQ=255)
                                                               ctTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGgtgt                                                                                       >  1:2826042/1‑150 (MQ=255)
                                                                tttCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGtt                                                                                      <  1:1785317/150‑1 (MQ=255)
                                                                  tctcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCa                                                                                    >  1:1933347/1‑150 (MQ=255)
                                                                    tcGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAg                                                                                  >  1:3319740/1‑150 (MQ=255)
                                                                       gAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCa                                                                               <  2:3586083/150‑1 (MQ=255)
                                                                        aGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaa                                                                              <  2:4312904/150‑1 (MQ=255)
                                                                         gTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACTCCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCaaa                                                                             >  2:1158081/1‑150 (MQ=255)
                                                                           aaCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAc                                                                           >  1:1097816/1‑150 (MQ=255)
                                                                           aaCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAc                                                                           <  1:1080112/150‑1 (MQ=255)
                                                                            aCTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAAcg                                                                          >  2:3535000/1‑150 (MQ=255)
                                                                               ataACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACgcgc                                                                       >  1:1681039/1‑150 (MQ=255)
                                                                                 aaCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCaa                                                                     >  2:1763571/1‑150 (MQ=255)
                                                                                  aCCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAAc                                                                    >  2:345569/1‑150 (MQ=255)
                                                                                    cGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACca                                                                  >  1:991269/1‑150 (MQ=255)
                                                                                    cGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACca                                                                  <  2:2213846/150‑1 (MQ=255)
                                                                                     gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcag                                                                 <  2:839426/150‑1 (MQ=255)
                                                                                     gTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcag                                                                 >  1:2226987/1‑150 (MQ=255)
                                                                                       aaTAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagca                                                               >  2:4584741/1‑150 (MQ=255)
                                                                                       aaTAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagca                                                               >  2:1873633/1‑150 (MQ=255)
                                                                                        aTAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACcagcag                                                              <  2:1867675/150‑1 (MQ=255)
                                                                                            ttATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTg                                                          >  2:1938426/1‑150 (MQ=255)
                                                                                             tataGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTgg                                                         <  2:1654938/150‑1 (MQ=255)
                                                                                               taGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTg                                                       >  2:1336681/1‑150 (MQ=255)
                                                                                                 gCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGcc                                                     <  2:504895/150‑1 (MQ=255)
                                                                                                  ccGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCa                                                    >  2:2947766/1‑150 (MQ=255)
                                                                                                   cGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCCAAACGCGCAACCAGCAGGCTGGTGCCaa                                                   >  2:2595527/1‑150 (MQ=255)
                                                                                                     tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                 >  2:2705034/1‑150 (MQ=255)
                                                                                                     tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                 >  2:1540642/1‑150 (MQ=255)
                                                                                                     tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                 >  2:1214625/1‑150 (MQ=255)
                                                                                                     tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                 >  1:1230281/1‑150 (MQ=255)
                                                                                                     tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                 <  2:698978/150‑1 (MQ=255)
                                                                                                     tAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACa                                                 >  1:866398/1‑150 (MQ=255)
                                                                                                      aaCTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAg                                                >  1:468502/1‑150 (MQ=255)
                                                                                                      aaCTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAg                                                <  1:3979166/150‑1 (MQ=255)
                                                                                                         tGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAAc                                             <  1:2138333/150‑1 (MQ=255)
                                                                                                          gTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACg                                            <  2:413738/150‑1 (MQ=255)
                                                                                                               cGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCa                                       <  2:371071/150‑1 (MQ=255)
                                                                                                                   gCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGAc                                   <  2:2644804/150‑1 (MQ=255)
                                                                                                                    cTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGAcc                                  <  1:2715011/150‑1 (MQ=255)
                                                                                                                     tGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACcg                                 <  2:4437881/150‑1 (MQ=255)
                                                                                                                       gcgcGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCTTTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACcgcg                               <  1:3094493/150‑1 (MQ=255)
                                                                                                                          cgTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCa                            >  1:54370/1‑150 (MQ=255)
                                                                                                                          cgTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCa                            >  1:3983527/1‑150 (MQ=255)
                                                                                                                                aTCACACCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCAtcactc                      <  1:2064834/150‑1 (MQ=255)
                                                                                                                                     acCATTGAGGATAGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGc                 >  1:1806230/1‑150 (MQ=255)
                                                                                                                                     acCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGc                 >  2:1790143/1‑150 (MQ=255)
                                                                                                                                      ccATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGcc                >  2:3689561/1‑150 (MQ=255)
                                                                                                                                      caaTTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTACAGACGCAGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAAACAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGcc                >  1:2537504/3‑150 (MQ=255)
                                                                                                                                       cATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCa               <  1:4541477/150‑1 (MQ=255)
                                                                                                                                       cATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCa               >  2:3247180/1‑150 (MQ=255)
                                                                                                                                        aTTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAg              >  1:2941497/1‑150 (MQ=255)
                                                                                                                                        aTTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAg              >  1:4424839/1‑150 (MQ=255)
                                                                                                                                             ggATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGAGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATcg         >  2:1513030/1‑150 (MQ=255)
                                                                                                                                                tAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATcggcg      <  1:3203001/150‑1 (MQ=255)
                                                                                                                                                tAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATcggcg      <  2:1286584/150‑1 (MQ=255)
                                                                                                                                                 aGCTCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATcggcgg     >  1:3817268/1‑150 (MQ=255)
                                                                                                                                                  gcgcCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCggcggc    <  2:652552/150‑1 (MQ=255)
                                                                                                                                                    gcCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCGt  >  2:3658882/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACTGTCTAATGCGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCCACTTTCAGTTTTACTCTTTCTCGGAGTAACTATAACCGTAATAGTTATAGCCGTAACTGTAAGCGGTGCTGGCGCGTTTAATCACACCATTGAGGATAGCGCCTTTAATATTGACGCCTGCCTGTTCCAGACGCTGCATTGACAAACTCACCTCTTTGGCGGTGTTCAAGCCAAAACGCGCAACCAGCAGGCTGGTGCCAACAGAACGCCCCACGACCGCGGCATCACTCACCGCCAGCATCGGCGGCGT  >  CP009273/1037428‑1037725

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