Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,904,230 A→T 5.0% 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 5.0% 815.4 / 21.7 238intergenic (+50/+18)PP_2552/PP_2553DOPA decarboxylase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (134/92);  new base T (7/5);  total (141/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 258 total aligned reads displayed.

CGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCAG  >  NC_002947/2904081‑2904379
                                                                                                                                                      |                                                                                                                                                     
cGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaaa                                                                                                                                                       <  1:1035740/150‑1 (MQ=255)
cGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCaaa                                                                                                                                                       <  2:5653522/150‑1 (MQ=255)
   gcagTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagca                                                                                                                                                    <  2:4417938/147‑1 (MQ=255)
        ggATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGcccc                                                                                                                                               <  2:2459925/150‑1 (MQ=255)
           tGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGc                                                                                                                                            <  2:464134/150‑1 (MQ=255)
             gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCaa                                                                                                                                          <  2:6304630/150‑1 (MQ=255)
             gTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCGGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGcgc                                                                                                                                          <  1:6471246/150‑3 (MQ=255)
                   gTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATc                                                                                                                                    >  2:4476224/1‑150 (MQ=255)
                   gTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATc                                                                                                                                    >  1:6171532/1‑150 (MQ=255)
                      tCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGt                                                                                                                                 <  2:5643216/150‑1 (MQ=255)
                        gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCa                                                                                                                                   >  2:3650386/1‑146 (MQ=255)
                        gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                               >  2:1937762/1‑150 (MQ=255)
                        gATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTgg                                                                                                                               <  1:4100675/150‑1 (MQ=255)
                             gTGCGCTGCCGACCGCGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGGATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCACTCTATCAGTGGCTCAc                                                                                                                          >  1:3033326/1‑150 (MQ=255)
                                cgcTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt                                                                                                                                         <  1:5640271/132‑1 (MQ=255)
                                  cTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCt                                                                                                                     <  2:3361713/150‑1 (MQ=255)
                                       gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt                                                                                                                                         <  2:3932280/125‑1 (MQ=255)
                                           gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                                                                                                            <  1:2340212/150‑1 (MQ=255)
                                            aGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTGGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAcc                                                                                                           <  1:1680954/150‑1 (MQ=255)
                                                gggggATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTa                                                                                                       >  2:5050890/1‑150 (MQ=255)
                                                    gATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggc                                                                                                   >  2:5060726/1‑150 (MQ=255)
                                                          cAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCtt                                                                                             >  2:1619915/1‑150 (MQ=255)
                                                          cAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCtt                                                                                             <  1:4681902/150‑1 (MQ=255)
                                                           aGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTc                                                                                            <  1:5202645/150‑1 (MQ=255)
                                                            gCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCa                                                                                           >  1:4951510/1‑150 (MQ=255)
                                                             cGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCAt                                                                                          <  1:6489370/150‑1 (MQ=255)
                                                               gCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAg                                                                                        >  2:4297627/1‑150 (MQ=255)
                                                                 tgtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGTGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                      <  1:5512216/150‑1 (MQ=255)
                                                                 tgtgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                      >  2:2444427/1‑150 (MQ=255)
                                                                 cttgatcACGTCCTGCAGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc                                                                                      >  2:6278370/7‑150 (MQ=255)
                                                                   tgGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTg                                                                                    >  1:2701510/1‑150 (MQ=255)
                                                                     ggCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGcc                                                                                  >  2:2729492/1‑150 (MQ=255)
                                                                        aCGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCttt                                                                               >  1:5733140/1‑150 (MQ=255)
                                                                              tGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCAt                                                                          >  2:6552637/1‑150 (MQ=255)
                                                                              tGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCAt                                                                          >  1:5073129/1‑150 (MQ=255)
                                                                               gCAGGTCTTGATCAAGGGCTTAAGCCTGTACGGTACCTGTATGAGGGGGCTTGCCCGCGATGGGGCTGCTAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGATAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                         <  1:770789/150‑1 (MQ=255)
                                                                               gCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                         <  1:1566672/150‑1 (MQ=255)
                                                                               gCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                         <  1:2223158/150‑1 (MQ=255)
                                                                               gCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATc                                                                         <  1:3866716/150‑1 (MQ=255)
                                                                                  ggACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGc                                                                      <  1:1777919/150‑1 (MQ=255)
                                                                                   gACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCt                                                                     >  1:6731271/1‑150 (MQ=255)
                                                                                    aCGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCtt                                                                    <  1:6019891/150‑1 (MQ=255)
                                                                                     cGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTg                                                                   <  1:5341250/150‑1 (MQ=255)
                                                                                       tGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGc                                                                 >  2:1067230/1‑150 (MQ=255)
                                                                                       tGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGc                                                                 <  2:4022088/150‑1 (MQ=255)
                                                                                        gATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCg                                                                >  2:2327585/1‑150 (MQ=255)
                                                                                        gATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCg                                                                >  1:4400372/1‑150 (MQ=255)
                                                                                          tCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTc                                                              >  1:4759251/1‑150 (MQ=255)
                                                                                            aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCg                                                            >  1:6353974/1‑150 (MQ=255)
                                                                                            aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCg                                                            >  1:5602663/1‑150 (MQ=255)
                                                                                            aaGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCg                                                            <  2:839636/150‑1 (MQ=255)
                                                                                               ggCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGt                                                         >  1:4339095/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                      >  1:6850616/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                      >  2:3740121/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                      >  2:5432733/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                      >  2:6231587/1‑150 (MQ=255)
                                                                                                  tGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGc                                                      <  2:18652/150‑1 (MQ=255)
                                                                                                    aaGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGccc                                                    >  1:1085674/1‑150 (MQ=255)
                                                                                                    aaGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGccc                                                    >  1:2133595/1‑150 (MQ=255)
                                                                                                      gCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAg                                                  >  2:1685935/1‑150 (MQ=255)
                                                                                                       ccTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAgg                                                 >  1:305853/1‑150 (MQ=255)
                                                                                                        cTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGt                                                >  1:1476411/1‑150 (MQ=255)
                                                                                                        cTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGt                                                <  1:908524/150‑1 (MQ=255)
                                                                                                          gTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCt                                              <  1:2896288/150‑1 (MQ=255)
                                                                                                           tACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCtt                                             <  2:2117969/150‑1 (MQ=255)
                                                                                                                tACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATc                                        <  1:1956339/150‑1 (MQ=255)
                                                                                                                   cTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGcc                                     >  2:3211116/1‑150 (MQ=255)
                                                                                                                     gTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAg                                   <  1:5041324/150‑1 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                  >  1:6111983/1‑150 (MQ=255)
                                                                                                                      tATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAgg                                  >  1:4010157/1‑150 (MQ=255)
                                                                                                                       cttcgCGGGCAAGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCGTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGc                                 >  1:26003/5‑150 (MQ=255)
                                                                                                                        ttcgCGGGCAAGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCg                                >  2:5482881/4‑150 (MQ=255)
                                                                                                                          aGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGt                              <  2:1535889/150‑1 (MQ=255)
                                                                                                                          aGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGt                              >  1:3677046/1‑150 (MQ=255)
                                                                                                                           gCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTg                             <  2:6577145/150‑1 (MQ=255)
                                                                                                                                cTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGttt                        <  1:3349793/150‑1 (MQ=255)
                                                                                                                                  tGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTAGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTcc                      >  2:5900454/1‑150 (MQ=255)
                                                                                                                                   gCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCg                     >  2:314174/1‑150 (MQ=255)
                                                                                                                                     ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                   >  2:6751262/1‑150 (MQ=255)
                                                                                                                                     ccGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCggg                   <  2:6353974/150‑1 (MQ=255)
                                                                                                                                       gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCgggtc                 <  2:4622419/150‑3 (MQ=255)
                                                                                                                                       gcgAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGc                 >  2:2007731/1‑150 (MQ=255)
                                                                                                                                         gAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCaa               >  2:4208094/1‑150 (MQ=255)
                                                                                                                                         gAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCaa               <  2:6172793/150‑1 (MQ=255)
                                                                                                                                           aGGGGCTGCTTTGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg             >  1:6658478/1‑150 (MQ=255)
                                                                                                                                           aGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATg             <  1:5576942/150‑1 (MQ=255)
                                                                                                                                            ggggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGt            <  2:2065162/150‑1 (MQ=255)
                                                                                                                                            ggggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGt            >  1:2288091/1‑150 (MQ=255)
                                                                                                                                             gggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa           <  2:6775425/150‑1 (MQ=255)
                                                                                                                                             gggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTa           <  2:2935967/150‑1 (MQ=255)
                                                                                                                                              ggCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAg          <  1:2313255/150‑1 (MQ=255)
                                                                                                                                               gCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAACTTTACGGCTGGCTTCATAGGCTGACTTTTCCATCGGCTTGGAGTCAGGGTTGCCCAGGTCTTCAATCGCCAGGCGGTGGGTTACCGGGGCAATGTAgg         >  1:2406826/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt       >  2:5971268/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt       >  2:4107067/1‑150 (MQ=255)
                                                                                                                                                 tGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCt       >  2:5565403/1‑150 (MQ=255)
                                                                                                                                                  gCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTg      >  1:3944674/1‑150 (MQ=255)
                                                                                                                                                    aaaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGcc    >  2:2857089/1‑150 (MQ=255)
                                                                                                                                                     aaGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCa   >  1:3801824/1‑150 (MQ=255)
                                                                                                                                                      agcagcCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCAg  <  2:5733140/150‑1 (MQ=255)
                                                                                                                                                      agcagcCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCAg  >  1:3588007/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGCTTGGCGTCCGGGTTGCCCAGGTCTTCCATCGCCAGGCGGTGGGTTTCCGGGGCAATGTAGGCTGCCAG  >  NC_002947/2904081‑2904379

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: