Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,969,710 G→T 100% G25C (GGT→TGT)  wzxE → O‑antigen translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,969,7100GT97.5% 941.9 / 13.1 283G25C (GGT→TGT) wzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (150/126);  total (152/131)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 283 total aligned reads displayed.

TGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCC  >  CP009273/3969563‑3969855
                                                                                                                                                   |                                                                                                                                                 
   cGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAg                                                                                                                                              <  2:1713204/150‑1 (MQ=255)
       gTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAATAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGc                                                                                                                                          >  2:2472864/1‑150 (MQ=255)
         tCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTg                                                                                                                                        <  1:4064322/150‑1 (MQ=255)
            aCAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCg                                                                                                                                     >  2:2086016/1‑150 (MQ=255)
            aCAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCg                                                                                                                                     >  2:4416734/1‑150 (MQ=255)
                ccTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGgtgt                                                                                                                                 <  2:4491720/150‑1 (MQ=255)
                  tGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCa                                                                                                                               >  2:3390761/1‑150 (MQ=255)
                   gTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAt                                                                                                                              <  2:3526776/150‑1 (MQ=255)
                    tCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAtt                                                                                                                             <  1:540258/150‑1 (MQ=255)
                     cGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAttt                                                                                                                            >  2:4309332/1‑150 (MQ=255)
                       cccGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTgg                                                                                                                          >  2:3561705/1‑150 (MQ=255)
                        ccGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGt                                                                                                                         >  2:512479/1‑150 (MQ=255)
                             aaTCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcggc                                                                                                                    >  2:3925142/1‑150 (MQ=255)
                             aaTCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcggc                                                                                                                    <  2:3430461/150‑1 (MQ=255)
                             aaTCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGACCGTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGATGGCGGTGTCAGTTGGTCaggc                                                                                                                    >  1:4301756/1‑150 (MQ=255)
                              aTCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcggcg                                                                                                                   >  1:3073722/1‑150 (MQ=255)
                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                 <  2:2684412/150‑1 (MQ=255)
                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                 <  1:1044767/150‑1 (MQ=255)
                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                 <  1:1396511/150‑1 (MQ=255)
                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                 <  1:3207689/150‑1 (MQ=255)
                                   cGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCt                                                                                                              <  1:3340583/150‑1 (MQ=255)
                                    gTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCtt                                                                                                             <  1:3495600/150‑1 (MQ=255)
                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                          <  1:2370334/150‑1 (MQ=255)
                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                          <  1:4309901/150‑1 (MQ=255)
                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                          >  1:64200/1‑150 (MQ=255)
                                          tAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTg                                                                                                       <  1:1229466/150‑1 (MQ=255)
                                          tAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTg                                                                                                       >  1:3073457/1‑150 (MQ=255)
                                            aTTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggc                                                                                                     <  1:2290902/150‑1 (MQ=255)
                                            aTTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggc                                                                                                     <  1:3678808/150‑1 (MQ=255)
                                                cGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                 >  1:184876/1‑150 (MQ=255)
                                                 gACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCaa                                                                                                <  2:4630726/150‑1 (MQ=255)
                                                  aCTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCaaa                                                                                               <  1:741699/150‑1 (MQ=255)
                                                  aCTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCaaa                                                                                               >  2:4641510/1‑150 (MQ=255)
                                                    tttGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAtt                                                                                             >  2:1409805/1‑150 (MQ=255)
                                                     ttgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAttt                                                                                            <  2:1822236/150‑1 (MQ=255)
                                                      tgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTc                                                                                           >  2:1281258/1‑150 (MQ=255)
                                                       gttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTcc                                                                                          >  1:283234/1‑150 (MQ=255)
                                                       gttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTcc                                                                                          <  2:2573459/150‑1 (MQ=255)
                                                         tgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGc                                                                                        >  1:4116736/1‑150 (MQ=255)
                                                          gAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGcg                                                                                       <  2:3344216/150‑2 (MQ=255)
                                                            actactTTTCCTGATATGTCTTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGATTGGGCTGGCGGCAAATTTCCGCCAg                                                                                     <  1:3070222/150‑1 (MQ=255)
                                                             ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGt                                                                                    >  2:3911573/1‑150 (MQ=255)
                                                               actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTg                                                                                  >  1:4185770/1‑150 (MQ=255)
                                                               actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTg                                                                                  >  2:1094348/1‑150 (MQ=255)
                                                               actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTg                                                                                  >  1:1830598/1‑150 (MQ=255)
                                                                ctTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGa                                                                                 >  1:3692131/1‑150 (MQ=255)
                                                                  tttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAtt                                                                               <  1:193520/150‑1 (MQ=255)
                                                                   ttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTa                                                                              >  1:3297015/1‑150 (MQ=255)
                                                                     ccTGATATGTCGTTGGCAAAAGCGTCCATGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                            >  2:3444353/1‑150 (MQ=255)
                                                                       tGATATGTCGTTGGCAAAAGCGTCCTTGTGGATGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGt                                                                          >  2:3310080/1‑150 (MQ=255)
                                                                         atatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGc                                                                        >  2:1600823/1‑150 (MQ=255)
                                                                          tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCt                                                                       <  2:3159561/150‑1 (MQ=255)
                                                                          tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCt                                                                       >  2:2914816/1‑150 (MQ=255)
                                                                           atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                      <  2:1938053/150‑1 (MQ=255)
                                                                           atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                      <  1:2539017/150‑1 (MQ=255)
                                                                           atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                      >  2:3722374/1‑150 (MQ=255)
                                                                            tGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCg                                                                     <  1:3145690/150‑1 (MQ=255)
                                                                                gTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGt                                                                 >  2:3946330/1‑150 (MQ=255)
                                                                                    gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                             <  1:1743572/150‑1 (MQ=255)
                                                                                     cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                            <  2:2987660/150‑1 (MQ=255)
                                                                                      aaaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGc                                                           >  2:2399544/1‑150 (MQ=255)
                                                                                      aaaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGc                                                           <  1:2465129/150‑1 (MQ=255)
                                                                                         aGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                                        <  1:1193494/150‑1 (MQ=255)
                                                                                         aGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                                        >  1:2801366/1‑150 (MQ=255)
                                                                                          gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                       <  2:1025872/150‑1 (MQ=255)
                                                                                           cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                      >  1:3765190/1‑150 (MQ=255)
                                                                                            gTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                     >  2:899442/1‑150 (MQ=255)
                                                                                            gTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                     >  2:227897/1‑150 (MQ=255)
                                                                                            gTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                     >  1:168413/1‑150 (MQ=255)
                                                                                             tCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                    <  1:2597803/150‑1 (MQ=255)
                                                                                                  gtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                               <  1:2008838/150‑1 (MQ=255)
                                                                                                      aCGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                           <  1:3040727/150‑1 (MQ=255)
                                                                                                         gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                        <  1:2195240/150‑1 (MQ=255)
                                                                                                         gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                        <  2:1082364/150‑1 (MQ=255)
                                                                                                         gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                        <  1:2674334/150‑1 (MQ=255)
                                                                                                         gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                        <  1:4528176/150‑1 (MQ=255)
                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                  >  2:3272736/1‑150 (MQ=255)
                                                                                                                   cacTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                              <  1:1114728/150‑1 (MQ=255)
                                                                                                                   cacTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                              <  2:3967235/150‑1 (MQ=255)
                                                                                                                   cacTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                              >  1:2226511/1‑150 (MQ=255)
                                                                                                                      tGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                           >  1:2385759/1‑150 (MQ=255)
                                                                                                                       ggTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                          >  1:908580/1‑150 (MQ=255)
                                                                                                                         tCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACg                        >  2:3547856/1‑150 (MQ=255)
                                                                                                                           aaGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt                      <  2:1191969/150‑1 (MQ=255)
                                                                                                                           aaGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt                      <  2:4038235/150‑1 (MQ=255)
                                                                                                                           aaGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTATATCGGTGTAACCAAATACGtt                      >  2:3813564/1‑150 (MQ=255)
                                                                                                                             gATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                    <  2:3189795/150‑1 (MQ=255)
                                                                                                                              aTTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGcc                   >  2:1497423/1‑150 (MQ=255)
                                                                                                                                tGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCa                 >  1:450525/1‑150 (MQ=255)
                                                                                                                                 ggTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg                >  1:1598577/1‑150 (MQ=255)
                                                                                                                                    gCCGGGTTACTGGTCTGTAGGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc             >  1:1959153/1‑150 (MQ=255)
                                                                                                                                    gCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc             <  1:4418536/150‑1 (MQ=255)
                                                                                                                                      cGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCa           >  2:2945528/1‑150 (MQ=255)
                                                                                                                                      cGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCa           >  1:2157426/1‑150 (MQ=255)
                                                                                                                                       gggTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCAt          >  2:2879360/1‑150 (MQ=255)
                                                                                                                                          ttACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt       <  2:799344/150‑1 (MQ=255)
                                                                                                                                           tACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATa      >  1:3588362/1‑150 (MQ=255)
                                                                                                                                            aCTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATaa     >  2:1585854/1‑150 (MQ=255)
                                                                                                                                              tGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATc   <  2:3614126/150‑1 (MQ=255)
                                                                                                                                               ggTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc  <  2:1957932/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCC  >  CP009273/3969563‑3969855

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: