Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,512,693 T→A 100% K407* (AAG→TAG)  yjhG ← putative dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,512,6930TA95.7% 1009.3 / 25.0 300K407* (AAG→TAG) yjhGputative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (3/10);  new base A (168/119);  total (171/129)
Fisher's exact test for biased strand distribution p-value = 1.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 300 total aligned reads displayed.

TGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAT  >  CP009273/4512546‑4512839
                                                                                                                                                   |                                                                                                                                                  
 ggCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGaa                                                                                                                                                 <  2:2215840/150‑1 (MQ=255)
   cGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAAcc                                                                                                                                               <  1:180971/150‑1 (MQ=255)
      tGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCt                                                                                                                                            >  1:500475/1‑150 (MQ=255)
        gATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGTGTTGCTAGAACCGCTGa                                                                                                                                          >  2:1561202/1‑150 (MQ=255)
          tttGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACg                                                                                                                                        <  1:2644016/150‑1 (MQ=255)
          tttGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACg                                                                                                                                        >  2:1124087/1‑150 (MQ=255)
           ttGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGc                                                                                                                                       >  2:3652713/1‑150 (MQ=255)
               tCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGCTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGtt                                                                                                                                   <  1:3993308/150‑1 (MQ=255)
                cACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTc                                                                                                                                  >  1:52425/1‑150 (MQ=255)
                    gAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAg                                                                                                                              <  1:2811122/150‑1 (MQ=255)
                     aaCCTTCTGGCGCAATATTGCCTACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGt                                                                                                                             >  1:1307799/1‑150 (MQ=255)
                     aaCCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGt                                                                                                                             >  2:3213156/1‑150 (MQ=255)
                          tCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTcc                                                                                                                        <  2:4411533/150‑1 (MQ=255)
                           cTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTccc                                                                                                                       <  1:33923/150‑1 (MQ=255)
                           cTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTccc                                                                                                                       >  1:3091141/1‑150 (MQ=255)
                            tGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCcca                                                                                                                      >  1:112167/1‑150 (MQ=255)
                              gcgcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCccacc                                                                                                                    >  2:4598557/1‑150 (MQ=255)
                              gcgcAATAGTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCccacc                                                                                                                    <  2:3932240/150‑1 (MQ=255)
                               cgcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCccacca                                                                                                                   <  1:3343976/150‑1 (MQ=255)
                                gcAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAg                                                                                                                  <  1:4075734/150‑1 (MQ=255)
                                 cAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGt                                                                                                                 >  1:1282237/1‑150 (MQ=255)
                                  aaTATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTc                                                                                                                >  1:1471898/1‑150 (MQ=255)
                                  aaTATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTc                                                                                                                >  1:3128756/1‑150 (MQ=255)
                                  aaTATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTACGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTc                                                                                                                >  2:4080399/1‑150 (MQ=255)
                                   atatTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCg                                                                                                               <  2:4169553/150‑1 (MQ=255)
                                    tatTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGa                                                                                                              >  1:1442990/1‑150 (MQ=255)
                                      ttGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAgg                                                                                                            >  2:1781437/1‑150 (MQ=255)
                                         cccACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGttt                                                                                                         >  2:3072194/1‑150 (MQ=255)
                                           cACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTtttc                                                                                                       >  2:2036594/1‑150 (MQ=255)
                                                ggAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCa                                                                                                  >  1:409492/1‑150 (MQ=255)
                                                 gAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAg                                                                                                 >  2:1874185/1‑150 (MQ=255)
                                                 gAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAg                                                                                                 <  1:3615149/150‑1 (MQ=255)
                                                  aaGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGc                                                                                                >  2:3395699/1‑150 (MQ=255)
                                                  aaGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGc                                                                                                >  1:354797/1‑150 (MQ=255)
                                                        aTAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGtgctg                                                                                          >  2:2908963/1‑150 (MQ=255)
                                                           gTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCg                                                                                       >  1:2160986/1‑150 (MQ=255)
                                                                ggTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTGAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAAc                                                                                  <  2:583173/150‑1 (MQ=255)
                                                                    aaTCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTc                                                                              >  2:1006526/1‑150 (MQ=255)
                                                                     aTCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCa                                                                             <  1:4481865/150‑1 (MQ=255)
                                                                       ccGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATa                                                                           >  2:1062804/1‑150 (MQ=255)
                                                                          cgcgcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACg                                                                        <  1:3795575/150‑1 (MQ=255)
                                                                            cgcTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTc                                                                      <  2:1382810/150‑1 (MQ=255)
                                                                              cTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCtt                                                                    <  2:2328942/150‑1 (MQ=255)
                                                                              cTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCtt                                                                    <  2:2814522/150‑1 (MQ=255)
                                                                              cTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCtt                                                                    >  1:3031617/1‑150 (MQ=255)
                                                                               ttttGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTc                                                                   >  2:3805908/1‑150 (MQ=255)
                                                                                tttGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCa                                                                  <  2:475214/150‑1 (MQ=255)
                                                                                tttGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCa                                                                  >  2:4459121/1‑150 (MQ=255)
                                                                                tttGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCa                                                                  <  1:2513674/150‑1 (MQ=255)
                                                                                 ttgcttgcTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCAt                                                                 >  1:3614191/1‑150 (MQ=255)
                                                                                   gcttgctGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGc                                                               >  1:659925/1‑150 (MQ=255)
                                                                                     ttgctGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaa                                                             >  1:2484484/1‑150 (MQ=255)
                                                                                      tgctgcGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGcaac                                                            <  1:153159/150‑1 (MQ=255)
                                                                                         tgcGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAAt                                                         <  2:3041140/150‑1 (MQ=255)
                                                                                            ggAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCg                                                      >  1:4515490/1‑150 (MQ=255)
                                                                                            ggAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCg                                                      >  2:235316/1‑150 (MQ=255)
                                                                                             gagaCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGa                                                     >  1:460501/1‑150 (MQ=255)
                                                                                             gagaCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGa                                                     >  1:3582273/1‑150 (MQ=255)
                                                                                              agaCATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAg                                                    >  2:3281104/1‑150 (MQ=255)
                                                                                                  aTGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTg                                                >  1:3564723/1‑150 (MQ=255)
                                                                                                  aTGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTg                                                >  2:2029677/1‑150 (MQ=255)
                                                                                                    gATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTgcg                                              >  1:2711005/1‑150 (MQ=255)
                                                                                                     aTCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTgcgc                                             <  2:650253/150‑1 (MQ=255)
                                                                                                      tCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCa                                            >  2:845941/1‑150 (MQ=255)
                                                                                                      tCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCa                                            >  1:3972283/1‑150 (MQ=255)
                                                                                                         cTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGAt                                         >  1:3420601/1‑150 (MQ=255)
                                                                                                             gTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCaa                                     >  1:1839217/1‑150 (MQ=255)
                                                                                                              tCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                    >  1:1714865/1‑150 (MQ=255)
                                                                                                              tCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                    >  2:4502934/1‑150 (MQ=255)
                                                                                                              tCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAAc                                    >  2:3317849/1‑150 (MQ=255)
                                                                                                                aGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACAt                                  <  2:1925762/150‑1 (MQ=255)
                                                                                                                  cgattaTTTGTTCCTGATCGAGCAGGAGTAGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTGTTCTTTCAGCGTGCTGCCTGTAACAGTCATACAGTCTTCATGCAACAATCCGAAGCTGCGCAGATGCAACATGa                                >  2:2075175/6‑150 (MQ=255)
                                                                                                                   gTTGATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGAc                               >  2:1160436/1‑150 (MQ=255)
                                                                                                                      gATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTc                            >  2:3829928/1‑150 (MQ=255)
                                                                                                                      gATTTGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTc                            >  1:3461520/1‑150 (MQ=255)
                                                                                                                        tttGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCg                          >  1:4327835/1‑150 (MQ=255)
                                                                                                                        tttGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCg                          <  2:4273065/150‑1 (MQ=255)
                                                                                                                         ttGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCgg                         >  1:2920114/1‑150 (MQ=255)
                                                                                                                          tGTTCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGc                        >  1:4191404/1‑150 (MQ=255)
                                                                                                                            ttCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGcac                      >  1:1162139/1‑150 (MQ=255)
                                                                                                                            ttCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGcac                      <  2:3567529/150‑1 (MQ=255)
                                                                                                                            ttCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGcac                      >  2:3443258/1‑150 (MQ=255)
                                                                                                                             tCCTGATCGAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGcaca                     >  1:1898446/1‑150 (MQ=255)
                                                                                                                             tCCTGATCGAGCAGGAGTTGCTAGAACCGCNGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGcaca                     >  1:4280398/1‑150 (MQ=255)
                                                                                                                                     gAGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGc             >  1:2052824/1‑150 (MQ=255)
                                                                                                                                      aGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGcc            >  1:2354868/1‑150 (MQ=255)
                                                                                                                                      aGCAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGcc            <  1:782439/150‑1 (MQ=255)
                                                                                                                                        cAGGAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCAt          <  2:3288225/150‑1 (MQ=255)
                                                                                                                                          ggAGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATaa        >  2:1153661/1‑150 (MQ=255)
                                                                                                                                            aGTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAg      <  2:430579/150‑1 (MQ=255)
                                                                                                                                             gTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGTGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAgg     <  1:469854/150‑1 (MQ=255)
                                                                                                                                             gTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAgg     <  1:2621383/150‑1 (MQ=255)
                                                                                                                                             gTTGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAgg     <  2:3976367/150‑1 (MQ=255)
                                                                                                                                              ttGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGc    <  1:2265064/150‑1 (MQ=255)
                                                                                                                                               tGCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCa   <  1:3148594/150‑1 (MQ=255)
                                                                                                                                                gCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAt  <  2:1443312/150‑1 (MQ=255)
                                                                                                                                                gCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAt  >  1:2008042/1‑150 (MQ=255)
                                                                                                                                                gCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAt  <  2:2738916/150‑1 (MQ=255)
                                                                                                                                                gCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAt  >  1:1394421/1‑150 (MQ=255)
                                                                                                                                                gCTAGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAt  <  2:932322/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGCGGTGGATTTGATCACCGAACCTTCTGGCGCAATATTGCCCACCGGGAAGGTGATAGTTGAGGTTAATCCGCGCGCTTTTGCTTGCTGCGGAGACATGATCACTTCGTCAGCGTTGATTTGTTCCTGATCGAGCAGGAGTTGCTTGAACCGCTGACGCCGTTCGGAGTGCTCCCACCAGTCGAGGTTTTCTTTCAGCGTGCTGCCGGTAACCGTCATAACGTCTTCATGCAACAATCCGAGGCTGCGCAGATGCAACATGACTTCCGGCACACCACCTGCCATAAAGGCAT  >  CP009273/4512546‑4512839

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