Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,215,252 G→T 14.2% E499D (GAG→GAT metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,215,2520GT14.2% 965.6 / 89.7 286E499D (GAG→GATmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (142/102);  new base T (18/23);  total (161/125)
Fisher's exact test for biased strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACAT  >  CP009273/4215103‑4215401
                                                                                                                                                     |                                                                                                                                                     
ggCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCgag                                                                                                                                                       <  1:2438732/150‑1 (MQ=255)
   aaaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAtt                                                                                                                                                    <  2:1307843/150‑1 (MQ=255)
   aaaGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGAtt                                                                                                                                                    <  2:1107431/150‑1 (MQ=255)
      ggCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGc                                                                                                                                                 >  2:3182720/1‑150 (MQ=255)
        cATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCg                                                                                                                                               >  1:3272582/1‑150 (MQ=255)
         aTTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGt                                                                                                                                              >  1:1017526/1‑150 (MQ=255)
         aTTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGt                                                                                                                                              >  2:3432726/1‑150 (MQ=255)
         aTTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGt                                                                                                                                              >  1:1668999/1‑150 (MQ=255)
          ttGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTc                                                                                                                                             <  2:3404437/150‑1 (MQ=255)
          ttGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTc                                                                                                                                             >  2:3510359/1‑150 (MQ=255)
          ttGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTc                                                                                                                                             >  2:2073038/1‑150 (MQ=255)
            gTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCgg                                                                                                                                           >  1:4429620/1‑150 (MQ=255)
              tAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                                         >  2:718431/1‑150 (MQ=255)
              tAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGc                                                                                                                                         >  2:3264047/1‑150 (MQ=255)
                aCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGt                                                                                                                                       <  1:4586995/150‑1 (MQ=255)
                aCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGt                                                                                                                                       <  2:1512706/150‑1 (MQ=255)
                aCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGt                                                                                                                                       >  1:81609/1‑150 (MQ=255)
                aCTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGt                                                                                                                                       <  2:293887/150‑1 (MQ=255)
                        tcGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATc                                                                                                                               <  1:3822185/150‑1 (MQ=255)
                         cGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATcc                                                                                                                              <  1:1674801/150‑1 (MQ=255)
                          gATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCt                                                                                                                             <  1:3766040/150‑1 (MQ=255)
                           aTGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTc                                                                                                                            >  1:2929501/1‑150 (MQ=255)
                           aTGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTc                                                                                                                            >  1:2187630/1‑150 (MQ=255)
                              aaaGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCAcc                                                                                                                         >  1:4297603/1‑150 (MQ=255)
                              aaaGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCAcc                                                                                                                         >  1:3191983/1‑150 (MQ=255)
                                 gagGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGaa                                                                                                                      <  2:467850/150‑1 (MQ=255)
                                  agGGCGTCGATGCCTTTATCCATCACGCGTAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAAAAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGatg                                                                                                                     <  1:180731/150‑3 (MQ=255)
                                  agGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAag                                                                                                                     >  1:2112163/1‑150 (MQ=255)
                                     gCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAgg                                                                                                                  >  1:197226/1‑150 (MQ=255)
                                     gCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAgg                                                                                                                  <  1:627407/150‑1 (MQ=255)
                                      cGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGt                                                                                                                 >  1:1314411/1‑150 (MQ=255)
                                        tCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTg                                                                                                               <  1:352453/150‑1 (MQ=255)
                                        tCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTg                                                                                                               >  1:933970/1‑150 (MQ=255)
                                         cGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTgg                                                                                                              <  1:125931/150‑1 (MQ=255)
                                            tGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCtt                                                                                                           >  2:2434230/1‑150 (MQ=255)
                                            tGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCtt                                                                                                           <  2:2529651/150‑1 (MQ=255)
                                               cTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTccc                                                                                                        >  1:1847048/1‑150 (MQ=255)
                                                tttATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCg                                                                                                       <  2:2771784/150‑1 (MQ=255)
                                                tttATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCg                                                                                                       <  2:2771463/150‑1 (MQ=255)
                                                  tATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGcc                                                                                                     >  1:4343397/1‑150 (MQ=255)
                                                    tCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCag                                                                                                   <  1:991356/150‑1 (MQ=255)
                                                      cATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                                                                 >  2:3719170/1‑150 (MQ=255)
                                                      cATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaa                                                                                                 >  1:1175640/1‑150 (MQ=255)
                                                       aTCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCagaag                                                                                                >  2:1376740/1‑150 (MQ=255)
                                                            gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                                           <  1:3173922/150‑1 (MQ=255)
                                                            gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                                           <  2:3647905/150‑1 (MQ=255)
                                                            gcgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatc                                                                                           <  2:4457318/150‑1 (MQ=255)
                                                             cgAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatca                                                                                          <  2:189429/150‑1 (MQ=255)
                                                              gAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatcat                                                                                         >  1:4035331/1‑150 (MQ=255)
                                                              gAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATatcat                                                                                         >  2:2314645/1‑150 (MQ=255)
                                                                aaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCt                                                                                       >  2:447166/1‑150 (MQ=255)
                                                                aaTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCt                                                                                       >  1:4248197/1‑150 (MQ=255)
                                                                 aTTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACCGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCtt                                                                                      >  2:2204197/1‑150 (MQ=255)
                                                                   tgttgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCg                                                                                    >  2:2248692/1‑150 (MQ=255)
                                                                      tgCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGAcc                                                                                 >  1:2038420/1‑150 (MQ=255)
                                                                            gCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACa                                                                           >  2:3484982/1‑150 (MQ=255)
                                                                            gCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACa                                                                           >  1:3604000/1‑150 (MQ=255)
                                                                               aCGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCt                                                                        <  1:308395/150‑1 (MQ=255)
                                                                               aCGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCt                                                                        >  2:4107783/1‑150 (MQ=255)
                                                                                cGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCtt                                                                       <  1:2854586/150‑1 (MQ=255)
                                                                                cGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCtt                                                                       <  1:3992916/150‑1 (MQ=255)
                                                                                 ggTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTc                                                                      >  2:607464/1‑150 (MQ=255)
                                                                                  gTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTcg                                                                     <  2:3584784/150‑1 (MQ=255)
                                                                                    gCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTcgcg                                                                   >  2:2318984/1‑150 (MQ=255)
                                                                                      ggCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGt                                                                 >  2:3154590/1‑150 (MQ=255)
                                                                                       gCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTc                                                                >  1:712137/1‑150 (MQ=255)
                                                                                         aGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGTTCGc                                                              >  1:1870365/1‑150 (MQ=255)
                                                                                             gtggtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACt                                                          <  1:1337770/150‑1 (MQ=255)
                                                                                                gtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                                       >  2:1278002/1‑150 (MQ=255)
                                                                                                gtAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGc                                                       <  2:2836198/150‑1 (MQ=255)
                                                                                                 tAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCa                                                      <  1:1147027/150‑1 (MQ=255)
                                                                                                 tAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCa                                                      <  1:2804560/150‑1 (MQ=255)
                                                                                                  aaTGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCAt                                                     >  2:2958334/1‑150 (MQ=255)
                                                                                                        cTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAaga                                               <  2:81609/150‑1 (MQ=255)
                                                                                                         tttGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAagag                                              <  2:3881/150‑1 (MQ=255)
                                                                                                               gAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCacaac                                        <  2:4334484/150‑1 (MQ=255)
                                                                                                               gAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCacaac                                        >  1:3686946/1‑150 (MQ=255)
                                                                                                                  cAGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACaacaac                                     >  1:1441111/1‑150 (MQ=255)
                                                                                                                   aGGGACAGGCCGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACt                                    <  1:3529102/150‑1 (MQ=255)
                                                                                                                   aGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACt                                    <  1:137275/150‑1 (MQ=255)
                                                                                                                   aGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACt                                    >  2:1154174/1‑150 (MQ=255)
                                                                                                                   aGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACt                                    >  1:2515982/1‑150 (MQ=255)
                                                                                                                   aGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACt                                    <  1:3330834/150‑1 (MQ=255)
                                                                                                                    gggACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTa                                   <  1:671863/150‑1 (MQ=255)
                                                                                                                     ggACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAc                                  <  2:1838666/150‑1 (MQ=255)
                                                                                                                      gACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAcg                                 >  1:1588699/1‑150 (MQ=255)
                                                                                                                       aCAGGCCGATACTCGCGCACGGAAAATCGACATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTACTCGGTCGCAACTGCCATTGAAGAGCACAACAACTAcgc                                >  1:4291190/1‑150 (MQ=255)
                                                                                                                          ggCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCa                             >  2:3003595/1‑150 (MQ=255)
                                                                                                                            ccGATACTCGCGCACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAgg                           <  2:1314447/150‑1 (MQ=255)
                                                                                                                            ccGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAgg                           <  1:1769648/150‑1 (MQ=255)
                                                                                                                                  cTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTa                     >  2:3053010/1‑150 (MQ=255)
                                                                                                                                   tCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTAt                    >  1:2317233/1‑150 (MQ=255)
                                                                                                                                   tCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTAt                    <  2:875998/150‑1 (MQ=255)
                                                                                                                                   tCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTAAGCGCAGGACTTTAt                    >  1:4528209/1‑150 (MQ=255)
                                                                                                                                    cgcgCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAACTCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATc                   <  2:4273501/150‑1 (MQ=255)
                                                                                                                                     gcgcACGGAAAATCGATATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCg                  <  2:3988948/150‑1 (MQ=255)
                                                                                                                                       gcACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGgc                >  1:1711266/1‑150 (MQ=255)
                                                                                                                                          cGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGcgcg             <  2:3367417/150‑1 (MQ=255)
                                                                                                                                             aaaaTCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCgtgt          <  2:2187630/150‑1 (MQ=255)
                                                                                                                                                  cGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCGTCCCGCCAGAAGATATCATCGTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGa     >  2:1645048/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAAAGGCATTGTTAACTCTATCTCGATGAAAGAGGGCGTCGATGCCTTTATCCATCACGCGAAATTGTTGCGTCGCTACGGTGCGGCAGTGGTGGTAATGGCCTTTGACGAACAGGGACAGGCCGATACTCGCGCACGGAAAATCGAGATTTGCCGTCGGGCGTACAAAATCCTCACCGAAGAGGTTGGCTTCCCGCCAGAAGATATCATCTTCGACCCAAACATCTTCGCGGTCGCAACTGGCATTGAAGAGCACAACAACTACGCGCAGGACTTTATCGGCGCGTGTGAAGACAT  >  CP009273/4215103‑4215401

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: