Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,956,165 C→A 18.3% E475D (GAG→GAT gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,956,1650CA18.3% 893.8 / 128.4 293E475D (GAG→GATgppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (129/110);  new base A (21/33);  total (150/143)
Fisher's exact test for biased strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 303 total aligned reads displayed.

GCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTT  >  CP009273/3956017‑3956313
                                                                                                                                                    |                                                                                                                                                    
 cGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAAtctc                                                                                                                                                    >  2:3418550/1‑150 (MQ=255)
 cGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCAACGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAAtctc                                                                                                                                                    >  1:1552629/1‑150 (MQ=255)
  gAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATctct                                                                                                                                                   <  1:3446531/150‑1 (MQ=255)
  gAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATctct                                                                                                                                                   <  2:3578168/150‑1 (MQ=255)
   aaTCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCtt                                                                                                                                                  <  1:631027/150‑1 (MQ=255)
    aTCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCttt                                                                                                                                                 <  1:1883027/150‑1 (MQ=255)
       tGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCCCTGGCTTTCCTGAGCAATAATCTCTTTAcc                                                                                                                                              <  2:2272002/150‑1 (MQ=255)
         aaaGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCa                                                                                                                                            <  1:844245/150‑1 (MQ=255)
         aaaGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCa                                                                                                                                            >  2:4528877/1‑150 (MQ=255)
         aaaGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCa                                                                                                                                            >  1:52339/1‑150 (MQ=255)
            gCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCg                                                                                                                                         <  2:2450619/150‑1 (MQ=255)
               atgttgGCGGCAAAACGTGTGAAATATTCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATTGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGAt                                                                                                                                      <  1:2344233/150‑1 (MQ=255)
                 gatgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGt                                                                                                                                    <  2:2849364/150‑1 (MQ=255)
                  atgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGtt                                                                                                                                   <  2:1199776/150‑1 (MQ=255)
                  atgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGtt                                                                                                                                   <  1:3437427/150‑1 (MQ=255)
                   tgGCGGCAACACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTg                                                                                                                                  <  1:1664635/150‑1 (MQ=255)
                        gcAAACCGCGGGAAATTAGCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTCATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGAAATAATCTCTTTACCCAGCGGATGTTGGGTTa                                                                                                                             <  1:199834/150‑1 (MQ=255)
                         cAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAg                                                                                                                            >  2:818628/1‑150 (MQ=255)
                           aaaCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGcc                                                                                                                          <  1:2289053/150‑1 (MQ=255)
                           aaaCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGcc                                                                                                                          >  1:3497940/1‑150 (MQ=255)
                            aaCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCa                                                                                                                         >  1:3768680/1‑150 (MQ=255)
                            aaCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCa                                                                                                                         <  2:3298218/150‑1 (MQ=255)
                               gAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAAcc                                                                                                                      >  1:2595568/1‑150 (MQ=255)
                                aGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCt                                                                                                                     <  1:664371/150‑1 (MQ=255)
                                     aataatCACGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACGTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCgg                                                                                                                <  2:1722170/150‑1 (MQ=255)
                                     aataatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCgg                                                                                                                <  2:1678532/150‑1 (MQ=255)
                                      ataatCAAGGCCGGAAAAGTTGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGAGGATGTTGGGTTAGCCAACCTTGCGGa                                                                                                               >  2:1550513/1‑150 (MQ=255)
                                      ataatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGa                                                                                                               <  2:4043433/150‑1 (MQ=255)
                                       taatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGaa                                                                                                              <  1:4586665/150‑1 (MQ=255)
                                          tCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCg                                                                                                           >  2:1192129/1‑150 (MQ=255)
                                            aaGGCCGGAAAAGATGCGTTAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTc                                                                                                         >  2:3469564/1‑150 (MQ=255)
                                              ggCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCaa                                                                                                       >  2:2797846/1‑150 (MQ=255)
                                                 cGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGt                                                                                                    >  2:1365339/1‑150 (MQ=255)
                                                   gAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACGTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCa                                                                                                  <  2:1933820/150‑1 (MQ=255)
                                                     aaaGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAAc                                                                                                >  1:169872/1‑150 (MQ=255)
                                                      aaGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACa                                                                                               <  1:2682958/150‑1 (MQ=255)
                                                              tCAGCATCGCATCCGGCACTTACTCTTTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATgg                                                                                       <  1:741342/150‑1 (MQ=255)
                                                              tCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATgg                                                                                       <  1:1377559/150‑1 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                     >  2:2908461/1‑150 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                     >  1:54162/1‑150 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                     >  2:1296665/1‑150 (MQ=255)
                                                                 gcatcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTa                                                                                    <  1:2127320/150‑1 (MQ=255)
                                                                 gcatcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTa                                                                                    >  1:1779456/1‑150 (MQ=255)
                                                                      gcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTg                                                                               <  1:2783009/150‑1 (MQ=255)
                                                                      gcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACAGAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTg                                                                               <  1:2656013/150‑1 (MQ=255)
                                                                       catcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGt                                                                              <  1:2445887/150‑1 (MQ=255)
                                                                        atcCGGCCCTTACTCATTAATGCACGTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTATACCCAGCGGAGGTTGGGTTAGCCAACCTTGAGGAAGCGTCAACGTCAACAGTGCATGGTTAGCCTGTa                                                                             >  2:770752/1‑150 (MQ=255)
                                                                        atcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTa                                                                             >  2:3148337/1‑150 (MQ=255)
                                                                        atcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTa                                                                             <  2:48499/150‑1 (MQ=255)
                                                                         tcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACGGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTaa                                                                            <  2:1404504/150‑1 (MQ=255)
                                                                           cGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATg                                                                          <  2:3074234/150‑1 (MQ=255)
                                                                            ggCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGt                                                                         >  2:3841264/1‑150 (MQ=255)
                                                                            ggCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGt                                                                         >  1:981779/1‑150 (MQ=255)
                                                                              cACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCa                                                                       >  1:2070318/1‑150 (MQ=255)
                                                                               aCTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAt                                                                      <  2:952794/150‑1 (MQ=255)
                                                                                cTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtc                                                                     <  1:1590953/150‑1 (MQ=255)
                                                                                 ttACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtct                                                                    <  2:1569950/150‑1 (MQ=255)
                                                                                  tACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtctc                                                                   >  1:1091652/1‑150 (MQ=255)
                                                                                  tACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTAAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtctc                                                                   >  2:4287536/1‑150 (MQ=255)
                                                                                    cTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTg                                                                 >  2:2301089/1‑150 (MQ=255)
                                                                                    cTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCACTGCCACAGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTg                                                                 >  1:4243232/1‑150 (MQ=255)
                                                                                     tCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTgg                                                                >  1:1634598/1‑150 (MQ=255)
                                                                                       aTTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCa                                                              >  1:4412485/1‑150 (MQ=255)
                                                                                       aTTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCa                                                              >  2:210082/1‑150 (MQ=255)
                                                                                         tAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACg                                                            >  1:2405848/1‑150 (MQ=255)
                                                                                         tAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACg                                                            <  1:2157805/150‑1 (MQ=255)
                                                                                          aaTGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACga                                                           >  2:2590548/1‑150 (MQ=255)
                                                                                           aTGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACgag                                                          >  1:3311131/1‑150 (MQ=255)
                                                                                            tGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGATAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACgaga                                                         >  2:4162720/1‑150 (MQ=255)
                                                                                            tGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACgaga                                                         <  2:3398121/150‑1 (MQ=255)
                                                                                                cgtCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGt                                                     <  2:555991/148‑1 (MQ=255)
                                                                                                   ccAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCAc                                                  <  2:701244/150‑1 (MQ=255)
                                                                                                   ccAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCAc                                                  <  1:1261376/150‑1 (MQ=255)
                                                                                                    cAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACg                                                 >  1:1471806/1‑150 (MQ=255)
                                                                                                      gCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGc                                               <  1:4101339/150‑1 (MQ=255)
                                                                                                              tGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTg                                       <  1:2939761/150‑1 (MQ=255)
                                                                                                              tGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTg                                       <  2:4018307/150‑1 (MQ=255)
                                                                                                                gACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGc                                     >  1:4291283/1‑150 (MQ=255)
                                                                                                                      gCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAat                               >  2:2621997/1‑150 (MQ=255)
                                                                                                                       ctctGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatg                              >  1:1564831/1‑150 (MQ=255)
                                                                                                                           gCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATgg                          <  2:1045571/150‑1 (MQ=255)
                                                                                                                           gCCACTGGCTTTCCTGAGCAATAATCTCTCTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATgg                          <  2:2416381/150‑1 (MQ=255)
                                                                                                                            ccACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGc                         <  2:3325547/150‑1 (MQ=255)
                                                                                                                              aCTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGNCACGGCGACGGCTGGCAAAAATGATGGCCa                       <  1:4284140/150‑1 (MQ=255)
                                                                                                                                 ggCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGc                    <  2:3880806/150‑1 (MQ=255)
                                                                                                                                  gCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCg                   >  2:1847928/1‑150 (MQ=255)
                                                                                                                                   cTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgt                  >  1:1003534/1‑150 (MQ=255)
                                                                                                                                    tttCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgta                 >  1:3818759/1‑150 (MQ=255)
                                                                                                                                      tCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtagt               >  1:2878614/1‑150 (MQ=255)
                                                                                                                                      tCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtagt               >  2:3274611/1‑150 (MQ=255)
                                                                                                                                       ccTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtagta              >  2:4365869/1‑150 (MQ=255)
                                                                                                                                         tGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCATCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTaaa            >  2:4073340/1‑150 (MQ=255)
                                                                                                                                           aGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGTTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACg          >  2:1832484/1‑150 (MQ=255)
                                                                                                                                             cAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGc        <  1:434305/150‑1 (MQ=255)
                                                                                                                                              aataatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCa       >  1:3946139/1‑150 (MQ=255)
                                                                                                                                              aataatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCa       <  1:4178499/150‑1 (MQ=255)
                                                                                                                                                 aatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCagag    >  1:523184/1‑150 (MQ=255)
                                                                                                                                                  atCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAgt   <  1:4528877/150‑1 (MQ=255)
                                                                                                                                                  atCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAgt   >  1:1015245/1‑150 (MQ=255)
                                                                                                                                                   tctcTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAgtt  <  1:3230870/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTT  >  CP009273/3956017‑3956313

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: