Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,618,693 G→T 100% E74* (GAA→TAA)  ppx → exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,618,6930GT96.2% 549.4 / 10.7 183E74* (GAA→TAA) ppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base T (85/91);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 183 total aligned reads displayed.

GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGA  >  CP009273/2618548‑2618835
                                                                                                                                                 |                                                                                                                                              
ggTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTaa                                                                                                                                              <  1:607948/148‑1 (MQ=255)
ggTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTaa                                                                                                                                              <  1:3530536/148‑1 (MQ=255)
   cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACgg                                                                                                                                           >  1:1070761/1‑148 (MQ=255)
     tAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCt                                                                                                                                         >  2:700073/1‑148 (MQ=255)
     tAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCt                                                                                                                                         <  2:1932841/148‑1 (MQ=255)
       gCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTAc                                                                                                                                       >  1:43525/1‑148 (MQ=255)
        cccGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACa                                                                                                                                      >  2:1063119/1‑148 (MQ=255)
        cccGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGAAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACa                                                                                                                                      <  1:2443487/148‑1 (MQ=255)
            tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAggg                                                                                                                                  <  1:1771380/148‑1 (MQ=255)
             gtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGt                                                                                                                                 >  1:1563085/1‑148 (MQ=255)
              tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtt                                                                                                                                <  2:952360/148‑1 (MQ=255)
                 tAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtc                                                                                                                             >  1:1644366/1‑148 (MQ=255)
                   gATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtctc                                                                                                                           >  1:1549911/1‑148 (MQ=255)
                      ggTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTg                                                                                                                        >  1:2519308/1‑148 (MQ=255)
                       gTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGc                                                                                                                       >  1:200204/1‑148 (MQ=255)
                        tGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGcc                                                                                                                      <  1:50809/148‑1 (MQ=255)
                        tGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGcc                                                                                                                      <  1:50601/148‑1 (MQ=255)
                            aTGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCg                                                                                                                  <  2:1478042/148‑1 (MQ=255)
                             tGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGt                                                                                                                 <  1:1861674/148‑1 (MQ=255)
                             tGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGt                                                                                                                 <  2:1754428/148‑1 (MQ=255)
                              gCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                                                                                <  1:1220089/148‑1 (MQ=255)
                              gCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTACTGCCAGCGTc                                                                                                                >  2:1465145/1‑148 (MQ=255)
                               cAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCt                                                                                                               <  2:3317694/148‑1 (MQ=255)
                               cAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCt                                                                                                               <  2:2862389/148‑1 (MQ=255)
                                  attattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGta                                                                                                            >  2:560932/1‑148 (MQ=255)
                                   ttattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                                                                                           <  1:2482476/148‑1 (MQ=255)
                                     attGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                                                                                         <  2:3472535/148‑1 (MQ=255)
                                     attGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                                                                                         <  2:3337417/148‑1 (MQ=255)
                                     attGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                                                                                         <  1:119959/148‑1 (MQ=255)
                                      ttGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                        <  1:2883126/148‑1 (MQ=255)
                                       tGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGtt                                                                                                       >  1:2320736/1‑148 (MQ=255)
                                        ggCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTg                                                                                                      <  1:1376191/148‑1 (MQ=255)
                                          ccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGt                                                                                                    <  1:1128077/148‑1 (MQ=255)
                                            gccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAc                                                                                                  <  1:1014636/148‑1 (MQ=255)
                                              cTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAccc                                                                                                >  1:1310924/1‑148 (MQ=255)
                                                gAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAt                                                                                              <  2:1460315/148‑1 (MQ=255)
                                                   aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                                                                                           >  2:3151135/1‑148 (MQ=255)
                                                   aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                                                                                           >  1:2848660/1‑148 (MQ=255)
                                                   aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                                                                                           >  1:2725787/1‑148 (MQ=255)
                                                   aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                                                                                           <  1:1736669/148‑1 (MQ=255)
                                                      gCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTg                                                                                        >  1:1397067/1‑148 (MQ=255)
                                                       cGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                                                                                       >  2:2643955/1‑148 (MQ=255)
                                                       cGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                                                                                       >  2:1979401/1‑148 (MQ=255)
                                                          gTGCATCTGGCGGACGGCCTGGGGCCAGAGAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                                                                                    <  2:2473217/148‑1 (MQ=255)
                                                                      gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACg                                                                        >  2:786549/1‑148 (MQ=255)
                                                                      gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACg                                                                        >  1:217724/1‑148 (MQ=255)
                                                                      gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACg                                                                        >  1:2016597/1‑148 (MQ=255)
                                                                       aCGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc                                                                       >  1:3332452/1‑148 (MQ=255)
                                                                       aCGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc                                                                       >  1:670443/1‑148 (MQ=255)
                                                                        cGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGcc                                                                      <  1:571271/148‑1 (MQ=255)
                                                                         ggCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCa                                                                     <  2:1099865/148‑1 (MQ=255)
                                                                           ccTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACt                                                                   >  2:1786958/1‑148 (MQ=255)
                                                                             tGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa                                                                 <  2:1841283/148‑1 (MQ=255)
                                                                                ggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACtt                                                              >  1:53667/1‑148 (MQ=255)
                                                                                 gCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGCTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACttt                                                             <  2:1183796/148‑1 (MQ=255)
                                                                                  ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTc                                                            <  2:1070761/148‑1 (MQ=255)
                                                                                  ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTc                                                            >  1:1187522/1‑148 (MQ=255)
                                                                                     gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGa                                                         <  1:2973062/148‑1 (MQ=255)
                                                                                     gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGa                                                         <  1:2973071/148‑1 (MQ=255)
                                                                                        aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGTCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                                                      >  1:965086/1‑148 (MQ=255)
                                                                                        aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                                                      >  1:331008/1‑148 (MQ=255)
                                                                                        aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGTCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                                                      >  2:3286487/1‑148 (MQ=255)
                                                                                         atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                     <  2:111667/148‑1 (MQ=255)
                                                                                         atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                     <  2:206862/148‑1 (MQ=255)
                                                                                          tatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgc                                                    >  1:3407989/1‑148 (MQ=255)
                                                                                           atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                   >  2:1583275/1‑148 (MQ=255)
                                                                                           atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                   >  2:90034/1‑148 (MQ=255)
                                                                                           atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                   >  1:1149552/1‑148 (MQ=255)
                                                                                            tGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcgc                                                  <  2:859931/148‑1 (MQ=255)
                                                                                              ttGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCgg                                                >  2:43671/1‑148 (MQ=255)
                                                                                              ttGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCgg                                                >  1:3102450/1‑148 (MQ=255)
                                                                                                 aGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaaa                                             >  1:354630/1‑148 (MQ=255)
                                                                                                   tGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAg                                           >  2:3005209/1‑148 (MQ=255)
                                                                                                     aaGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGt                                         >  1:2215486/1‑148 (MQ=255)
                                                                                                     aaGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGt                                         <  2:258854/148‑1 (MQ=255)
                                                                                                     aaGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGt                                         >  2:3276489/1‑148 (MQ=255)
                                                                                                           cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                   >  1:3405398/1‑148 (MQ=255)
                                                                                                              tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                >  1:3500014/1‑148 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                              <  2:2503092/148‑1 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                              >  1:723805/1‑148 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                              >  1:617476/1‑148 (MQ=255)
                                                                                                                 cgcgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAccc                             >  2:2569768/1‑148 (MQ=255)
                                                                                                                  gcgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCg                            <  2:2263026/148‑1 (MQ=255)
                                                                                                                          tAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAat                    >  2:2534720/1‑148 (MQ=255)
                                                                                                                           aaaCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAatt                   <  2:2320736/148‑1 (MQ=255)
                                                                                                                             aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattat                 >  2:2970185/1‑148 (MQ=255)
                                                                                                                             aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattat                 <  2:42045/148‑1 (MQ=255)
                                                                                                                              ctgtctgtCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                >  2:3080044/1‑148 (MQ=255)
                                                                                                                              ctgtctgtCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                >  1:2449415/1‑148 (MQ=255)
                                                                                                                               tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt               <  2:2254522/148‑1 (MQ=255)
                                                                                                                                     tcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTa         <  2:653879/148‑1 (MQ=255)
                                                                                                                                       gCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAAt       <  2:2010921/148‑1 (MQ=255)
                                                                                                                                        cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATg      <  2:853295/148‑1 (MQ=255)
                                                                                                                                         tGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATga     <  2:1257752/148‑1 (MQ=255)
                                                                                                                                         tGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATga     <  2:3513274/148‑1 (MQ=255)
                                                                                                                                         tGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATga     <  2:3512383/148‑1 (MQ=255)
                                                                                                                                            ttGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaaga  <  1:3151135/148‑1 (MQ=255)
                                                                                                                                            ttGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaaga  >  2:2122098/1‑148 (MQ=255)
                                                                                                                                            ttGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaaga  <  1:90034/148‑1 (MQ=255)
                                                                                                                                            ttGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaaga  <  1:1583275/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGA  >  CP009273/2618548‑2618835

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