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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,618,6930GT63.7% 180.9 / inf 273E74* (GAA→TAA) ppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (48/51);  new base T (77/97);  total (125/148)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 273 total aligned reads displayed.

CATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGCC  >  CP009273/2618545‑2618839
                                                                                                                                                    |                                                                                                                                                  
caTGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGa                                                                                                                                                   >  1:548402/1‑150 (MQ=255)
 aTGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTaa                                                                                                                                                  >  2:3752893/1‑150 (MQ=255)
   ggTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACg                                                                                                                                                >  2:3000412/1‑150 (MQ=255)
        tAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTAc                                                                                                                                           >  1:3731648/1‑150 (MQ=255)
             cGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAggg                                                                                                                                      <  1:2219166/150‑1 (MQ=255)
                 tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtttt                                                                                                                                  <  1:2788274/150‑1 (MQ=255)
                   gTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtct                                                                                                                                <  2:1515557/150‑1 (MQ=255)
                    tAGATGGTTCCATGAAGATTATTGGCCTCCTGTAACAGCGGGTGCATCTGGCGGACGGCCTGGGGACAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTtctc                                                                                                                               <  2:4317405/150‑1 (MQ=255)
                    tAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTtctc                                                                                                                               <  2:2072114/150‑1 (MQ=255)
                       aTGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTg                                                                                                                            <  2:3457338/150‑1 (MQ=255)
                         ggTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGcc                                                                                                                          >  2:2831311/1‑150 (MQ=255)
                          gTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCa                                                                                                                         >  1:456148/1‑150 (MQ=255)
                           tGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAg                                                                                                                        <  1:129081/150‑1 (MQ=255)
                            gCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGc                                                                                                                       <  1:2073779/150‑1 (MQ=255)
                               aTGCAGATTTTTGGCCGACTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGTGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCTGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCATCGTc                                                                                                                    <  2:284062/150‑1 (MQ=255)
                               aTGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                                                                                    <  1:826052/150‑1 (MQ=255)
                               aTGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                                                                                    <  1:3108958/150‑1 (MQ=255)
                                tGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCt                                                                                                                   <  2:1132207/150‑1 (MQ=255)
                                  cAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGt                                                                                                                 <  1:4507705/150‑1 (MQ=255)
                                  cAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGt                                                                                                                 >  2:2131856/1‑150 (MQ=255)
                                    gATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                                                                                               >  2:22358/1‑150 (MQ=255)
                                       tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                            >  1:4011122/1‑150 (MQ=255)
                                       tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                            >  1:325494/1‑150 (MQ=255)
                                       tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                            >  1:2326912/1‑150 (MQ=255)
                                            gccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTa                                                                                                       <  2:4177530/150‑1 (MQ=255)
                                             ccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAc                                                                                                      <  1:3447963/150‑1 (MQ=255)
                                                       cAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTg                                                                                            <  1:244877/150‑1 (MQ=255)
                                                        aGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGc                                                                                           <  2:4062443/150‑1 (MQ=255)
                                                         gCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCg                                                                                          <  2:2309566/150‑1 (MQ=255)
                                                         gCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCg                                                                                          <  2:2309597/150‑1 (MQ=255)
                                                             gTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAg                                                                                      <  2:1596694/150‑1 (MQ=255)
                                                             gTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAg                                                                                      <  1:2915782/150‑1 (MQ=255)
                                                               gCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgc                                                                                    <  2:1608643/150‑1 (MQ=255)
                                                                  tCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCt                                                                                 <  2:896141/150‑1 (MQ=255)
                                                                     ggCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                                                              >  2:3218533/1‑150 (MQ=255)
                                                                     ggCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                                                              >  2:3292567/1‑150 (MQ=255)
                                                                     ggCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                                                              >  1:1970079/1‑150 (MQ=255)
                                                                     ggCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                                                              >  2:2376532/1‑150 (MQ=255)
                                                                        ggACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc                                                                           >  1:572576/1‑150 (MQ=255)
                                                                         gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGcc                                                                          <  2:4373277/150‑1 (MQ=255)
                                                                            ggCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACt                                                                       >  1:2375811/1‑150 (MQ=255)
                                                                             gCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTg                                                                      <  1:4347058/150‑1 (MQ=255)
                                                                             gCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTg                                                                      <  1:4347108/150‑1 (MQ=255)
                                                                             gCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTg                                                                      <  1:3441817/150‑1 (MQ=255)
                                                                              ccTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa                                                                     <  2:1819388/150‑1 (MQ=255)
                                                                              ccTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa                                                                     <  2:1819343/150‑1 (MQ=255)
                                                                              ccTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa                                                                     <  1:1997873/150‑1 (MQ=255)
                                                                               cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGAc                                                                    >  2:687523/1‑150 (MQ=255)
                                                                                 ggggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACtt                                                                  >  1:3802538/1‑150 (MQ=255)
                                                                                  gggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACttt                                                                 >  1:1076696/1‑150 (MQ=255)
                                                                                   ggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTc                                                                <  2:2092412/150‑1 (MQ=255)
                                                                                      cAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGa                                                             >  1:578568/1‑150 (MQ=255)
                                                                                        gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaaa                                                           >  1:1315336/1‑150 (MQ=255)
                                                                                        gATAATATGTTGAGTGAAGAGGCAATGACGCGCGGATTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaaa                                                           <  2:253932/150‑1 (MQ=255)
                                                                                         ataataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                                                          >  1:388697/1‑150 (MQ=255)
                                                                                          taataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                         <  2:3487487/150‑1 (MQ=255)
                                                                                           aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgc                                                        <  2:3403066/150‑1 (MQ=255)
                                                                                           aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgc                                                        <  1:1789254/150‑1 (MQ=255)
                                                                                           aataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgc                                                        >  1:3905039/1‑150 (MQ=255)
                                                                                            atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                       <  2:3622886/150‑1 (MQ=255)
                                                                                             tatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcgc                                                      >  1:1157696/1‑150 (MQ=255)
                                                                                               tGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCgg                                                    >  1:3510446/1‑150 (MQ=255)
                                                                                                  tGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaaa                                                 <  2:468854/150‑1 (MQ=255)
                                                                                                   gAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTTCCCATAGGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaaaa                                                <  1:2940360/150‑1 (MQ=255)
                                                                                                      tGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTGTCTGAAACGCGCGGAAAAGGt                                             >  1:3059085/1‑150 (MQ=255)
                                                                                                          gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCAtt                                         <  1:306676/150‑1 (MQ=255)
                                                                                                            ggCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                       <  1:3597275/150‑1 (MQ=255)
                                                                                                             gCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTccc                                      <  2:3549098/150‑1 (MQ=255)
                                                                                                              cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCt                                     >  2:3110257/1‑150 (MQ=255)
                                                                                                              cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCt                                     <  2:4257974/150‑1 (MQ=255)
                                                                                                               aaTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTTGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                    <  2:1291573/150‑1 (MQ=255)
                                                                                                                aTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAc                                   >  2:3925013/1‑150 (MQ=255)
                                                                                                                aTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAc                                   >  2:3920109/1‑150 (MQ=255)
                                                                                                                aTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAc                                   >  2:2715182/1‑150 (MQ=255)
                                                                                                                 tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                                  >  1:664944/1‑150 (MQ=255)
                                                                                                                  gACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGGAAAGGTCATTCCCTAccc                                 >  1:1957924/1‑150 (MQ=255)
                                                                                                                  gACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAccc                                 >  2:738731/1‑150 (MQ=255)
                                                                                                                  gACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAccc                                 >  1:677369/1‑150 (MQ=255)
                                                                                                                    ctcgatgTTTAAACTGACAGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCAGTATAGTTGGTACCCAAACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACACGa                               >  2:1417720/7‑150 (MQ=255)
                                                                                                                       gcgGTTTAAACTGTCTGTGGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTg                            <  1:4275357/150‑1 (MQ=255)
                                                                                                                       gcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTg                            >  2:1065562/1‑150 (MQ=255)
                                                                                                                        cgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGa                           <  2:2607553/150‑1 (MQ=255)
                                                                                                                        cgGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGa                           >  1:4170716/1‑150 (MQ=255)
                                                                                                                         ggTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                          >  2:3317242/1‑150 (MQ=255)
                                                                                                                          gTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaaa                         <  2:2043309/150‑1 (MQ=255)
                                                                                                                          gTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAAGCGCGGAAAAGGTCATTCCCTACCCGATTGaaa                         >  2:4148814/1‑150 (MQ=255)
                                                                                                                           tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAat                        >  1:1687392/1‑150 (MQ=255)
                                                                                                                            ttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAatt                       >  1:2141165/1‑150 (MQ=255)
                                                                                                                               aaCTGTCTGTCGCTGTTTGACTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                    >  2:1785111/1‑150 (MQ=255)
                                                                                                                                aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                   <  2:1003114/150‑1 (MQ=255)
                                                                                                                                aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                   >  1:1594829/1‑150 (MQ=255)
                                                                                                                                aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                   <  1:2582622/150‑1 (MQ=255)
                                                                                                                                  tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc                 <  1:3591601/150‑1 (MQ=255)
                                                                                                                                  tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc                 <  2:572576/150‑1 (MQ=255)
                                                                                                                                  tgtctgtcGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc                 >  1:4420778/1‑150 (MQ=255)
                                                                                                                                   gtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCg                <  1:3218121/150‑1 (MQ=255)
                                                                                                                                    tctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCgg               <  1:2831311/150‑1 (MQ=255)
                                                                                                                                        tcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAAt           <  1:3398904/150‑1 (MQ=255)
                                                                                                                                           cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaa        >  1:2543523/1‑150 (MQ=255)
                                                                                                                                              tttGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaagaa     >  2:2372087/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGCC  >  CP009273/2618545‑2618839

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