Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,945,306 C→G 100% G172A (GGC→GCC)  recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,945,3060CG98.6% 1090.2 / 4.8 283G172A (GGC→GCC) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base G (132/146);  total (135/148)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 292 total aligned reads displayed.

CGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTG  >  CP009273/2945159‑2945452
                                                                                                                                                   |                                                                                                                                                  
cGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGggcg                                                                                                                                                  >  2:2869491/1‑150 (MQ=255)
     aGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGa                                                                                                                                             <  2:4232378/150‑1 (MQ=255)
     aGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGa                                                                                                                                             >  1:2134666/1‑150 (MQ=255)
        cTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                          <  1:1968062/150‑1 (MQ=255)
        cTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                          >  1:3245648/1‑150 (MQ=255)
        cTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                          <  1:1874080/150‑1 (MQ=255)
         ttGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATc                                                                                                                                         <  1:2504827/150‑1 (MQ=255)
          tGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCa                                                                                                                                        <  2:1765306/150‑1 (MQ=255)
           gCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCAc                                                                                                                                       <  1:3367030/150‑1 (MQ=255)
             cGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                     >  1:1449272/1‑150 (MQ=255)
             cGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                     >  1:5094185/1‑150 (MQ=255)
              gagagaGTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                    <  1:4676030/150‑1 (MQ=255)
               agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCgag                                                                                                                                   >  1:4515034/1‑150 (MQ=255)
                   aTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATcc                                                                                                                               <  1:4378460/150‑1 (MQ=255)
                    gtCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccg                                                                                                                              <  1:4151570/149‑1 (MQ=255)
                      tggTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgcc                                                                                                                            >  1:4442766/2‑150 (MQ=255)
                      cGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgcc                                                                                                                            >  2:2304272/1‑150 (MQ=255)
                       ggTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgccg                                                                                                                           <  1:3639939/150‑1 (MQ=255)
                            aGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCa                                                                                                                      >  2:1280738/1‑150 (MQ=255)
                                cgcgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAgcgc                                                                                                                  <  1:3349648/150‑1 (MQ=255)
                                 gcgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGcc                                                                                                                 >  2:3057847/1‑150 (MQ=255)
                                  cgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCa                                                                                                                >  1:2191423/1‑150 (MQ=255)
                                  cgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCa                                                                                                                >  1:3042477/1‑150 (MQ=255)
                                    gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACt                                                                                                              >  2:1043252/1‑150 (MQ=255)
                                     cAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTg                                                                                                             >  2:3793855/1‑150 (MQ=255)
                                      aGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgc                                                                                                            >  1:1337905/1‑150 (MQ=255)
                                       gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                           <  1:580685/150‑1 (MQ=255)
                                       gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                           <  2:3874948/150‑1 (MQ=255)
                                       gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                           <  1:1376352/150‑1 (MQ=255)
                                       gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                           <  1:1594963/150‑1 (MQ=255)
                                       gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                           >  2:2751175/1‑150 (MQ=255)
                                             cccGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCa                                                                                                     >  1:847639/1‑150 (MQ=255)
                                              ccGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCaa                                                                                                    <  1:4648451/150‑1 (MQ=255)
                                               cGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                   <  2:1155991/150‑1 (MQ=255)
                                                gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                                  >  1:993610/1‑150 (MQ=255)
                                                   ggTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACtttt                                                                                               <  1:4828995/150‑1 (MQ=255)
                                                   ggTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACtttt                                                                                               >  1:4625117/1‑150 (MQ=255)
                                                    gTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACttttt                                                                                              <  1:654080/150‑1 (MQ=255)
                                                       cAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGcc                                                                                           <  2:319585/150‑1 (MQ=255)
                                                       cAGCCAGACGGAGACGGCAGCGTTCGCCGTCGGCCATGTGAATTAACGCTGCCAGCAACTTCGCGACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGcc                                                                                           <  2:4392127/150‑1 (MQ=255)
                                                         gCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAg                                                                                         <  1:1262556/150‑1 (MQ=255)
                                                            aGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTa                                                                                      >  1:4863510/1‑150 (MQ=255)
                                                             gACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTaa                                                                                     <  1:3021473/150‑1 (MQ=255)
                                                              aCGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAt                                                                                    <  1:3666128/150‑1 (MQ=255)
                                                              aCGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAt                                                                                    <  1:2829119/150‑1 (MQ=255)
                                                               cGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAtt                                                                                   <  2:2149207/150‑1 (MQ=255)
                                                                 gATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTc                                                                                 <  2:1371767/150‑1 (MQ=255)
                                                                  aTACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                                >  1:3828299/1‑150 (MQ=255)
                                                                  aTACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                                >  1:3257402/1‑150 (MQ=255)
                                                                   tACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCAt                                                                               <  1:3914798/150‑1 (MQ=255)
                                                                    aCGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATc                                                                              >  1:1726426/1‑150 (MQ=255)
                                                                         aGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                         >  1:3127412/1‑150 (MQ=255)
                                                                         aGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                         <  1:1823256/150‑1 (MQ=255)
                                                                         aGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                         >  2:2540657/1‑150 (MQ=255)
                                                                          gCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTAc                                                                        <  2:836216/150‑1 (MQ=255)
                                                                           cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACt                                                                       >  2:3995848/1‑150 (MQ=255)
                                                                           cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACt                                                                       >  1:1906756/1‑150 (MQ=255)
                                                                           cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACt                                                                       >  2:1906562/1‑150 (MQ=255)
                                                                            gTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                      <  1:1038484/150‑1 (MQ=255)
                                                                            gTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                      >  2:295090/1‑150 (MQ=255)
                                                                            gTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                      <  1:2555652/150‑1 (MQ=255)
                                                                             ttCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTgg                                                                     <  2:2268140/150‑1 (MQ=255)
                                                                              tCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGa                                                                    <  2:17997/150‑1 (MQ=255)
                                                                                 ccGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGaaaa                                                                 <  1:3494125/150‑1 (MQ=255)
                                                                                      ggCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGttt                                                            >  1:1217491/1‑150 (MQ=255)
                                                                                      ggCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGttt                                                            <  2:1625916/150‑1 (MQ=255)
                                                                                       gCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTg                                                           >  2:4721590/1‑150 (MQ=255)
                                                                                       gCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTg                                                           <  2:3347178/150‑1 (MQ=255)
                                                                                         cATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                         >  2:4080020/1‑150 (MQ=255)
                                                                                          aTTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTcc                                                        <  1:67233/150‑1 (MQ=255)
                                                                                          aTTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTcc                                                        <  1:1153089/150‑1 (MQ=255)
                                                                                             tGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAgg                                                     >  1:4607276/1‑150 (MQ=255)
                                                                                             tGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAgg                                                     >  1:2934589/1‑150 (MQ=255)
                                                                                                  tAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTg                                                <  1:4226201/150‑1 (MQ=255)
                                                                                                  tAACGCTGCCAGCAACGTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTg                                                <  2:2144480/150‑1 (MQ=255)
                                                                                                   aaCGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTgc                                               <  2:4357225/150‑1 (MQ=255)
                                                                                                      gCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGcc                                            <  1:705315/150‑1 (MQ=255)
                                                                                                        tGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAg                                          >  1:5139338/1‑150 (MQ=255)
                                                                                                        tGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAg                                          >  2:4853477/1‑150 (MQ=255)
                                                                                                          ccAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                        >  1:4283872/1‑150 (MQ=255)
                                                                                                          ccAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                        <  1:2773294/150‑1 (MQ=255)
                                                                                                            aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTaga                                      >  1:2411985/1‑150 (MQ=255)
                                                                                                            aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTaga                                      >  1:3622838/1‑150 (MQ=255)
                                                                                                             gCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTagag                                     >  2:4108932/1‑150 (MQ=255)
                                                                                                               aaCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCt                                   >  2:3313660/1‑150 (MQ=255)
                                                                                                                  ttCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtca                                <  1:1280322/150‑1 (MQ=255)
                                                                                                                    cGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcatc                              >  1:2567642/1‑150 (MQ=255)
                                                                                                                        aCGGTGGTCGTTTTACCGGTGCCAGGGGAGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAAcc                          <  2:5069956/150‑1 (MQ=255)
                                                                                                                         cGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCt                         <  2:196314/150‑1 (MQ=255)
                                                                                                                          ggtggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTc                        <  2:586945/150‑1 (MQ=255)
                                                                                                                             ggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAAt                     <  2:1228354/150‑1 (MQ=255)
                                                                                                                                 gTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCa                 >  1:4334218/1‑150 (MQ=255)
                                                                                                                                     tACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAt             >  2:1435515/1‑150 (MQ=255)
                                                                                                                                      aCCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAtt            >  2:121959/1‑150 (MQ=255)
                                                                                                                                      aCCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAtt            >  1:233980/1‑150 (MQ=255)
                                                                                                                                       ccGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTc           >  1:3553574/1‑150 (MQ=255)
                                                                                                                                       ccGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTc           >  2:2076800/1‑150 (MQ=255)
                                                                                                                                        cGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCa          <  2:3672888/150‑1 (MQ=255)
                                                                                                                                         ggTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCAc         >  2:3107687/1‑150 (MQ=255)
                                                                                                                                               aGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGtt   <  1:4974400/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTTG  >  CP009273/2945159‑2945452

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