Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,899,237 A→T 8.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,899,2370AT 8.1% 180.9 / 5.9 62intergenic (+49/+19)PP_2552/PP_2553DOPA decarboxylase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (3/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCT  >  NC_002947/2899095‑2899382
                                                                                                                                               |                                                                                                                                                 
gTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagcagc                                                                                                                                             <  3:209938/149‑1 (MQ=255)
 tGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGcc                                                                                                                                            >  7:78169/1‑149 (MQ=255)
  ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGccc                                                                                                                                           >  4:290772/1‑149 (MQ=255)
       gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTTCAGCGGCTGTGGGCACGTCTGCATGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCa                                                                                                                                      <  2:42503/149‑1 (MQ=255)
       gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCa                                                                                                                                      <  8:110148/149‑1 (MQ=255)
        tGCGGGTTTCGATTGGTGCGCTGCTGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCaa                                                                                                                                     <  1:193627/149‑1 (MQ=255)
          cGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGCAGGGGCTGCTTTGCAGCCCCttcgcgg                                                                                                                                   <  5:12634/149‑7 (MQ=255)
          cGGGTGTCGATTGGTGCCCGGCCGCCCGAGCGGGGGGATGGGCAGCGGCTGTGGCCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATc                                                                                                                                   <  6:253606/149‑1 (MQ=255)
             gTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTAt                                                                                                                                >  4:56344/1‑149 (MQ=255)
              tttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATc                                                                                                                               <  8:78169/149‑1 (MQ=255)
               ttCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTAt                                                                                                                                >  6:210965/1‑147 (MQ=255)
                  gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagc                                                                                                                                                >  3:105828/1‑128 (MQ=255)
                  gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagc                                                                                                                                                <  4:105828/128‑1 (MQ=255)
                  gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGt                                                                                                                            <  5:210965/148‑1 (MQ=255)
                   aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                          <  3:142209/149‑1 (MQ=255)
                          cgcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCg                                                                                                                   >  7:169399/1‑149 (MQ=255)
                          cgcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCg                                                                                                                   >  5:358082/1‑149 (MQ=255)
                                 gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCttcgcg                                                                                                                                    <  4:245293/125‑6 (MQ=39)
                                 gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCttcgcg                                                                                                                                    >  3:245293/1‑120 (MQ=39)
                                 gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa                                                                                                            <  3:182465/149‑1 (MQ=255)
                                 gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGa                                                                                                              >  3:246540/1‑147 (MQ=255)
                                 gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCAGCGAAGGGGCTGCAAAGCAGCCACTGCAATCAATCAGTGGCTCACCGCCTGa                                                                                                              >  2:233968/1‑147 (MQ=255)
                                    cgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGc                                                                                                          >  5:19090/1‑148 (MQ=255)
                                     gaGCGGGGGTATGTGCAGCGGCTTTGGGCTCGTCTGCAGGACGTTATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTa                                                                                                        <  6:19090/149‑1 (MQ=255)
                                            gggATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTAc                                                                                                 <  4:246540/149‑1 (MQ=255)
                                            gggATGGGCAGCGGCTGTGGGCACGTCTGCAGGACGGGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTAc                                                                                                 <  6:283025/149‑1 (MQ=255)
                                               aTGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggc                                                                                              >  3:59233/1‑149 (MQ=255)
                                                   gCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc                                                                                          >  1:63879/1‑149 (MQ=255)
                                                       cGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATc                                                                                                                               >  8:248871/1‑108 (MQ=255)
                                                       cGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATc                                                                                                                               <  7:248871/108‑1 (MQ=255)
                                                         gCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATa                                                                                    <  4:316782/149‑1 (MQ=255)
                                                         gCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATa                                                                                    <  4:107146/149‑1 (MQ=255)
                                                              gatcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTg                                                                               >  4:305314/4‑149 (MQ=255)
                                                                   cGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCttt                                                                          <  3:56344/149‑1 (MQ=255)
                                                                    gTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCtttt                                                                         >  3:306745/1‑149 (MQ=255)
                                                                     tCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTGCCTTTACGGCTGGCTTCATAGGCTGCCTTTTc                                                                        <  6:358082/149‑1 (MQ=255)
                                                                     tCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTc                                                                        <  8:169399/149‑1 (MQ=255)
                                                                          cAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                    >  5:20111/1‑149 (MQ=255)
                                                                          cAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                    <  4:306745/149‑1 (MQ=255)
                                                                             gACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGc                                                                 >  7:306695/1‑149 (MQ=255)
                                                                                  gATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATCGGCTGCCTTTTCCATCGGCTTGGc                                                            >  6:259299/1‑149 (MQ=255)
                                                                                        gggCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCgg                                                      >  5:237506/1‑149 (MQ=255)
                                                                                              aaGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGcc                                                <  4:59233/149‑1 (MQ=255)
                                                                                               aGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGccc                                               <  5:259299/149‑1 (MQ=255)
                                                                                                  cTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCtt                                                                                                    >  7:103077/1‑93 (MQ=255)
                                                                                                  cTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCtt                                                                                                    <  8:103077/93‑1 (MQ=255)
                                                                                                          tACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCAt                                    <  3:162659/149‑1 (MQ=255)
                                                                                                            ccTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATc                                   <  8:306695/148‑1 (MQ=255)
                                                                                                               gTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCg                                                                   <  6:234941/113‑1 (MQ=255)
                                                                                                               gTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCg                                                                   >  5:234941/1‑113 (MQ=255)
                                                                                                                  tGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                             <  3:290772/148‑1 (MQ=255)
                                                                                                                    aGCGGGCTTGCGCGCGAAGGGTCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCgg                          <  1:233968/149‑1 (MQ=255)
                                                                                                                         gCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGt                     >  8:282333/1‑149 (MQ=255)
                                                                                                                              cccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCg                >  5:337252/1‑149 (MQ=255)
                                                                                                                                 gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCgggg             <  6:20111/149‑1 (MQ=255)
                                                                                                                                     aGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACGTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAAt         >  4:63675/1‑149 (MQ=255)
                                                                                                                                       gggCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg        >  7:52790/1‑148 (MQ=255)
                                                                                                                                       gggCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg        >  7:284186/1‑148 (MQ=255)
                                                                                                                                        ggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa      <  3:305314/149‑1 (MQ=255)
                                                                                                                                          cTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAgg    >  5:248671/1‑149 (MQ=255)
                                                                                                                                           tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGc   >  1:291197/1‑149 (MQ=255)
                                                                                                                                            gCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt  >  8:128401/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCT  >  NC_002947/2899095‑2899382

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