Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,259,883 G→A 6.8% intergenic (‑112/+32) PP_4634 ← / ← PP_4635 alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,259,8830GA 6.8% 318.2 / 8.8 103intergenic (‑112/+32)PP_4634/PP_4635alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base G (57/39);  new base A (4/3);  total (61/42)
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.94e-01

Only 100 of 103 total aligned reads displayed.

GGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGTCAG  >  NC_002947/5259735‑5260030
                                                                                                                                                    |                                                                                                                                                   
ggCGACATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGc                                                                                                                                                    <  1:3782698/150‑1 (MQ=255)
    acaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCGTGCCcgc                                                                                                                                                <  2:3641828/148‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                             <  2:3202624/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                             <  2:13926/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                             <  1:167834/150‑1 (MQ=255)
       tGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGaa                                                                                                                                             <  1:2489848/150‑1 (MQ=255)
        ggggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAag                                                                                                                                            >  1:5023068/1‑150 (MQ=255)
        ggggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAag                                                                                                                                            >  1:5023024/1‑150 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                           >  1:1784277/1‑150 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                           <  1:3409969/150‑1 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCAACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAaga                                                                                                                                           >  1:755287/1‑150 (MQ=255)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                           <  2:4382773/150‑1 (MQ=17)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                           <  2:1348832/150‑1 (MQ=17)
         gggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                           <  2:1086435/150‑1 (MQ=17)
          ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                          <  1:3319393/150‑1 (MQ=255)
          ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                          <  1:747479/150‑1 (MQ=255)
          ggTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                          <  1:3879889/150‑1 (MQ=255)
            tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        >  2:2336135/1‑150 (MQ=255)
            tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                        <  1:333952/150‑1 (MQ=255)
                      gATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGc                                                                                                                              >  1:3641593/1‑150 (MQ=255)
                         gAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTa                                                                                                                           <  2:3658032/150‑1 (MQ=255)
                              tACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGa                                                                                                                      <  2:1461603/150‑1 (MQ=255)
                                 tCATGGCGCCTGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATGTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCGGGGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                   >  2:2422195/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                   >  1:5034913/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                   >  1:3621559/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                   >  2:3224517/1‑150 (MQ=255)
                                 tCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                                                                   <  1:3824067/150‑1 (MQ=255)
                                   aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACa                                                                                                                 >  1:2043127/1‑150 (MQ=255)
                                       cgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCGGTCCGGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTc                                                                                                             <  1:3349877/150‑1 (MQ=255)
                                       cgcATGTCCGTTTGTTGTGCGCAAAAACTACAAAAGATTTTGTAGCCTGTACTGGCCTCTTCGGGGGTGAACCCGCTCCTACCGGGGTAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTc                                                                                                             <  1:4344492/150‑1 (MQ=38)
                                              ccTGTTTTTGTGCGCAAAAACTAGAAAAGTTTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTGGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAg                                                                                                      <  1:4806616/150‑1 (MQ=38)
                                              ccTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAg                                                                                                      <  2:4336192/150‑1 (MQ=255)
                                                tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTc                                                                                                    >  2:2713691/1‑150 (MQ=255)
                                                tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              <  2:5119364/108‑3 (MQ=255)
                                                tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                              >  1:5119364/1‑106 (MQ=255)
                                                 gttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCg                                                                                                   >  1:1189256/1‑150 (MQ=255)
                                                 gttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGGGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCg                                                                                                   >  1:1189261/1‑150 (MQ=255)
                                                    gttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAAc                                                                                                <  2:1566889/150‑1 (MQ=255)
                                                             cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCg                                                                                       >  1:4978378/1‑150 (MQ=255)
                                                                       aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                             >  1:2805371/1‑150 (MQ=255)
                                                                       aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                             >  2:3970816/1‑150 (MQ=255)
                                                                       aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                             >  2:4198754/1‑150 (MQ=255)
                                                                        aaGATTTTGTAGCCTGTACTGGCCTCTTAGCGGGTGAACCCGCACCTACAGGGGGAGAACCTGCGGGAGCGGGCTTGCCCGCGCAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTc                                                                            >  1:866859/1‑150 (MQ=255)
                                                                         aGATTTGGTAGCCTGTAGTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGGGCGACCCTGTGGGCGCGGGCTTGCCCGCGACGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCt                                                                           <  1:654641/150‑1 (MQ=39)
                                                                           aTTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                                                         <  2:1001024/150‑1 (MQ=255)
                                                                              ttGTAGCCTGTACTGGCCTCTTCGCGGGTGACCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTc                                                                      <  2:4831023/150‑1 (MQ=255)
                                                                                           tGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCa                                                         >  2:371975/1‑150 (MQ=255)
                                                                                             gCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGGGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgc                                                       >  2:4390519/1‑150 (MQ=255)
                                                                                              ccTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg                                                      >  1:3077469/1‑150 (MQ=255)
                                                                                              ccTCTTCGCGGGTGAACCCGCTCCTACAGCGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg                                                      >  1:3734721/1‑150 (MQ=255)
                                                                                                tctTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                    >  2:1941199/1‑150 (MQ=255)
                                                                                                tctTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTATGTCGTAGTCAACGTCAgcgcg                                                    >  1:3346310/1‑150 (MQ=255)
                                                                                                tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                    >  2:366698/1‑150 (MQ=17)
                                                                                                tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                    >  2:2924204/1‑150 (MQ=17)
                                                                                                tctTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcg                                                    >  1:108873/1‑150 (MQ=17)
                                                                                                  ttCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCt                                                  >  2:121699/1‑150 (MQ=255)
                                                                                                     gcgGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGa                                               >  1:656526/1‑150 (MQ=255)
                                                                                                      cgGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGaa                                              <  2:3638605/150‑1 (MQ=255)
                                                                                                        ggTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGc                                            >  1:1050213/1‑150 (MQ=255)
                                                                                                         gTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc                                           >  2:4019164/1‑150 (MQ=255)
                                                                                                         gTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc                                           >  2:4019482/1‑150 (MQ=255)
                                                                                                             aCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAAc                                       <  2:5023024/150‑1 (MQ=255)
                                                                                                             aCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAAc                                       <  2:5023068/150‑1 (MQ=255)
                                                                                                              cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACa                                      >  1:1313923/1‑150 (MQ=255)
                                                                                                              cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACa                                      >  1:3893932/1‑150 (MQ=255)
                                                                                                               ccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAg                                             <  2:1008500/142‑1 (MQ=255)
                                                                                                               ccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAg                                             >  1:1008500/1‑142 (MQ=255)
                                                                                                                 gCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGt                                   >  2:913028/1‑150 (MQ=255)
                                                                                                                   tCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTg                                 >  1:717336/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                               >  2:897331/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                               >  1:856864/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                               >  1:856723/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                               >  1:4938462/1‑150 (MQ=255)
                                                                                                                     cTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAg                               >  1:4938550/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGTGCGCTCGAAGCCGAACAGGTTGAgg                              >  2:2002544/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                              >  2:4145282/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                              >  1:1312958/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                              >  1:1313223/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                              >  2:4646865/1‑150 (MQ=255)
                                                                                                                      tACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAgg                              >  2:4646812/1‑150 (MQ=255)
                                                                                                                        cAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGaa                            >  1:1052957/1‑150 (MQ=255)
                                                                                                                          gggggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACt                          <  2:3445535/150‑1 (MQ=255)
                                                                                                                           ggggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTg                         >  1:3754823/1‑150 (MQ=255)
                                                                                                                             ggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGct                       >  1:466542/1‑150 (MQ=255)
                                                                                                                               tcccccTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttc                     >  1:5120712/5‑150 (MQ=38)
                                                                                                                               agaACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGCTGCCAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttc                     <  2:3864998/150‑1 (MQ=255)
                                                                                                                               agaACCTGGGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGGAGTCAACGGCTGCGCGCTCGAAGCCGCCCAGGTTGAGGAACTGcttc                     >  1:325485/1‑150 (MQ=39)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGTAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                   <  2:2714981/150‑1 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                   >  1:2447783/1‑150 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                   <  1:2046787/150‑1 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                   <  2:5034913/150‑1 (MQ=255)
                                                                                                                                 aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttcta                   <  1:1671638/150‑2 (MQ=255)
                                                                                                                                  aCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGCACTGCTTCTTg                  <  2:1207333/150‑1 (MQ=255)
                                                                                                                                         ggAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCgg           >  1:3310947/1‑150 (MQ=255)
                                                                                                                                             cGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTa       <  2:3621559/150‑1 (MQ=255)
                                                                                                                                             cGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTa       >  2:1923020/1‑150 (MQ=255)
                                                                                                                                               ggCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTagt     <  2:755287/150‑1 (MQ=255)
                                                                                                                                                gCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                                                             <  1:3398174/93‑1 (MQ=255)
                                                                                                                                                gCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                                                             >  2:3398174/1‑93 (MQ=255)
                                                                                                                                                  ttGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTagtcag  >  2:694961/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGTCAG  >  NC_002947/5259735‑5260030

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: