Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,257,232 G→T 100% R34S (CGC→AGC)  prs ← phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,257,2320GT97.3% 686.9 / 8.7 221R34S (CGC→AGC) prsphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (120/95);  total (123/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

TGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGA  >  CP009273/1257086‑1257377
                                                                                                                                                  |                                                                                                                                                 
tGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAc                                                                                                                                                 <  1:3455946/149‑1 (MQ=255)
tGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAc                                                                                                                                                 <  1:2515310/149‑1 (MQ=255)
 gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTAcc                                                                                                                                                >  2:2256346/1‑149 (MQ=255)
    gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                             >  2:1124310/1‑149 (MQ=255)
    gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                             >  2:709579/1‑149 (MQ=255)
    gAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGctacctac                                                                                                                                             >  2:1493942/1‑149 (MQ=255)
     aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACa                                                                                                                                            >  2:1755802/1‑149 (MQ=255)
           gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCAGTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAgcggcg                                                                                                                                      >  1:2370462/1‑149 (MQ=255)
             cgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTc                                                                                                                                    <  1:3046686/149‑1 (MQ=255)
               cAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGc                                                                                                                                  <  1:2747610/149‑1 (MQ=255)
                     aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTg                                                                                                                            >  1:3170019/1‑149 (MQ=255)
                     aTCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTg                                                                                                                            >  2:1350378/1‑149 (MQ=255)
                      tCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGa                                                                                                                           >  2:3163229/1‑149 (MQ=255)
                       cAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGaa                                                                                                                          >  1:2344272/1‑149 (MQ=255)
                         aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAgt                                                                                                                        >  2:3477206/1‑149 (MQ=255)
                          ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAgtg                                                                                                                       <  1:1146057/149‑1 (MQ=255)
                          ccATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAgtg                                                                                                                       >  2:1575723/1‑149 (MQ=255)
                            aTAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTa                                                                                                                     <  1:1270841/149‑1 (MQ=255)
                             tAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTAc                                                                                                                    >  2:1540445/1‑149 (MQ=255)
                              aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACa                                                                                                                   <  2:481172/149‑1 (MQ=255)
                                cGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAgg                                                                                                                 >  2:3924953/1‑149 (MQ=255)
                                    tAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGt                                                                                                             <  1:1337616/149‑1 (MQ=255)
                                      aTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTg                                                                                                           >  2:2392851/1‑149 (MQ=255)
                                        tCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGc                                                                                                         >  1:2039779/1‑149 (MQ=255)
                                          cATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCaa                                                                                                       <  1:1875678/149‑1 (MQ=255)
                                           aTCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAt                                                                                                      >  1:3075741/1‑149 (MQ=255)
                                            tCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATa                                                                                                     <  1:2644995/149‑1 (MQ=255)
                                             cAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTACGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATAc                                                                                                    <  2:1125901/149‑1 (MQ=255)
                                              aGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACg                                                                                                   <  1:4025534/149‑1 (MQ=255)
                                                 ttGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                                >  1:828779/1‑149 (MQ=255)
                                                 ttGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                                <  1:637731/149‑1 (MQ=255)
                                                 ttGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtg                                                                                                <  1:3215190/149‑1 (MQ=255)
                                                  tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTtgt                                                                                               <  1:1022183/149‑1 (MQ=255)
                                                    tCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGc                                                                                             <  1:3095047/149‑1 (MQ=255)
                                                           tagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTcc                                                                                      <  2:2590773/149‑1 (MQ=255)
                                                             ggggCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCgg                                                                                    <  1:3330221/149‑1 (MQ=255)
                                                                gCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGt                                                                                 >  2:1200149/1‑149 (MQ=255)
                                                                 cacaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTg                                                                                <  1:1607322/149‑1 (MQ=255)
                                                                  acaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAGGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTgg                                                                               >  2:3605301/1‑149 (MQ=255)
                                                                   caAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGc                                                                              >  2:38608/1‑149 (MQ=255)
                                                                     aGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGt                                                                            >  1:2627200/1‑149 (MQ=255)
                                                                      gTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGtt                                                                           >  2:3329212/1‑149 (MQ=255)
                                                                       tGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTa                                                                          <  2:3594703/149‑1 (MQ=255)
                                                                       tGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTa                                                                          <  1:3477206/149‑1 (MQ=255)
                                                                          aCTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCa                                                                       <  2:2362911/149‑1 (MQ=255)
                                                                            tGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGc                                                                     >  2:3999723/1‑149 (MQ=255)
                                                                             ggATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCa                                                                    >  2:2007968/1‑149 (MQ=255)
                                                                                tgatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaa                                                                 >  1:2915914/1‑149 (MQ=255)
                                                                                tgatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaa                                                                 <  1:454752/149‑1 (MQ=255)
                                                                                 gatgaAAATATCACCACCGCGTACCTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaaa                                                                >  2:3041003/1‑149 (MQ=255)
                                                                                 gatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCaaaaa                                                                <  2:3748706/149‑1 (MQ=255)
                                                                                  atgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAg                                                               >  1:3698450/1‑149 (MQ=255)
                                                                                     aaaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCtt                                                            <  1:697116/149‑1 (MQ=255)
                                                                                      aaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                           <  1:1028917/149‑1 (MQ=255)
                                                                                      aaaTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                           <  1:3498311/149‑1 (MQ=255)
                                                                                         tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCata                                                        >  1:357727/1‑149 (MQ=255)
                                                                                          atCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCatat                                                       <  1:3550557/149‑1 (MQ=255)
                                                                                           tCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATc                                                      <  2:1993596/149‑1 (MQ=255)
                                                                                            caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCa                                                     <  1:2712356/149‑1 (MQ=255)
                                                                                              ccaccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAgg                                                   >  2:2885353/1‑149 (MQ=255)
                                                                                              ccaccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAgg                                                   >  2:824739/1‑149 (MQ=255)
                                                                                               caccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGc                                                  >  1:1116433/1‑149 (MQ=255)
                                                                                                accGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCa                                                 >  1:221162/1‑149 (MQ=255)
                                                                                                accGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCa                                                 <  2:140237/149‑1 (MQ=255)
                                                                                                 ccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCAc                                                >  2:432421/1‑149 (MQ=255)
                                                                                                   gcgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACga                                              <  2:20503/149‑1 (MQ=255)
                                                                                                      tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaa                                           >  1:4009277/1‑149 (MQ=255)
                                                                                                       aCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaag                                          >  1:1233309/1‑149 (MQ=255)
                                                                                                        cATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaaga                                         >  2:1819418/1‑149 (MQ=255)
                                                                                                         aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGAAAAAAGCTTCATATCAGGCACGAgaagaa                                        <  2:3520051/149‑1 (MQ=255)
                                                                                                              cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgag                                             >  1:1235948/1‑139 (MQ=255)
                                                                                                              cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACgag                                             <  2:1235948/139‑1 (MQ=255)
                                                                                                              cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCa                                   >  1:1446717/1‑149 (MQ=255)
                                                                                                              cATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCa                                   <  2:1268481/149‑1 (MQ=255)
                                                                                                               aTTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAg                                  <  1:1124353/149‑1 (MQ=255)
                                                                                                               aTTAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAg                                  >  2:3328042/1‑149 (MQ=255)
                                                                                                                ttAATTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAgg                                 <  2:107053/149‑1 (MQ=255)
                                                                                                                  aaTTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCa                               >  2:1525921/1‑149 (MQ=255)
                                                                                                                   aTTTGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCAt                              >  2:4045141/1‑149 (MQ=255)
                                                                                                                     ttGTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGc                            <  1:3420398/149‑1 (MQ=255)
                                                                                                                       gTACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGt                          >  1:1739842/1‑149 (MQ=255)
                                                                                                                        tACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTc                         >  2:246609/1‑149 (MQ=255)
                                                                                                                        tACGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCACAGGCATGCGTc                         >  1:721009/1‑149 (MQ=255)
                                                                                                                         aCGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTcc                        <  1:3643823/149‑1 (MQ=255)
                                                                                                                          cGCTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCa                       <  1:1575723/149‑1 (MQ=255)
                                                                                                                            cTGACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAtt                     >  1:3399572/1‑149 (MQ=255)
                                                                                                                              gACTTCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtgg                   >  1:562379/1‑149 (MQ=255)
                                                                                                                                catCGCCATCGCTAAAGCTACCTACAGCGGCGTCGCCCACTGAAGTGTACACGCGGATGGCAATACGTTGTGCTAGTTCCGGGGAGGCGTTACCAGCAAAAAGCTTCATATCAGGCAAGAGAAGAACCTAAGGCATGCGTCCAGaggtg                 >  1:3433025/3‑149 (MQ=255)
                                                                                                                                    gCCATCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGaaa             <  2:2933379/149‑1 (MQ=255)
                                                                                                                                       aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa          <  2:1282787/149‑1 (MQ=255)
                                                                                                                                       aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa          <  2:646777/149‑1 (MQ=255)
                                                                                                                                       aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGaa          >  1:3555958/1‑149 (MQ=255)
                                                                                                                                       aTCGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCAGTGGTGGAAAGaa          >  1:981246/1‑149 (MQ=255)
                                                                                                                                         cGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATc        <  1:3232160/149‑1 (MQ=255)
                                                                                                                                         cGCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATc        >  2:2737941/1‑149 (MQ=255)
                                                                                                                                          gCTAAAGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCt       <  2:2071687/149‑1 (MQ=255)
                                                                                                                                             aaaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGcc    >  2:383047/1‑149 (MQ=255)
                                                                                                                                             aaaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGcc    >  1:1794884/1‑149 (MQ=255)
                                                                                                                                              aaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCg   >  2:3492170/1‑149 (MQ=255)
                                                                                                                                              aaGCTACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCg   >  2:1860762/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGA  >  CP009273/1257086‑1257377

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: