Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,904,229 A→T 6.1% intergenic (+49/+19) 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,2290AT 6.1% 1085.3 / 37.6 329intergenic (+49/+19)PP_2552/PP_2553DOPA decarboxylase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (177/129);  new base T (5/16);  total (183/146)
Fisher's exact test for biased strand distribution p-value = 2.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

GACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCA  >  NC_002947/2904079‑2904378
                                                                                                                                                       |                                                                                                                                                     
  cGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaaa                                                                                                                                                      <  2:8647887/150‑1 (MQ=255)
     gCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                                                                                                   <  2:2679430/150‑1 (MQ=255)
     gCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                                                                                                   >  2:6166556/1‑150 (MQ=255)
       ggTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagcagc                                                                                                                                                 <  2:2459884/150‑1 (MQ=255)
        gTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGcc                                                                                                                                                <  2:109224/150‑1 (MQ=255)
         tGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCACGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGccc                                                                                                                                               <  1:6885349/150‑1 (MQ=255)
          ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTCCCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGcccc                                                                                                                                              <  2:630998/150‑1 (MQ=255)
          ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGcccc                                                                                                                                              <  1:5399580/150‑1 (MQ=255)
          ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGcccc                                                                                                                                              <  2:6539983/150‑1 (MQ=255)
            aTGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCtt                                                                                                                                            <  2:3835893/150‑2 (MQ=255)
               gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCaa                                                                                                                                         >  2:7432807/1‑150 (MQ=255)
                tGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt                                                                                                                                        <  1:3307342/150‑1 (MQ=255)
                 gCGGGTTTCGATTGGTGCGCTGCCGACCGTGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATc                                                                                                                                       <  1:4632335/150‑1 (MQ=255)
                  cGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATTTGCAGCGGCTGTGGGCACGTCTGCTGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAATGGGCTGCAAAGCAGCCCCTGCAATCt                                                                                                                                      <  1:2784694/150‑1 (MQ=255)
                   gggTTTCGATTGGTGCGCTGCCGACCGAGGGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTa                                                                                                                                     <  1:7495166/150‑1 (MQ=255)
                      tttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTGCGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                                                                                                  <  2:2564343/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)
                          gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                              <  2:4159183/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)
                                   gcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGcc                                                                                                                     >  1:7564154/1‑150 (MQ=255)
                                   gcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGcc                                                                                                                     >  1:5769536/1‑150 (MQ=255)
                                      gCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGa                                                                                                                  <  2:1392286/150‑1 (MQ=255)
                                       ccgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaa                                                                                                                 <  1:1279212/150‑1 (MQ=255)
                                        cgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                                                                                                <  1:8588202/150‑1 (MQ=255)
                                         gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAg                                                                                                               >  2:3247387/1‑150 (MQ=255)
                                             gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                                                                                                           <  1:157810/150‑1 (MQ=255)
                                                      gATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggc                                                                                                  <  2:872373/150‑1 (MQ=255)
                                                       aTGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggct                                                                                                 <  1:4022347/150‑1 (MQ=255)
                                                       aTGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggct                                                                                                 >  1:4962325/1‑150 (MQ=255)
                                                              gCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCa                                                                                          >  1:3071978/1‑150 (MQ=255)
                                                              gCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCa                                                                                          >  2:4536446/1‑150 (MQ=255)
                                                                   tgtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                     >  1:651349/1‑150 (MQ=255)
                                                                    ttgatcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAACGCTACCTTTACGGCTGGCTTCATAGGCt                                                                                    >  2:5545395/6‑150 (MQ=255)
                                                                    ttgatcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCt                                                                                    >  2:8407075/6‑150 (MQ=255)
                                                                       ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                 <  1:1161143/150‑1 (MQ=255)
                                                                       ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                 <  1:2230477/150‑1 (MQ=255)
                                                                        tcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCt                                                                                >  2:3308994/2‑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)
                                                                                  cAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCg                                                                       >  2:5417085/1‑150 (MQ=255)
                                                                                   aGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCgg                                                                      >  1:2134264/1‑150 (MQ=255)
                                                                                      aCGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCtt                                                                   <  1:8119853/150‑1 (MQ=255)
                                                                                       cGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTg                                                                  <  1:8398831/150‑1 (MQ=255)
                                                                                           aTCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGt                                                              >  1:3703519/1‑150 (MQ=255)
                                                                                            tCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTc                                                             >  1:6823512/1‑150 (MQ=255)
                                                                                              aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCAGCGGCTTGGCGTCCg                                                           >  1:506871/1‑150 (MQ=255)
                                                                                               aGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCgg                                                          >  2:7711663/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)
                                                                                                      aaGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGccc                                                   >  2:817694/1‑150 (MQ=255)
                                                                                                        gCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAg                                                 <  1:1434898/150‑1 (MQ=255)
                                                                                                           tGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTc                                              >  1:6012280/1‑150 (MQ=255)
                                                                                                               cGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTGTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTcc                                          >  2:2516175/1‑150 (MQ=255)
                                                                                                                  tACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATc                                       <  1:8721509/150‑1 (MQ=255)
                                                                                                                  tACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATc                                       >  1:2941806/1‑150 (MQ=255)
                                                                                                                    ccTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGc                                     <  1:7027912/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)
                                                                                                                        tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTAt                                                                                                                                    <  1:3401615/51‑1 (MQ=255)
                                                                                                                        tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                 >  1:7073753/1‑150 (MQ=255)
                                                                                                                        tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                 >  2:6011749/1‑150 (MQ=255)
                                                                                                                        tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                 >  1:1983617/1‑150 (MQ=255)
                                                                                                                            aGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGt                             <  2:2007706/150‑1 (MQ=255)
                                                                                                                             gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            <  1:4284783/150‑1 (MQ=255)
                                                                                                                             gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            <  2:1299624/150‑1 (MQ=255)
                                                                                                                             gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                            >  2:4987498/1‑150 (MQ=255)
                                                                                                                              cGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                           >  2:6050912/1‑150 (MQ=255)
                                                                                                                              cGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                           <  1:7350276/150‑1 (MQ=255)
                                                                                                                                 gCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGtt                        >  2:1705875/1‑150 (MQ=255)
                                                                                                                                       ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  >  1:910312/1‑150 (MQ=255)
                                                                                                                                       ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  >  2:3385684/1‑150 (MQ=255)
                                                                                                                                       ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  >  2:6450204/1‑150 (MQ=255)
                                                                                                                                       ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  <  2:8046078/150‑1 (MQ=255)
                                                                                                                                         ttgcaGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGc                >  1:7983674/5‑150 (MQ=255)
                                                                                                                                         gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGc                >  1:3142149/1‑150 (MQ=255)
                                                                                                                                          cgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCa               <  1:4545187/150‑1 (MQ=255)
                                                                                                                                             aGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg            >  2:1753243/1‑150 (MQ=255)
                                                                                                                                              ggggCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                           <  1:6239730/134‑1 (MQ=255)
                                                                                                                                               gggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa          <  1:5559229/150‑1 (MQ=255)
                                                                                                                                                 gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg        >  2:8797723/1‑150 (MQ=255)
                                                                                                                                                 gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg        >  2:216113/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)
                                                                                                                                                   tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCAATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt      >  1:7802581/1‑150 (MQ=255)
                                                                                                                                                     cAAAGCAGCCCCTGCAATCTATCAGCGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                  >  2:476137/1‑136 (MQ=255)
                                                                                                                                                      aaaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGcc   >  1:4598346/1‑150 (MQ=255)
                                                                                                                                                      aaaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGcc   >  2:3986596/1‑150 (MQ=255)
                                                                                                                                                       aaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCa  <  1:2508237/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCA  >  NC_002947/2904079‑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: