Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,901,202 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,901,2020AT 6.1% 546.5 / 17.0 164intergenic (+49/+19)PP_2552/PP_2553DOPA decarboxylase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (94/60);  new base T (4/6);  total (98/66)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 168 total aligned reads displayed.

ACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCC  >  NC_002947/2901053‑2901350
                                                                                                                                                     |                                                                                                                                                    
aCGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaa                                                                                                                                                      <  1:3093600/150‑1 (MQ=255)
  gggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAag                                                                                                                                                    <  1:1886293/150‑1 (MQ=255)
   ggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagc                                                                                                                                                   <  2:2800499/150‑1 (MQ=255)
    gCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                                                                                                  >  1:3177960/1‑150 (MQ=255)
      ggTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagcagc                                                                                                                                                <  1:3718640/150‑1 (MQ=255)
         ggATGGTGCGGGTTTCGATTTGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGcccc                                                                                                                                             <  1:2514617/150‑1 (MQ=255)
          gATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGCAGGGGCTGCTTTGCAGCCCCt                                                                                                                                            <  2:3756928/150‑1 (MQ=255)
          gATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCt                                                                                                                                            >  2:1325456/1‑150 (MQ=255)
            tGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCttc                                                                                                                                          <  2:2200555/150‑3 (MQ=255)
            tGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGc                                                                                                                                          >  1:3481116/1‑150 (MQ=255)
            tGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGc                                                                                                                                          <  2:3828335/150‑1 (MQ=255)
               tGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt                                                                                                                                       >  1:3242213/1‑150 (MQ=255)
                  gggTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTa                                                                                                                                    >  2:1669591/1‑150 (MQ=255)
                  gggTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTa                                                                                                                                    >  2:3995791/1‑150 (MQ=255)
                    gTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATc                                                                                                                                  >  2:2670415/1‑150 (MQ=255)
                     tttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                                                                                                 >  2:3496862/1‑150 (MQ=255)
                         gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                             >  2:4775071/1‑150 (MQ=255)
                          aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                                                                                                 >  2:4481063/1‑145 (MQ=255)
                          aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                                                                                                 <  1:4481063/145‑1 (MQ=255)
                          aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGc                                                                                                                            <  1:4087318/150‑1 (MQ=255)
                               tGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCAcc                                                                                                                       >  1:3682749/1‑150 (MQ=255)
                                   cTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCt                                                                                                                   <  2:4329434/150‑1 (MQ=255)
                                      ccgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaa                                                                                                                >  1:1322292/1‑150 (MQ=255)
                                       cgaccgcGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                                                                                               <  1:4967158/150‑1 (MQ=255)
                                          ccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTg                                                                                                                                           <  2:5047405/119‑1 (MQ=255)
                                            gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                                                                                                          <  1:4467323/150‑1 (MQ=255)
                                            gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                                                                                                          <  1:1223873/150‑1 (MQ=255)
                                            gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                                                                                                          <  1:1122052/150‑1 (MQ=255)
                                                  ggggATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTAc                                                                                                    >  1:3599433/1‑150 (MQ=255)
                                                      aTGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggct                                                                                                <  1:55870/150‑1 (MQ=255)
                                                       tgtgCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctg                                                                                               >  1:1774130/1‑150 (MQ=255)
                                                       cctgCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctg                                                                                               >  1:2003282/3‑150 (MQ=255)
                                                        gtgCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctgg                                                                                              >  2:1980004/1‑150 (MQ=255)
                                                         tgCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGACGCCTGTACGGTACCTGTATGAGCGGGCTGGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc                                                                                             <  1:2289712/150‑1 (MQ=255)
                                                         tgCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc                                                                                             >  2:1508749/1‑150 (MQ=255)
                                                         tgCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc                                                                                             >  1:2038385/1‑150 (MQ=255)
                                                                gCTGTGGGCTCGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGCGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAg                                                                                      <  1:455572/150‑1 (MQ=255)
                                                                gCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAg                                                                                      <  2:5128915/150‑1 (MQ=255)
                                                                  tgtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                    <  2:4508454/150‑1 (MQ=255)
                                                                  tgtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                    <  1:2673616/150‑1 (MQ=255)
                                                                  cttgatcACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTAGCTTCATAGGc                                                                                    >  2:1299893/7‑150 (MQ=255)
                                                                   gtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCt                                                                                   <  2:2632568/150‑1 (MQ=255)
                                                                   gtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCt                                                                                   >  2:5095383/1‑150 (MQ=255)
                                                                     gggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGc                                                                                 >  1:235239/1‑150 (MQ=255)
                                                                      ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                >  1:4750428/1‑150 (MQ=255)
                                                                      ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                >  2:4898971/1‑150 (MQ=255)
                                                                         aCGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCttt                                                                             >  1:266444/1‑150 (MQ=255)
                                                                           gTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTc                                                                           >  2:3216689/1‑150 (MQ=255)
                                                                              tGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCAt                                                                        >  2:4918177/1‑150 (MQ=255)
                                                                                cAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCg                                                                      >  2:4519356/1‑150 (MQ=255)
                                                                                   gACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCt                                                                   >  1:2802359/1‑150 (MQ=255)
                                                                                    aCGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCtt                                                                  >  1:685954/1‑150 (MQ=255)
                                                                                     cGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTg                                                                 <  2:2230092/150‑1 (MQ=255)
                                                                                      gTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTgg                                                                >  2:10384/1‑150 (MQ=255)
                                                                                       tGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGACGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGc                                                               <  1:1854167/150‑1 (MQ=255)
                                                                                       tGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGc                                                               >  2:2016760/1‑150 (MQ=255)
                                                                                       tGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGc                                                               >  2:2016629/1‑150 (MQ=255)
                                                                                        gATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCg                                                              >  2:1768317/1‑150 (MQ=255)
                                                                                           cAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTcc                                                           >  1:4390023/1‑150 (MQ=255)
                                                                                           cAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTcc                                                           <  1:472716/150‑1 (MQ=255)
                                                                                            aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGTGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCg                                                          >  2:1215815/1‑150 (MQ=255)
                                                                                            aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCg                                                          >  1:2863364/1‑150 (MQ=255)
                                                                                             aGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCgg                                                         <  2:3639671/150‑1 (MQ=255)
                                                                                              gggCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCggg                                                        <  2:1831721/150‑1 (MQ=255)
                                                                                               ggCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGt                                                       >  1:807238/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGACTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                    >  2:3045575/1‑150 (MQ=255)
                                                                                                    aaGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGccc                                                  <  1:4932302/150‑1 (MQ=255)
                                                                                                     aGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCa                                                 >  2:2535759/1‑150 (MQ=255)
                                                                                                      gCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAg                                                <  2:3786809/150‑1 (MQ=255)
                                                                                                        cTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGt                                              >  2:549763/1‑150 (MQ=255)
                                                                                                        cTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGt                                              >  1:3473004/1‑150 (MQ=255)
                                                                                                         tGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTc                                             >  2:2202890/1‑150 (MQ=255)
                                                                                                          gTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCt                                            <  2:2486476/150‑1 (MQ=255)
                                                                                                              ggTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCa                                        >  1:3800345/1‑150 (MQ=255)
                                                                                                                    ccccttcgCGGGCAAGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCa                                  >  1:3349038/8‑150 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                >  2:1408920/1‑150 (MQ=255)
                                                                                                                       cttcgCGGGCAAGCCCGCGAAGGGGCTGCAAAGCAGCCCCGGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATAGCCAGGc                               >  1:2157620/5‑150 (MQ=255)
                                                                                                                        ttcgCGGGCAAGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCg                              >  2:1907619/4‑150 (MQ=255)
                                                                                                                        tGAGCGGGCTTGCGCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTGCGGGTTGCCCAGGTCTTCCATCGCCTGGCg                              >  2:2153567/1‑150 (MQ=255)
                                                                                                                          aGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGt                            <  2:1701064/150‑1 (MQ=255)
                                                                                                                           gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                           >  1:2972165/1‑150 (MQ=255)
                                                                                                                            cGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                          >  1:811318/1‑150 (MQ=255)
                                                                                                                            cGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTgg                          <  2:3242213/150‑1 (MQ=255)
                                                                                                                             gggCTTGCCCGCGAAGGGGCTGCAAAGCAGCCACTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAACTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGTGTTGCCCAGGTCTTCCATCGCCAGGCGGTggg                         <  1:3440394/150‑1 (MQ=255)
                                                                                                                              ggcaagCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGt                        >  1:2223875/6‑150 (MQ=255)
                                                                                                                                 ttGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTggaagac                     <  1:5095383/150‑6 (MQ=255)
                                                                                                                                 ttGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTc                     >  2:4383274/1‑150 (MQ=255)
                                                                                                                                 ttGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTc                     >  1:5244377/1‑150 (MQ=255)
                                                                                                                                 ttGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTc                     >  1:2111025/1‑150 (MQ=255)
                                                                                                                                   gCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                                                                                               <  2:3446870/58‑1 (MQ=255)
                                                                                                                                       gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGc               >  2:2769306/1‑150 (MQ=255)
                                                                                                                                           aGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAg                                 >  1:4788141/1‑128 (MQ=255)
                                                                                                                                           aGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAg                                 <  2:4788141/128‑1 (MQ=255)
                                                                                                                                           aGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg           >  1:2462405/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt     >  1:3934868/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt     >  1:261273/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt     >  1:4928705/1‑150 (MQ=255)
                                                                                                                                                   cAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGc   >  1:882049/1‑150 (MQ=255)
                                                                                                                                                    aaaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGcc  >  1:1072188/1‑150 (MQ=255)
                                                                                                                                                    aaaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGcc  >  1:1071948/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCC  >  NC_002947/2901053‑2901350

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