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,6930GT98.3% 926.1 / 4.5 285E74* (GAA→TAA) ppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (124/156);  total (126/159)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 285 total aligned reads displayed.

ATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGC  >  CP009273/2618546‑2618838
                                                                                                                                                   |                                                                                                                                                 
aTGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTaa                                                                                                                                                 >  1:4479087/1‑150 (MQ=255)
 tGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCCGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAAc                                                                                                                                                <  2:3937250/150‑1 (MQ=255)
   gTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACgg                                                                                                                                              <  1:832178/150‑1 (MQ=255)
    tCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGc                                                                                                                                             <  2:4186225/150‑1 (MQ=255)
     cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCt                                                                                                                                            >  1:4539686/1‑150 (MQ=255)
     cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCt                                                                                                                                            <  2:3856186/150‑1 (MQ=255)
          cccGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAg                                                                                                                                       <  2:4328247/150‑1 (MQ=255)
          cacGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAG‑GGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAgg                                                                                                                                      <  2:2962492/148‑1 (MQ=255)
            cGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAggg                                                                                                                                     >  2:4103122/1‑150 (MQ=255)
              tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtt                                                                                                                                   <  2:4669049/150‑1 (MQ=255)
              tgtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtt                                                                                                                                   <  2:3960723/150‑1 (MQ=255)
               gtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGttt                                                                                                                                  <  1:3805474/150‑1 (MQ=255)
                tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtttt                                                                                                                                 >  2:745415/1‑150 (MQ=255)
                tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtttt                                                                                                                                 >  1:4898096/1‑150 (MQ=255)
                tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtttt                                                                                                                                 >  2:223107/1‑150 (MQ=255)
                 ggTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtc                                                                                                                                <  2:1684341/150‑1 (MQ=255)
                 ggTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtc                                                                                                                                <  1:1471316/150‑1 (MQ=255)
                 ggTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtc                                                                                                                                >  2:4393612/1‑150 (MQ=255)
                     gATGGTGCCATGCAGATTATTGGTCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCt                                                                                                                            <  1:1743217/150‑1 (MQ=255)
                      aTGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTg                                                                                                                           <  2:16262/150‑1 (MQ=255)
                        ggTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGcc                                                                                                                         <  1:4770700/150‑1 (MQ=255)
                            ccATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCg                                                                                                                     >  1:143361/1‑150 (MQ=255)
                                   gATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                                                                                              <  1:3539603/150‑1 (MQ=255)
                                   gATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                                                                                              <  2:1263060/150‑1 (MQ=255)
                                    attattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                                                             <  1:3128093/150‑1 (MQ=255)
                                      tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                           >  1:3745751/1‑150 (MQ=255)
                                      tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                           >  1:5152847/1‑150 (MQ=255)
                                         tGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                                                                                        <  2:1536313/150‑1 (MQ=255)
                                         tGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTgg                                                                                                        >  1:1088535/1‑150 (MQ=255)
                                          ggCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGt                                                                                                       >  2:2460639/1‑150 (MQ=255)
                                          ggCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGt                                                                                                       <  2:765956/150‑1 (MQ=255)
                                           gccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTa                                                                                                      >  1:4355893/1‑150 (MQ=255)
                                            ccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAc                                                                                                     <  2:2380835/150‑1 (MQ=255)
                                             cgccTGTCACCGCGGGTGCTTCTGGCGGACGGCCTGGCGCGAGATAATATGTTGAGTGAAGAGGCAATGACGCGCTGTTGAATCTGTCTTTCGCGGTTTGACTAACGGCTACTAGGGTTTTCTCCTGCCAGCGTCTGTTTAGTTGGTaac                                                                                                    <  2:2638751/150‑3 (MQ=255)
                                             cgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAcc                                                                                                    >  1:4378941/1‑150 (MQ=255)
                                               ccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCa                                                                                                  >  1:1913112/1‑150 (MQ=255)
                                                cTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGTTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAt                                                                                                 >  1:1551126/1‑150 (MQ=255)
                                                   aaaCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACg                                                                                              >  2:844227/1‑150 (MQ=255)
                                                     aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCt                                                                                            <  1:3576809/150‑1 (MQ=255)
                                                     aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCt                                                                                            <  1:4073327/150‑1 (MQ=255)
                                                     aCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCt                                                                                            >  1:2191558/1‑150 (MQ=255)
                                                          gggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTc                                                                                       >  2:1855551/1‑150 (MQ=255)
                                                           ggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCa                                                                                      <  1:3629044/150‑1 (MQ=255)
                                                              gCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgc                                                                                   <  2:1516127/150‑1 (MQ=255)
                                                                aTCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGgcgc                                                                                 >  1:907506/1‑150 (MQ=255)
                                                                 tCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGACTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCt                                                                                <  1:1775207/150‑1 (MQ=255)
                                                                   tGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGa                                                                              <  1:1868935/150‑1 (MQ=255)
                                                                   tGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGa                                                                              <  2:3806249/150‑1 (MQ=255)
                                                                    ggCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGaa                                                                             <  1:947677/150‑1 (MQ=255)
                                                                     gCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAAc                                                                            <  1:204815/150‑1 (MQ=255)
                                                                        gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGcc                                                                         <  1:648401/150‑1 (MQ=255)
                                                                        gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGcc                                                                         <  1:1998541/150‑1 (MQ=255)
                                                                          cGGCCTGGGGCCAGATAATATGTTGAGTCAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCAc                                                                       >  1:3868441/1‑150 (MQ=255)
                                                                              cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGAc                                                                   >  2:4555628/1‑150 (MQ=255)
                                                                              cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGAc                                                                   >  1:5061767/1‑150 (MQ=255)
                                                                              cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAACCATACGCTGCGTCAGGCGCTGAACGCCACTGAc                                                                   >  2:585175/1‑150 (MQ=255)
                                                                                ggggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACtt                                                                 >  2:1020522/1‑150 (MQ=255)
                                                                                  ggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTc                                                               <  2:2124810/150‑1 (MQ=255)
                                                                                    ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTg                                                             >  2:3517519/1‑150 (MQ=255)
                                                                                    ccAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTg                                                             >  2:924779/1‑150 (MQ=255)
                                                                                      aGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaa                                                           <  2:205479/150‑1 (MQ=255)
                                                                                        ataataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAc                                                         >  1:3318890/1‑150 (MQ=255)
                                                                                           atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                      >  2:2167773/1‑150 (MQ=255)
                                                                                           atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcg                                                      <  2:1406416/150‑1 (MQ=255)
                                                                                            tatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcgc                                                     >  2:1601024/1‑150 (MQ=255)
                                                                                             atGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTGGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgcgcg                                                    <  2:137486/150‑1 (MQ=255)
                                                                                              tGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCgg                                                   <  1:2252496/150‑1 (MQ=255)
                                                                                                   aGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAg                                              <  2:2701345/150‑1 (MQ=255)
                                                                                                   aGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAg                                              >  2:4835607/1‑150 (MQ=255)
                                                                                                      gAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTc                                           >  2:658820/1‑150 (MQ=255)
                                                                                                      gAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTc                                           <  1:2500197/150‑1 (MQ=255)
                                                                                                         gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCAtt                                        <  2:3546071/150‑1 (MQ=255)
                                                                                                          agGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTc                                       <  1:639794/150‑1 (MQ=255)
                                                                                                           ggCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                      <  2:3496011/150‑1 (MQ=255)
                                                                                                             cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCt                                    <  2:3752720/150‑1 (MQ=255)
                                                                                                             cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCt                                    <  2:569812/150‑1 (MQ=255)
                                                                                                              aaTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                   <  2:3330629/150‑1 (MQ=255)
                                                                                                                    gcgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGAt                             <  2:516109/150‑1 (MQ=255)
                                                                                                                    gcgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGAt                             <  2:5037667/150‑1 (MQ=255)
                                                                                                                      gcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTg                           <  1:4810969/150‑1 (MQ=255)
                                                                                                                      gcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTg                           <  2:4884838/150‑1 (MQ=255)
                                                                                                                      gcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTg                           >  2:2890661/1‑150 (MQ=255)
                                                                                                                       cgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGa                          <  1:157514/150‑1 (MQ=255)
                                                                                                                       cgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGa                          <  1:4631198/150‑1 (MQ=255)
                                                                                                                       cgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGa                          <  2:202117/150‑1 (MQ=255)
                                                                                                                        ggTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATGGaa                         >  1:2631094/1‑150 (MQ=255)
                                                                                                                              aaCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                   <  1:4296319/150‑1 (MQ=255)
                                                                                                                              aaCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                   >  1:4519689/1‑150 (MQ=255)
                                                                                                                               aCTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                  >  1:2138079/1‑150 (MQ=255)
                                                                                                                                 tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc                <  1:844660/150‑1 (MQ=255)
                                                                                                                                 tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc                <  2:1168216/150‑1 (MQ=255)
                                                                                                                                  gtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCg               >  2:4583131/1‑150 (MQ=255)
                                                                                                                                    ctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGt             <  1:2771132/150‑1 (MQ=255)
                                                                                                                                       tcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAAt          <  1:3798306/150‑1 (MQ=255)
                                                                                                                                       tcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTGTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAAt          <  2:3905598/150‑1 (MQ=255)
                                                                                                                                        cGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATg         <  1:3896754/150‑1 (MQ=255)
                                                                                                                                        cGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATg         >  2:3049596/1‑150 (MQ=255)
                                                                                                                                          cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaa       >  1:591605/1‑150 (MQ=255)
                                                                                                                                          cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaa       <  1:1166519/150‑1 (MQ=255)
                                                                                                                                           tGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaag      <  2:1991767/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGC  >  CP009273/2618546‑2618838

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