Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,424,234 C→A 100% S310* (TCA→TAA)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,424,2340CA97.4% 1021.8 / 15.3 307S310* (TCA→TAA) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base A (162/137);  total (165/142)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 313 total aligned reads displayed.

TTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCGCAG  >  CP009273/1424088‑1424377
                                                                                                                                                  |                                                                                                                                               
ttCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTaaac                                                                                                                                              <  2:3094467/150‑1 (MQ=255)
  ccTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTaaacaa                                                                                                                                            >  1:2720297/1‑150 (MQ=255)
    tCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGt                                                                                                                                          <  2:922626/150‑1 (MQ=255)
        cAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCg                                                                                                                                      <  1:2679962/150‑1 (MQ=255)
        cAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCg                                                                                                                                      >  1:2004247/1‑150 (MQ=255)
         aaCGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCgg                                                                                                                                     <  1:2910543/150‑1 (MQ=255)
          aCGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCggg                                                                                                                                    >  2:2988711/1‑150 (MQ=255)
            ggcggcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCa                                                                                                                                  >  1:3573588/1‑150 (MQ=255)
               ggcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGc                                                                                                                               >  1:4467556/1‑150 (MQ=255)
                gcAGGAAATTCCGCGAAGGCGGCAACCACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGcc                                                                                                                              <  2:4188375/150‑1 (MQ=255)
                gcAGGAAATTCCGCGAAGGCGGCAAAAACGTCAGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAATCAAGTGCCGGGCAGGcc                                                                                                                              <  2:2857723/150‑1 (MQ=255)
                   ggAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCa                                                                                                                           >  2:1710632/1‑150 (MQ=255)
                   ggAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCa                                                                                                                           <  1:3127485/150‑1 (MQ=255)
                   ggAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCa                                                                                                                           <  1:1844043/150‑1 (MQ=255)
                     aaaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGc                                                                                                                         >  2:4267668/1‑150 (MQ=255)
                     aaaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGc                                                                                                                         >  2:2120160/1‑150 (MQ=255)
                     aaaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGc                                                                                                                         >  2:3740317/1‑150 (MQ=255)
                      aaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGcc                                                                                                                        <  2:4578781/150‑1 (MQ=255)
                          ccGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTg                                                                                                                    <  2:1276578/150‑1 (MQ=255)
                            gcgAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCACAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGcc                                                                                                                  <  2:443694/150‑1 (MQ=255)
                                        aaaaCGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGa                                                                                                      >  2:1203417/1‑150 (MQ=255)
                                         aaaCGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGaa                                                                                                     >  1:3743179/1‑150 (MQ=255)
                                          aaCGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGaaa                                                                                                    >  1:3017263/1‑150 (MQ=255)
                                            cGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAAt                                                                                                  <  2:1800632/150‑1 (MQ=255)
                                            cGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAAt                                                                                                  >  1:476504/1‑150 (MQ=255)
                                             gTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATc                                                                                                 <  2:1876538/150‑1 (MQ=255)
                                                cGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGc                                                                                              >  1:1572101/1‑150 (MQ=255)
                                                 gagaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCa                                                                                             <  2:3877064/150‑1 (MQ=255)
                                                   gaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGa                                                                                           <  2:3670757/150‑1 (MQ=255)
                                                   gaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGa                                                                                           <  1:3265058/150‑1 (MQ=255)
                                                     cGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGaaa                                                                                         >  2:2606400/1‑150 (MQ=255)
                                                      gAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAg                                                                                        >  2:1452131/1‑150 (MQ=255)
                                                      gAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAg                                                                                        >  1:2080743/1‑150 (MQ=255)
                                                        aCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcg                                                                                      >  1:2792849/1‑150 (MQ=255)
                                                        aCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcg                                                                                      <  1:3411370/150‑1 (MQ=255)
                                                        aCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcg                                                                                      <  2:3779396/150‑1 (MQ=255)
                                                         cGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcgc                                                                                     >  2:2820429/1‑150 (MQ=255)
                                                             aGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCa                                                                                 >  2:3010712/1‑150 (MQ=255)
                                                             aGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCa                                                                                 <  2:525723/150‑1 (MQ=255)
                                                              ggTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCAt                                                                                >  1:2589144/1‑150 (MQ=255)
                                                               gTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATc                                                                               >  1:326848/1‑150 (MQ=255)
                                                                 cttctGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGt                                                                             <  1:1508840/150‑1 (MQ=255)
                                                                   tctGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCt                                                                           <  1:1042954/150‑1 (MQ=255)
                                                                     tGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGc                                                                         >  1:1277222/1‑150 (MQ=255)
                                                                        aaCGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTc                                                                      >  1:1658712/1‑150 (MQ=255)
                                                                         aCGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCa                                                                     >  2:1016040/1‑150 (MQ=255)
                                                                         aCGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTATGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCa                                                                     <  2:554796/150‑1 (MQ=255)
                                                                            gCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACa                                                                  >  2:715517/1‑150 (MQ=255)
                                                                             cAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAg                                                                 >  2:3751588/1‑150 (MQ=255)
                                                                              aGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGc                                                                >  1:4367924/1‑150 (MQ=255)
                                                                               gCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGcc                                                               <  2:413360/150‑1 (MQ=255)
                                                                               gCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGcc                                                               >  1:1450145/1‑150 (MQ=255)
                                                                                cGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGTCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCca                                                              <  1:1236378/150‑1 (MQ=255)
                                                                                   gACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACaa                                                           >  1:570970/1‑150 (MQ=255)
                                                                                     cAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACACGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACg                                                         <  2:2110233/150‑1 (MQ=255)
                                                                                      agagCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGa                                                        <  1:2250538/150‑1 (MQ=255)
                                                                                       gagCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGaa                                                       <  2:1407747/150‑1 (MQ=255)
                                                                                         gcgcCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAgg                                                     <  1:3642494/150‑1 (MQ=255)
                                                                                          cgcCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                                    <  2:1499282/150‑1 (MQ=255)
                                                                                          cgcCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                                    <  2:3226858/150‑1 (MQ=255)
                                                                                           gcCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCt                                                   >  1:2242707/1‑150 (MQ=255)
                                                                                            ccTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTg                                                  <  2:4328810/150‑1 (MQ=255)
                                                                                             cTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTgg                                                 >  2:3404495/1‑150 (MQ=255)
                                                                                                 gCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGaa                                             >  1:2393300/1‑150 (MQ=255)
                                                                                                  cTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAg                                            >  2:3083933/1‑150 (MQ=255)
                                                                                                      ggCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCAc                                        >  1:1344081/1‑150 (MQ=255)
                                                                                                       gCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACt                                       <  1:163094/150‑1 (MQ=255)
                                                                                                        cAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTg                                      >  2:1747176/1‑150 (MQ=255)
                                                                                                         aGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGa                                     <  1:3559255/150‑1 (MQ=255)
                                                                                                          ggCTCAAAAACAGAGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCATGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGaa                                    >  2:1980069/1‑150 (MQ=255)
                                                                                                           gCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAAc                                   <  1:2334586/150‑1 (MQ=255)
                                                                                                            cTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACa                                  <  2:4206872/150‑1 (MQ=255)
                                                                                                            cTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACa                                  >  1:383318/1‑150 (MQ=255)
                                                                                                              cAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAgg                                <  2:1658712/150‑1 (MQ=255)
                                                                                                                aaaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGcc                              >  1:4095450/1‑150 (MQ=255)
                                                                                                                   aCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgc                           >  2:347182/1‑150 (MQ=255)
                                                                                                                    cAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgcg                          >  1:808104/1‑150 (MQ=255)
                                                                                                                    cAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgcg                          >  1:4325341/1‑150 (MQ=255)
                                                                                                                     aGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCCGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAgcgc                         >  1:3943231/1‑150 (MQ=255)
                                                                                                                      gCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgca                        <  1:1938560/150‑1 (MQ=255)
                                                                                                                        ggCTGCGTCGTCTGCCAGTGCAGCGTAATCAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACCACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcagc                      >  2:2467913/1‑150 (MQ=255)
                                                                                                                        ggCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcagc                      >  1:170201/1‑150 (MQ=255)
                                                                                                                         gCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcagca                     <  1:4575418/150‑1 (MQ=255)
                                                                                                                           tGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCagcagc                   >  2:4104255/1‑150 (MQ=255)
                                                                                                                           tGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCagcagc                   >  2:290653/1‑150 (MQ=255)
                                                                                                                             cgtcgtCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGa                 >  2:2689633/1‑150 (MQ=255)
                                                                                                                              gtcgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAg                >  1:1533062/1‑150 (MQ=255)
                                                                                                                               tcgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGNAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAgg               <  1:3471973/150‑1 (MQ=255)
                                                                                                                                cgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGt              >  2:3446398/1‑150 (MQ=255)
                                                                                                                                cgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGt              >  2:2559429/1‑150 (MQ=255)
                                                                                                                                 gtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTc             <  1:2531964/150‑1 (MQ=255)
                                                                                                                                  tcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt            <  2:2345435/150‑1 (MQ=255)
                                                                                                                                     gCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCt         >  1:4234736/1‑150 (MQ=255)
                                                                                                                                     gCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCt         >  2:2012655/1‑150 (MQ=255)
                                                                                                                                      ccAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCtt        <  1:3037635/150‑1 (MQ=255)
                                                                                                                                       cAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTc       >  1:1847473/1‑150 (MQ=255)
                                                                                                                                        aGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTcc      >  2:4411606/1‑150 (MQ=255)
                                                                                                                                        aGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTcc      <  1:676178/150‑1 (MQ=255)
                                                                                                                                        aGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTcc      <  2:4452344/150‑1 (MQ=255)
                                                                                                                                        aGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTcc      <  2:3792888/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCGCAG  >  CP009273/1424088‑1424377

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