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.3% 929.6 / 4.6 230G172A (GGC→GCC) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base G (103/123);  total (106/124)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 238 total aligned reads displayed.

GCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTT  >  CP009273/2945160‑2945451
                                                                                                                                                  |                                                                                                                                                 
   aaGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGgcgccg                                                                                                                                              >  1:2460442/1‑149 (MQ=255)
    aGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCgg                                                                                                                                             >  1:3556170/1‑149 (MQ=255)
      ccTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGaa                                                                                                                                           <  2:550501/149‑1 (MQ=255)
        ttGCCGAGAGATTCGGTTAAGTGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                         >  2:3264017/1‑149 (MQ=255)
           ccGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCAc                                                                                                                                      <  2:898434/149‑1 (MQ=255)
           ccGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCAc                                                                                                                                      <  1:990471/149‑1 (MQ=255)
             gagagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                    >  1:853089/1‑149 (MQ=255)
              agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                   >  1:2813167/1‑149 (MQ=255)
                 gaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGAt                                                                                                                                >  1:3182732/1‑149 (MQ=255)
                     cGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgc                                                                                                                            <  2:3779930/149‑1 (MQ=255)
                      ggTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgcc                                                                                                                           >  1:2861418/1‑149 (MQ=255)
                          aaGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGt                                                                                                                       >  2:2646240/1‑149 (MQ=255)
                          aaGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGt                                                                                                                       >  2:2161593/1‑149 (MQ=255)
                           aGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTc                                                                                                                      <  1:3086186/149‑1 (MQ=255)
                           aGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTc                                                                                                                      >  1:1099636/1‑149 (MQ=255)
                                gcgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAgcgc                                                                                                                 <  1:2894019/149‑1 (MQ=255)
                                 cgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGcc                                                                                                                <  1:3171293/149‑1 (MQ=255)
                                   gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCAc                                                                                                              <  1:1565387/149‑1 (MQ=255)
                                   gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCAc                                                                                                              <  1:1373298/149‑1 (MQ=255)
                                   gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCAc                                                                                                              <  1:3502541/149‑1 (MQ=255)
                                   gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCAc                                                                                                              <  2:1581507/149‑1 (MQ=255)
                                   gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCAc                                                                                                              >  2:1903433/1‑149 (MQ=255)
                                     aGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTg                                                                                                            >  1:1764179/1‑149 (MQ=255)
                                      gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgc                                                                                                           >  2:380065/1‑149 (MQ=255)
                                      gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgc                                                                                                           >  1:4000976/1‑149 (MQ=255)
                                       cTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                          <  2:2751522/149‑1 (MQ=255)
                                          tACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccgcc                                                                                                       >  2:2121861/1‑149 (MQ=255)
                                               gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                  <  1:1545646/149‑1 (MQ=255)
                                               gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                  <  1:3177512/149‑1 (MQ=255)
                                               gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                  <  1:2713631/149‑1 (MQ=255)
                                               gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                  <  2:2192184/149‑1 (MQ=255)
                                               gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                  <  1:2711223/149‑1 (MQ=255)
                                               gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                  >  1:647674/1‑149 (MQ=255)
                                                ttGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                                 <  1:2142261/149‑1 (MQ=255)
                                                ttGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                                 <  2:708571/149‑1 (MQ=255)
                                                   gTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACtttt                                                                                              <  2:681858/149‑1 (MQ=255)
                                                   gTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACtttt                                                                                              <  2:2526126/149‑1 (MQ=255)
                                                   gTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACtttt                                                                                              <  2:1945364/149‑1 (MQ=255)
                                                   gTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACtttt                                                                                              <  2:1486198/149‑1 (MQ=255)
                                                    tGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACttttt                                                                                             <  1:3780830/149‑1 (MQ=255)
                                                          cAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGt                                                                                       <  1:840633/149‑1 (MQ=255)
                                                          cAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGt                                                                                       <  2:1220220/149‑1 (MQ=255)
                                                            gACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTa                                                                                     >  2:2278709/1‑149 (MQ=255)
                                                            gACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTa                                                                                     >  2:2280091/1‑149 (MQ=255)
                                                              cGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAt                                                                                   >  1:1756043/1‑149 (MQ=255)
                                                               ggATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAtt                                                                                  >  2:1492049/1‑149 (MQ=255)
                                                                  tACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                               <  2:3161842/149‑1 (MQ=255)
                                                                  tACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                               <  1:3309099/149‑1 (MQ=255)
                                                                    cGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATc                                                                             <  1:1218313/149‑1 (MQ=255)
                                                                       cAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCt                                                                          <  2:1761145/149‑1 (MQ=255)
                                                                         gCGTTCGCCGTCGGCCATTTGAATTAACGCTGACAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGAGCCGGAAATCACCGAGATCCGCCGCGTCAGAGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                        <  2:1110555/149‑1 (MQ=255)
                                                                          cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTAc                                                                       <  1:2259641/149‑1 (MQ=255)
                                                                          cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTAc                                                                       >  2:429937/1‑149 (MQ=255)
                                                                            ttCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                     >  2:1004130/1‑149 (MQ=255)
                                                                            ttCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                     <  2:3932616/149‑1 (MQ=255)
                                                                            ttCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                     >  2:3060375/1‑149 (MQ=255)
                                                                              cGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGa                                                                   <  2:1854877/149‑1 (MQ=255)
                                                                                     ggCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGtt                                                            >  1:3216535/1‑149 (MQ=255)
                                                                                      gCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGttt                                                           <  2:2226308/149‑1 (MQ=255)
                                                                                       ccATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTg                                                          >  2:1849572/1‑149 (MQ=255)
                                                                                        cATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGt                                                         >  2:1535827/1‑149 (MQ=255)
                                                                                        cATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGt                                                         >  2:1535487/1‑149 (MQ=255)
                                                                                        cATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGt                                                         >  2:1303175/1‑149 (MQ=255)
                                                                                         aTTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                        >  2:1828722/1‑149 (MQ=255)
                                                                                            tGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAg                                                     <  2:3515159/149‑1 (MQ=255)
                                                                                              aaTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                   >  2:2943328/1‑149 (MQ=255)
                                                                                               aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGt                                                  <  1:3935708/149‑1 (MQ=255)
                                                                                               aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGt                                                  >  2:845729/1‑149 (MQ=255)
                                                                                                ttAACGCTGCCAGCAACGTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGtt                                                 <  2:261755/149‑1 (MQ=255)
                                                                                                     gCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTgcgc                                            >  1:484861/1‑149 (MQ=255)
                                                                                                      cTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGcc                                           <  2:3325511/149‑1 (MQ=255)
                                                                                                        gCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAg                                         >  1:922413/1‑149 (MQ=255)
                                                                                                          cAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                       <  2:2714124/149‑1 (MQ=255)
                                                                                                          cAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                       <  2:1670223/149‑1 (MQ=255)
                                                                                                           aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTag                                      >  2:1976866/1‑149 (MQ=255)
                                                                                                           aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTag                                      <  2:3619948/149‑1 (MQ=255)
                                                                                                                cTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCtt                                 <  2:1508534/149‑1 (MQ=255)
                                                                                                                cTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCtt                                 >  2:3400327/1‑149 (MQ=255)
                                                                                                                 ttCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtc                                <  1:3554122/149‑1 (MQ=255)
                                                                                                                   cGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGTCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcat                              <  1:3836152/149‑1 (MQ=255)
                                                                                                                   cGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGTCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcat                              <  1:86559/149‑1 (MQ=255)
                                                                                                                    gCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcatc                             >  2:1585979/1‑149 (MQ=255)
                                                                                                                    gCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcatc                             >  1:3649348/1‑149 (MQ=255)
                                                                                                                         ggtggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCt                        <  1:2551041/149‑1 (MQ=255)
                                                                                                                            ggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCaa                     <  2:2861418/149‑1 (MQ=255)
                                                                                                                             gtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAAt                    <  1:1902989/149‑1 (MQ=255)
                                                                                                                              tCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATg                   >  1:3444724/1‑149 (MQ=255)
                                                                                                                                gTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGc                 >  1:1950540/1‑149 (MQ=255)
                                                                                                                                    tACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGa             >  1:3438100/1‑149 (MQ=255)
                                                                                                                                     aCCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAt            <  2:1473323/149‑1 (MQ=255)
                                                                                                                                      ccGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAtt           <  2:3556170/149‑1 (MQ=255)
                                                                                                                                       cGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTc          <  2:661789/149‑1 (MQ=255)
                                                                                                                                       cGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTc          <  1:2010075/149‑1 (MQ=255)
                                                                                                                                       cGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGAGCCACTGCCGCCGCAACTTTTTGCCTGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTc          <  2:2724547/149‑1 (MQ=255)
                                                                                                                                        ggTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCa         <  2:488186/149‑1 (MQ=255)
                                                                                                                                         gTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCAc        >  1:2970742/1‑149 (MQ=255)
                                                                                                                                             cAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAtt           >  1:1145913/1‑142 (MQ=255)
                                                                                                                                             cAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCg    >  2:3654191/1‑149 (MQ=255)
                                                                                                                                               ggggCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGtt  >  1:2721824/1‑149 (MQ=255)
                                                                                                                                               ggggCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGtt  >  1:2320615/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGTT  >  CP009273/2945160‑2945451

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: