Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,698,343 C→A 16.9% D75Y (GAT→TAT)  ydgJ ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,698,3430CA16.9% 840.3 / 104.3 266D75Y (GAT→TAT) ydgJputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (104/117);  new base A (19/26);  total (123/143)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 268 total aligned reads displayed.

CGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCG  >  CP009273/1698195‑1698492
                                                                                                                                                    |                                                                                                                                                     
  gTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGtt                                                                                                                                                    <  2:2301283/150‑1 (MQ=255)
   ttATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTg                                                                                                                                                   <  1:3286968/150‑1 (MQ=255)
    tATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTgg                                                                                                                                                  >  2:785200/1‑150 (MQ=255)
      tGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGgt                                                                                                                                                <  2:2606023/150‑1 (MQ=255)
      tGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGgt                                                                                                                                                >  2:3400324/1‑150 (MQ=255)
       ggAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGgtg                                                                                                                                               <  1:1510234/150‑1 (MQ=255)
        gAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCCCATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGgtgt                                                                                                                                              <  2:2068952/150‑1 (MQ=255)
          aTACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAg                                                                                                                                            <  1:3837724/150‑1 (MQ=255)
           tACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAgg                                                                                                                                           >  2:1667581/1‑150 (MQ=255)
            acagacagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGa                                                                                                                                          >  1:1952878/1‑150 (MQ=255)
             cagacagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGaa                                                                                                                                         <  2:4333392/150‑1 (MQ=255)
               gacagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATg                                                                                                                                       <  1:4337781/150‑1 (MQ=255)
                acagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGa                                                                                                                                      <  2:3440093/150‑1 (MQ=255)
                acagCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGa                                                                                                                                      <  2:1116213/150‑1 (MQ=255)
                    cacaCGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAAt                                                                                                                                  <  2:3995673/150‑1 (MQ=255)
                       acGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAg                                                                                                                               <  2:4115000/150‑1 (MQ=255)
                       acGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAg                                                                                                                               >  1:4077302/1‑150 (MQ=255)
                       acGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAg                                                                                                                               <  1:429093/150‑1 (MQ=255)
                        cGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAgg                                                                                                                              <  2:473181/150‑1 (MQ=255)
                          ccccAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTc                                                                                                                            <  1:1873400/150‑1 (MQ=255)
                           cccAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCt                                                                                                                           >  2:1315373/1‑150 (MQ=255)
                             cAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTAt                                                                                                                         <  2:1280764/150‑1 (MQ=255)
                             cAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTAt                                                                                                                         <  1:3670303/150‑1 (MQ=255)
                              aGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATg                                                                                                                        >  1:269174/1‑150 (MQ=255)
                               ggCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGt                                                                                                                       >  1:3850705/1‑150 (MQ=255)
                                     ttGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGAt                                                                                                                 <  2:3601427/150‑1 (MQ=255)
                                      tGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATc                                                                                                                <  2:3231945/150‑1 (MQ=255)
                                       gCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCg                                                                                                               <  1:2400519/150‑1 (MQ=255)
                                       gCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTGTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCg                                                                                                               >  2:4266464/1‑150 (MQ=255)
                                            cgcATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTaaa                                                                                                          >  1:2115869/1‑150 (MQ=255)
                                              cATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACa                                                                                                        <  2:2184638/150‑1 (MQ=255)
                                                 ccAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGAt                                                                                                     >  2:2807467/1‑150 (MQ=255)
                                                 ccAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGAt                                                                                                     <  1:1133552/150‑1 (MQ=255)
                                                 ccAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGAt                                                                                                     >  2:893929/1‑150 (MQ=255)
                                                   aGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGc                                                                                                   >  1:2302402/1‑150 (MQ=255)
                                                    gCTCTCGCGCGTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCt                                                                                                  <  2:3895542/150‑1 (MQ=255)
                                                      tctcGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTc                                                                                                <  2:3652348/150‑1 (MQ=255)
                                                        tcGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCgg                                                                                              <  1:4429286/150‑1 (MQ=255)
                                                            gcTTGTGACAGTGTCACGGTAAAGGGTTTAACAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCa                                                                                          >  1:3931595/1‑150 (MQ=255)
                                                             cTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCag                                                                                         >  2:1727834/1‑150 (MQ=255)
                                                             cTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGAGCGTTAAACAGATGCTTCGGCTCag                                                                                         >  1:951975/1‑150 (MQ=255)
                                                              ttGTGACAGTGTCACGGTAAAGTGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCaga                                                                                        <  1:3994024/150‑1 (MQ=255)
                                                               tgtgACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCagag                                                                                       >  1:1011145/1‑150 (MQ=255)
                                                               tgtgACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCagag                                                                                       >  1:4074972/1‑150 (MQ=255)
                                                                  gACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACa                                                                                    >  1:1875478/1‑150 (MQ=255)
                                                                       tgtCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgt                                                                               <  2:2524857/150‑1 (MQ=255)
                                                                         tCACGGTAACGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACATTCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaa                                                                             <  2:1625985/150‑1 (MQ=255)
                                                                          cACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaac                                                                            <  1:3551756/150‑1 (MQ=255)
                                                                           aCGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaacc                                                                           <  2:310037/150‑1 (MQ=255)
                                                                            cGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaaccg                                                                          >  1:6254/1‑150 (MQ=255)
                                                                            cGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaaccg                                                                          >  2:1184774/1‑150 (MQ=255)
                                                                             ggTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaaccgt                                                                         >  2:4396843/1‑150 (MQ=255)
                                                                             ggTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACaaccgtaaccgt                                                                         >  1:1673308/1‑150 (MQ=255)
                                                                              gTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGtt                                                                        <  2:2414559/150‑1 (MQ=255)
                                                                                aaaGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTgg                                                                      >  1:2800758/1‑150 (MQ=255)
                                                                                  aGGGTTTATCAACGACCACATGTTTACCCGCCTCACGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTACCCGTTGGcc                                                                    <  1:4372011/150‑1 (MQ=255)
                                                                                  aGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGcc                                                                    <  1:785200/150‑1 (MQ=255)
                                                                                  aGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGcc                                                                    <  1:3017831/150‑1 (MQ=255)
                                                                                  aGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGcc                                                                    <  2:577099/150‑1 (MQ=255)
                                                                                   gggTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCa                                                                   <  2:4321334/150‑1 (MQ=255)
                                                                                     gTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGt                                                                 <  2:1318969/150‑1 (MQ=255)
                                                                                         aTCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGc                                                             >  1:3532153/1‑150 (MQ=255)
                                                                                         aTCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGc                                                             >  1:1209631/1‑150 (MQ=255)
                                                                                            aaCGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCttt                                                          <  1:426326/150‑1 (MQ=255)
                                                                                             aCGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCtttt                                                         <  1:974433/150‑1 (MQ=255)
                                                                                             aCGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCtttt                                                         <  1:2586719/150‑1 (MQ=255)
                                                                                             aCGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCtttt                                                         >  1:4186090/1‑150 (MQ=255)
                                                                                              cGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTa                                                        >  1:2034616/1‑150 (MQ=255)
                                                                                               gACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTAc                                                       >  2:1942550/1‑150 (MQ=255)
                                                                                               gACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTAc                                                       >  2:122183/1‑150 (MQ=255)
                                                                                                      tGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGtt                                                <  1:4623618/150‑1 (MQ=255)
                                                                                                       gTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGttt                                               <  1:2662178/150‑1 (MQ=255)
                                                                                                       gTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGttt                                               >  1:4605465/1‑150 (MQ=255)
                                                                                                         ttACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTtca                                             <  1:378665/150‑1 (MQ=255)
                                                                                                          tACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTtcat                                            >  1:2578707/1‑150 (MQ=255)
                                                                                                            cccGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCGGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTtcatca                                          <  2:2826504/150‑1 (MQ=255)
                                                                                                            cccGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTtcatca                                          <  2:2310657/150‑1 (MQ=255)
                                                                                                              cGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCAct                                        >  1:1693829/1‑150 (MQ=255)
                                                                                                                 cTCAAGCGCCGCTTTGGCTAACGGGAATTGGGTATCGTTGGGTGTATGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCActgct                                     <  2:4119004/150‑1 (MQ=255)
                                                                                                                    aaGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGga                                  >  2:2612804/1‑150 (MQ=255)
                                                                                                                    aaGCGCCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGga                                  >  2:2854760/1‑150 (MQ=255)
                                                                                                                      gcgcCGCTTTGGCTAACGGGAATTGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGgaga                                <  1:228333/150‑1 (MQ=255)
                                                                                                                      gcgcCGCTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGgaga                                <  2:2708194/150‑1 (MQ=255)
                                                                                                                        gccgcTTTGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATTCTTCGGCTCAGAGACAAACTTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGAtt                              >  1:2948919/1‑150 (MQ=255)
                                                                                                                             tttGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGc                         >  1:1702420/1‑150 (MQ=255)
                                                                                                                             tttGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGc                         >  1:4469070/1‑150 (MQ=255)
                                                                                                                             tttGGCTAACGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGc                         >  1:2925491/1‑150 (MQ=255)
                                                                                                                                 gCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGt                     <  2:1472431/150‑1 (MQ=255)
                                                                                                                                  cTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGtt                    <  1:3028644/150‑1 (MQ=255)
                                                                                                                                  cTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGtt                    <  1:781467/150‑1 (MQ=255)
                                                                                                                                    aaCGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTcc                  >  1:3031670/1‑150 (MQ=255)
                                                                                                                                      cGGGAAATGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTg                <  1:3260628/150‑1 (MQ=255)
                                                                                                                                       gggAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCAGc               >  1:4272198/1‑150 (MQ=255)
                                                                                                                                         gAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGccc             >  1:1618402/1‑150 (MQ=255)
                                                                                                                                          aaaTGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTAGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCg            <  1:242775/150‑1 (MQ=255)
                                                                                                                                          aaaTGGGTATAGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCg            <  2:2907666/150‑1 (MQ=255)
                                                                                                                                             tGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGc         <  1:541199/150‑1 (MQ=255)
                                                                                                                                               ggTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGt       <  1:4004509/150‑1 (MQ=255)
                                                                                                                                                   tCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTATTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGccc   >  1:250819/1‑150 (MQ=255)
                                                                                                                                                    cGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCg  <  1:3400324/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGGTAAAGGGTTTATCAACGACCACATGTTTACCCGCCTCAAGCGCCGCTTTGGCTAACGGGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCG  >  CP009273/1698195‑1698492

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: