Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,904,230 A→T 6.5% intergenic (+50/+18) PP_2552 → / ← PP_2553 DOPA decarboxylase/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,904,2300AT 6.5% 1052.6 / 43.4 322intergenic (+50/+18)PP_2552/PP_2553DOPA decarboxylase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (177/124);  new base T (5/16);  total (182/140)
Fisher's exact test for biased strand distribution p-value = 2.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

Only 100 of 342 total aligned reads displayed.

CGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCA  >  NC_002947/2904081‑2904378
                                                                                                                                                      |                                                                                                                                                    
cGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaaa                                                                                                                                                      <  2:8647887/150‑1 (MQ=255)
   gCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                                                                                                   >  2:6166556/1‑150 (MQ=255)
       tGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCACGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGccc                                                                                                                                               <  1:6885349/150‑1 (MQ=255)
       tGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGccc                                                                                                                                               <  2:2599959/150‑1 (MQ=255)
        ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTCCCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGcccc                                                                                                                                              <  2:630998/150‑1 (MQ=255)
        ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGcccc                                                                                                                                              <  1:5399580/150‑1 (MQ=255)
          aTGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCtt                                                                                                                                            <  2:3835893/150‑2 (MQ=255)
           tGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGCAGGGGCTGCTTTGCAGCCCCttc                                                                                                                                           <  1:777052/150‑3 (MQ=255)
             gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCaa                                                                                                                                         >  2:7432807/1‑150 (MQ=255)
             gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGGACCGGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCttcgc                                                                                                                                         <  2:7712523/150‑5 (MQ=255)
              tGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt                                                                                                                                        <  1:3307342/150‑1 (MQ=255)
               gCGGGTTTCGATTGGTGCGCTGCCGACCGTGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATc                                                                                                                                       <  1:4632335/150‑1 (MQ=255)
               gCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATc                                                                                                                                       <  2:5614215/150‑1 (MQ=255)
               gCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGGACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATc                                                                                                                                       <  2:659528/150‑1 (MQ=255)
                cGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATTTGCAGCGGCTGTGGGCACGTCTGCTGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAATGGGCTGCAAAGCAGCCCCTGCAATCt                                                                                                                                      <  1:2784694/150‑1 (MQ=255)
                 gggTTTCGATTGGTGCGCTGCCGACCGAGGGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTa                                                                                                                                     <  1:7495166/150‑1 (MQ=255)
                    tttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                                                                                                  >  2:6809358/1‑150 (MQ=255)
                     ttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAg                                                                                                                                 >  2:598312/1‑150 (MQ=255)
                      tCGATTGGTGGGCTGGCGACCGAGCGGGGGGTTGTGCTGCGGCTGTGGGCACGTCTGCTGGTCGTGTTCAAGGGCTGTTGCGTGTAGGGTACCTGTATGTGCGGGGTTGCCCGCGAAGGGGCTGCTTTGCTGCCGCTGCAATGTATCAGt                                                                                                                                <  1:1882095/150‑1 (MQ=255)
                        gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                              >  1:4439504/1‑150 (MQ=255)
                        gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                              >  1:3859915/1‑150 (MQ=255)
                         aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGc                                                                                                                             >  1:4831295/1‑150 (MQ=255)
                           tGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTCGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTc                                                                                                                           <  2:1319514/150‑1 (MQ=255)
                           tGGTGCGCGGCCGACCGAGCGGGGGGAGGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTc                                                                                                                           <  1:6653202/150‑1 (MQ=255)
                            ggTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCa                                                                                                                          >  1:1006008/1‑150 (MQ=255)
                                cgcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGc                                                                                                                      <  1:3001255/150‑1 (MQ=255)
                                cgcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGc                                                                                                                      <  2:3364287/150‑1 (MQ=255)
                                 gcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGcc                                                                                                                     >  1:7564154/1‑150 (MQ=255)
                                     ccgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaa                                                                                                                 <  1:6199139/150‑1 (MQ=255)
                                     ccgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaa                                                                                                                 <  1:1279212/150‑1 (MQ=255)
                                      cgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                                                                                                <  1:4584550/150‑1 (MQ=255)
                                      cgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                                                                                                <  1:8588202/150‑1 (MQ=255)
                                           gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                                                                                                           <  1:157810/150‑1 (MQ=255)
                                                   ggATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACgg                                                                                                   >  2:2388382/1‑150 (MQ=255)
                                                    gATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggc                                                                                                  <  2:872373/150‑1 (MQ=255)
                                                    gATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggc                                                                                                  >  2:4716147/1‑150 (MQ=255)
                                                                cTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAgg                                                                                      <  1:2086004/150‑1 (MQ=255)
                                                                 tgtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                     >  1:651349/1‑150 (MQ=255)
                                                                  ttgatcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCt                                                                                    >  2:8407075/6‑150 (MQ=255)
                                                                     ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                 <  1:1161143/150‑1 (MQ=255)
                                                                     ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                 >  1:5647624/1‑150 (MQ=255)
                                                                     ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                 <  1:2230477/150‑1 (MQ=255)
                                                                      tcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCt                                                                                >  2:3308994/2‑150 (MQ=255)
                                                                      gCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCt                                                                                >  1:7014355/1‑150 (MQ=255)
                                                                        acgtccTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCttt                                                                              >  1:8203430/6‑150 (MQ=255)
                                                                          gTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGa                                                                                                                  <  2:2731186/112‑1 (MQ=255)
                                                                          gTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGa                                                                                                                  >  1:2731186/1‑112 (MQ=255)
                                                                           tCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTcc                                                                           >  2:5161212/1‑150 (MQ=255)
                                                                               gCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                        >  1:3049734/1‑150 (MQ=255)
                                                                                cAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCg                                                                       >  2:5417085/1‑150 (MQ=255)
                                                                                    aCGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCtt                                                                   <  1:8119628/150‑1 (MQ=255)
                                                                                            aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCAGCGGCTTGGCGTCCg                                                           >  1:506871/1‑150 (MQ=255)
                                                                                             aGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCgg                                                          <  2:686856/150‑1 (MQ=255)
                                                                                             aGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTAATAGGCTGCCTTTTCCATCGGCTTGGCGTCCgg                                                          >  2:3426138/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                     >  1:5142970/1‑150 (MQ=255)
                                                                                                   gAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGcc                                                    <  2:1989157/150‑1 (MQ=255)
                                                                                                    aaGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGccc                                                   >  2:817694/1‑150 (MQ=255)
                                                                                                      gCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAg                                                 <  1:1434898/150‑1 (MQ=255)
                                                                                                         tGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTc                                              >  1:6012280/1‑150 (MQ=255)
                                                                                                         tGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTc                                              >  1:835226/1‑150 (MQ=255)
                                                                                                             cGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTGTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTcc                                          >  2:2516175/1‑150 (MQ=255)
                                                                                                                tACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATc                                       >  1:2941806/1‑150 (MQ=255)
                                                                                                                tACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATc                                       <  1:8721509/150‑1 (MQ=255)
                                                                                                                   cTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCACAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGcc                                    >  2:7114349/1‑150 (MQ=255)
                                                                                                                   cTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGcc                                    >  1:4129392/1‑150 (MQ=255)
                                                                                                                    tGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCa                                   >  2:7527054/1‑150 (MQ=255)
                                                                                                                     gTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAg                                  >  2:1080069/1‑150 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTAt                                                                                                                                    <  1:3401615/51‑1 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                 >  1:3627569/1‑150 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                 >  1:7073753/1‑150 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                 >  1:1983617/1‑150 (MQ=255)
                                                                                                                        ttcgCGGGCAAGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCg                               >  1:7261282/4‑150 (MQ=255)
                                                                                                                           gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            >  2:4987498/1‑150 (MQ=255)
                                                                                                                           gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            <  1:4284783/150‑1 (MQ=255)
                                                                                                                           gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            <  2:1299624/150‑1 (MQ=255)
                                                                                                                            cGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                           >  2:6050912/1‑150 (MQ=255)
                                                                                                                                     ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  <  2:8046078/150‑1 (MQ=255)
                                                                                                                                     ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  >  2:6450204/1‑150 (MQ=255)
                                                                                                                                       ttgcaGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGc                >  1:7983674/5‑150 (MQ=255)
                                                                                                                                       gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                           <  2:3137121/139‑1 (MQ=255)
                                                                                                                                       gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGc                >  1:3142149/1‑150 (MQ=255)
                                                                                                                                        cgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCa               <  1:4545187/150‑1 (MQ=255)
                                                                                                                                           aGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg            >  2:1753243/1‑150 (MQ=255)
                                                                                                                                             gggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa          <  2:4320444/150‑1 (MQ=255)
                                                                                                                                             gggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa          <  1:6809358/150‑1 (MQ=255)
                                                                                                                                             gggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa          <  1:5559229/150‑1 (MQ=255)
                                                                                                                                              ggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAg         >  1:907067/1‑150 (MQ=255)
                                                                                                                                               gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg        >  2:216113/1‑150 (MQ=255)
                                                                                                                                               gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg        >  2:1607106/1‑150 (MQ=255)
                                                                                                                                               gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg        >  2:6030278/1‑150 (MQ=255)
                                                                                                                                               gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg        >  2:8797723/1‑150 (MQ=255)
                                                                                                                                                cTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCGTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGc       <  1:33700/150‑1 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt      >  2:7298281/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt      >  2:102649/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt      >  2:4488320/1‑150 (MQ=255)
                                                                                                                                                   cAAAGCAGCCCCTGCAATCTATCAGCGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  >  2:476137/1‑136 (MQ=255)
                                                                                                                                                    aaaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGcc   >  2:3986596/1‑150 (MQ=255)
                                                                                                                                                     aaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            <  2:2527848/124‑1 (MQ=255)
                                                                                                                                                     aaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCa  <  1:2508237/150‑1 (MQ=255)
                                                                                                                                                     aaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCa  <  1:3967857/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCA  >  NC_002947/2904081‑2904378

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: