Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,020,480 A→C 5.2% Y185S (TAT→TCT)  helD → DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,020,4800AC 5.2% 832.2 / 31.7 268Y185S (TAT→TCT) helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (121/133);  new base C (8/6);  total (129/139)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 270 total aligned reads displayed.

TTTGTCGGCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGTTGCT  >  CP009273/1020331‑1020627
                                                                                                                                                     |                                                                                                                                                   
tttGTCGGCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTa                                                                                                                                                     >  2:410499/1‑150 (MQ=255)
   gTCGGCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGc                                                                                                                                                  <  2:2920396/150‑1 (MQ=255)
   gTCGGCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGc                                                                                                                                                  >  1:3155016/1‑150 (MQ=255)
    tCGGCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTCTGCg                                                                                                                                                 >  2:1171249/1‑150 (MQ=255)
      ggCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGa                                                                                                                                               >  2:858215/1‑150 (MQ=255)
      ggCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGa                                                                                                                                               >  2:856725/1‑150 (MQ=255)
       gCGTTGCCGTTGCCGGTTAACCGACTGGCAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATAGTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAt                                                                                                                                              >  2:1116520/1‑150 (MQ=255)
       gCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAATATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAt                                                                                                                                              >  2:4281524/1‑150 (MQ=255)
       gCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAt                                                                                                                                              >  2:4281797/1‑150 (MQ=255)
        cgttgccgttgcCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAtt                                                                                                                                             <  2:572410/150‑1 (MQ=255)
        cgttgccgttgcCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAtt                                                                                                                                             <  2:2855434/150‑1 (MQ=255)
         gttgccgttgccGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAttt                                                                                                                                            >  2:2162095/1‑150 (MQ=255)
          ttgccgttgccgGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGAtttt                                                                                                                                           <  1:3551003/150‑1 (MQ=255)
             ccgttgccgGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTc                                                                                                                                        >  2:3914202/1‑150 (MQ=255)
             ccgttgccgGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTc                                                                                                                                        >  2:324154/1‑150 (MQ=255)
              cgttgccgGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTcc                                                                                                                                       >  2:198928/1‑150 (MQ=255)
                    cgGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAgg                                                                                                                                 <  2:417189/150‑1 (MQ=255)
                     ggTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGt                                                                                                                                >  1:2428776/1‑150 (MQ=255)
                         aaCCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAg                                                                                                                            <  2:1237646/150‑1 (MQ=255)
                           ccGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTc                                                                                                                          <  1:2458196/150‑1 (MQ=255)
                             gACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCtt                                                                                                                        <  1:1981887/150‑1 (MQ=255)
                              aCTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTc                                                                                                                       <  2:1132832/150‑1 (MQ=255)
                               cTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCa                                                                                                                      <  1:4252755/150‑1 (MQ=255)
                                  gaagaaTTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCg                                                                                                                   <  1:390615/150‑1 (MQ=255)
                                    agaaTTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTAACCGCt                                                                                                                 >  2:616438/1‑150 (MQ=255)
                                           gATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCg                                                                                                          >  2:4050133/1‑150 (MQ=255)
                                           gATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCg                                                                                                          >  2:3048777/1‑150 (MQ=255)
                                           gATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCg                                                                                                          >  2:3049431/1‑150 (MQ=255)
                                                 tGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAg                                                                                                    <  2:17145/150‑1 (MQ=255)
                                                     gTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGccc                                                                                                <  1:4058729/150‑1 (MQ=255)
                                                          gcgtggcgtAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCa                                                                                           <  1:3548613/150‑1 (MQ=255)
                                                             tggcgtAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTc                                                                                        <  1:3618372/150‑1 (MQ=255)
                                                                 gtAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTa                                                                                    >  2:1279348/1‑150 (MQ=255)
                                                                 gtAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTa                                                                                    <  1:1033174/150‑1 (MQ=255)
                                                                 gtAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTa                                                                                    <  1:3310453/150‑1 (MQ=255)
                                                                 gtAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTa                                                                                    <  2:665176/150‑1 (MQ=255)
                                                                   aaaTGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAAt                                                                                  <  1:834227/150‑1 (MQ=255)
                                                                      tGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGc                                                                               >  2:3177083/1‑150 (MQ=255)
                                                                       gTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCg                                                                              <  2:1784379/150‑1 (MQ=255)
                                                                        tCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGa                                                                             <  1:1273871/150‑1 (MQ=255)
                                                                         cAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAg                                                                            <  2:728051/150‑1 (MQ=255)
                                                                            gCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTCTGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCAt                                                                         >  1:3420944/1‑150 (MQ=255)
                                                                               tGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtct                                                                      >  2:1473475/1‑150 (MQ=255)
                                                                               tGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtct                                                                      <  1:2740142/150‑1 (MQ=255)
                                                                               tGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtct                                                                      <  1:3048777/150‑1 (MQ=255)
                                                                                ggTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtctg                                                                     <  2:3924709/150‑2 (MQ=255)
                                                                                ggTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtctc                                                                     <  1:3576193/150‑1 (MQ=255)
                                                                                ggTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtctc                                                                     <  1:2192192/150‑1 (MQ=255)
                                                                                ggTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTACACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATtctc                                                                     <  1:815/150‑1 (MQ=255)
                                                                                  ttGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTg                                                                   >  1:3995125/1‑150 (MQ=255)
                                                                                    gAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGtt                                                                 >  2:2610972/1‑150 (MQ=255)
                                                                                      aaGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAg                                                               <  1:985208/150‑1 (MQ=255)
                                                                                      aaGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAg                                                               <  1:3798157/150‑1 (MQ=255)
                                                                                      aaGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAg                                                               <  1:985686/150‑1 (MQ=255)
                                                                                            ttGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTgg                                                         <  2:3511172/150‑1 (MQ=255)
                                                                                              gAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCa                                                       >  1:710156/1‑150 (MQ=255)
                                                                                               aaaGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAg                                                      >  1:2927213/1‑150 (MQ=255)
                                                                                                  gcgcTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTg                                                   <  1:236167/150‑1 (MQ=255)
                                                                                                    gcTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCa                                                 >  1:134954/1‑150 (MQ=255)
                                                                                                      tCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAgg                                               <  2:3559027/150‑1 (MQ=255)
                                                                                                         gTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAg                                            <  2:3656408/150‑1 (MQ=255)
                                                                                                          ttGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCGGGTGCAGGAAGc                                           >  1:3132280/1‑150 (MQ=255)
                                                                                                           tGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCg                                          >  2:3594651/1‑150 (MQ=255)
                                                                                                           tGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCg                                          >  2:2080296/1‑150 (MQ=255)
                                                                                                           tGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCg                                          >  2:2080422/1‑150 (MQ=255)
                                                                                                            gcagcaTAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCgg                                         <  2:4492433/150‑1 (MQ=255)
                                                                                                             cagcaTAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGTATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGa                                        >  1:4475255/1‑150 (MQ=255)
                                                                                                               gcaTAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGaaa                                      >  2:4497823/1‑150 (MQ=255)
                                                                                                                  tAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAAc                                   >  1:3503960/1‑150 (MQ=255)
                                                                                                                   aaCCAGGCGTATACCGAAGCCATGCTTACCGAGTCTGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACg                                  >  2:1894531/1‑150 (MQ=255)
                                                                                                                   aaCCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACg                                  >  1:1098141/1‑150 (MQ=255)
                                                                                                                    aCCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGt                                 >  1:530932/1‑150 (MQ=255)
                                                                                                                      cAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCg                               <  1:2108688/150‑1 (MQ=255)
                                                                                                                        ggCGTATACCGAAGCCATGCTTACCGAGTCTGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGt                             <  1:2614157/150‑1 (MQ=255)
                                                                                                                        ggCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGt                             <  2:3369437/150‑1 (MQ=255)
                                                                                                                              taCCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCtggt                       <  1:813777/150‑1 (MQ=255)
                                                                                                                                ccGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCtggtgg                     >  2:968693/1‑150 (MQ=255)
                                                                                                                                ccGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCtggtgg                     <  1:1091169/150‑1 (MQ=255)
                                                                                                                                ccGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCtggtgg                     >  2:2269210/1‑150 (MQ=255)
                                                                                                                                ccGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCtggtgg                     >  1:3488015/1‑150 (MQ=255)
                                                                                                                                ccGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCtggtgg                     >  1:3488205/1‑150 (MQ=255)
                                                                                                                                  gAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCTAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGcc                   >  1:1893794/1‑150 (MQ=255)
                                                                                                                                   aaGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGccc                  <  2:1156193/150‑1 (MQ=255)
                                                                                                                                   aaGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGccc                  <  2:3427860/150‑1 (MQ=255)
                                                                                                                                   aaGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGccc                  <  1:3292875/150‑1 (MQ=255)
                                                                                                                                    aGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCg                 <  2:751855/150‑1 (MQ=255)
                                                                                                                                     gCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGt                <  1:4266654/150‑1 (MQ=255)
                                                                                                                                       cATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGc              <  2:1530890/150‑1 (MQ=255)
                                                                                                                                       cATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGc              <  2:1530883/150‑1 (MQ=255)
                                                                                                                                        aTGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCa             >  1:985520/1‑150 (MQ=255)
                                                                                                                                          gCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAgg           <  1:600882/150‑1 (MQ=255)
                                                                                                                                            ttACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCt         <  1:2680106/150‑1 (MQ=255)
                                                                                                                                             tACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTg        <  1:3602244/150‑1 (MQ=255)
                                                                                                                                                cGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGtt     >  1:1816933/1‑150 (MQ=255)
                                                                                                                                                cGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGtt     >  1:2904985/1‑150 (MQ=255)
                                                                                                                                                cGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGtt     >  2:1175529/1‑150 (MQ=255)
                                                                                                                                                 gAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGTTg    >  1:1776721/1‑150 (MQ=255)
                                                                                                                                                 gAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGTTg    >  1:3433311/1‑150 (MQ=255)
                                                                                                                                                  aGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGTTGc   <  1:2279802/150‑1 (MQ=255)
                                                                                                                                                  aGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGTTGc   <  1:1513111/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTTGTCGGCGTTGCCGTTGCCGGTTAACCGACTGGAAGAATTCGATAACTGCCGTGAGGCGTGGCGTAAATGTCAGGCCTGGTTGAAAGATATTGAAAGCGCTCGGTTGCAGCATAACCAGGCGTATACCGAAGCCATGCTTACCGAGTATGCGGATTTTTTCCGCCAGGTCGAGTCTTCACCGCTGAATCCGGCGCAGGCCCGGGCAGTCGTTAATGGCGAGCATTCTCTGTTAGTGCTGGCAGGTGCAGGAAGCGGAAAAACGTCGGTGCTGGTGGCCCGTGCAGGCTGGTTGCT  >  CP009273/1020331‑1020627

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: