Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,711,943 G→A 9.4% A27A (GCC→GCT paaB ← ring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,711,9430GA 9.4% 258.4 / 15.2 96A27A (GCC→GCTpaaBring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base A (4/5);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

CTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCA  >  NC_002947/3711796‑3712089
                                                                                                                                                   |                                                                                                                                                  
cTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCggcggc                                                                                                                                                  >  2:1210298/1‑150 (MQ=255)
 ttGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGgcggcg                                                                                                                                                 >  1:2207628/1‑150 (MQ=255)
    tcgtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGc                                                                                                                                              >  1:920212/1‑150 (MQ=255)
           aGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGc                                                                                                                                       >  1:1448143/1‑150 (MQ=255)
              gggCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGa                                                                                                                                    <  1:844866/150‑1 (MQ=255)
              gggCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGa                                                                                                                                    <  1:545799/150‑1 (MQ=255)
                gCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACa                                                                                                                                  >  1:1026311/1‑150 (MQ=255)
                          ggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGt                                                                                                                        <  2:1761388/150‑1 (MQ=255)
                          ggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGt                                                                                                                        <  2:330171/150‑1 (MQ=255)
                               ttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                         <  2:2163518/144‑1 (MQ=255)
                               ttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                         >  1:2163518/1‑144 (MQ=255)
                               ttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACGCACACCTTCGCTACGGCGGGTGGAGAGCGCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAgg                                                                                                                   <  2:971860/150‑1 (MQ=255)
                                tctcGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGat                                                                                                                                   >  2:1591076/1‑132 (MQ=255)
                                tctcGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGat                                                                                                                                   <  1:1591076/133‑2 (MQ=255)
                                 ctcGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGcc                                                                                                                 <  1:863605/150‑1 (MQ=255)
                                   cgtcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGt                                                                                                               <  2:1026311/150‑1 (MQ=255)
                                    gtcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGtg                                                                                                              >  2:1434658/1‑150 (MQ=255)
                                       gggAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGttt                                                                                                                             >  2:1803002/1‑132 (MQ=255)
                                       gggAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGttt                                                                                                                             <  1:1803002/132‑1 (MQ=255)
                                           gaggagGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTtgctg                                                                                                       >  2:1315606/1‑150 (MQ=255)
                                            aggaggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTtgctgc                                                                                                      >  2:1491907/1‑150 (MQ=255)
                                             ggaggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTgcg                                                                                                     >  1:287037/1‑150 (MQ=255)
                                                ggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCAc                                                                                                  <  1:52889/150‑1 (MQ=255)
                                                   ggTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGaa                                                                                               <  1:1315606/150‑1 (MQ=255)
                                                       aTCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTg                                                                                                                                               >  2:1322632/1‑98 (MQ=255)
                                                       aTCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTg                                                                                                                                               <  1:1322632/98‑1 (MQ=255)
                                                       aTCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGAGGGCGTTCTCGAGTGGCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACt                                                                                           >  1:1498404/1‑150 (MQ=255)
                                                        tCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTCCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACtt                                                                                          >  2:2124262/1‑150 (MQ=255)
                                                            cgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGATGCCGACATGTTTGTGCTTAAGGCCGTGCTTGCTCCGCACGAACGCTTCGTa                                                                                      <  2:1900097/150‑1 (MQ=255)
                                                             gcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTag                                                                                     >  1:1181378/1‑150 (MQ=255)
                                                              ccGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                    >  1:1780776/1‑150 (MQ=255)
                                                              ccGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                    <  1:1491907/150‑1 (MQ=255)
                                                                gAAGGCACCCCCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCGCGCGGGCGTTCTCGATGGCCATGGCGGCAGCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAgg                                                                                  <  2:2102365/150‑1 (MQ=255)
                                                                   ggCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTc                                                                               <  1:1210298/150‑1 (MQ=255)
                                                                     caccacCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgt                                                                                                                           >  1:795043/1‑104 (MQ=255)
                                                                     caccacCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgt                                                                                                                           <  2:795043/104‑1 (MQ=255)
                                                                      accaccCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcag                                                                            >  1:616194/1‑150 (MQ=255)
                                                                        cacccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTTCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagac                                                                          >  2:534985/1‑150 (MQ=255)
                                                                        cacccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagac                                                                          >  1:1193942/1‑150 (MQ=255)
                                                                        cacccacAGACTCACACCTTCGCTAAGGAGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagac                                                                          >  1:677279/1‑150 (MQ=255)
                                                                           ccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagacaga                                                                       >  2:1818343/1‑150 (MQ=255)
                                                                           ccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATAGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagacaga                                                                       >  2:1406544/1‑150 (MQ=255)
                                                                            cacaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCTGGCGTTCTCGATGGCCATGGTGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCagacagac                                                                      >  2:1393797/1‑150 (MQ=255)
                                                                             acaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACa                                                                                                                                  <  2:1011897/89‑1 (MQ=255)
                                                                             acaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACa                                                                                                                                  >  1:1011897/1‑89 (MQ=255)
                                                                                  cTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCTGGCGTTCTCGATGGCCATGGTGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                <  1:1393797/150‑1 (MQ=255)
                                                                                  cTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                <  2:200658/150‑1 (MQ=255)
                                                                                   tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAac                                                                                              >  2:995838/1‑119 (MQ=255)
                                                                                   tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAac                                                                                              <  1:995838/119‑1 (MQ=255)
                                                                                   tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCg                                                               >  1:850003/1‑150 (MQ=255)
                                                                                     acacCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAt                                                             <  2:1498404/150‑1 (MQ=255)
                                                                                     acacCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCAc                                                                                                                                            <  2:2059849/71‑1 (MQ=255)
                                                                                     acacCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCAc                                                                                                                                            >  1:2059849/1‑71 (MQ=255)
                                                                                      cacCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAtc                                                            >  1:2272533/1‑150 (MQ=255)
                                                                                       acCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAtct                                                           >  1:1516019/1‑150 (MQ=255)
                                                                                         cTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                         >  1:1578405/1‑86 (MQ=255)
                                                                                         cTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                         <  2:1578405/86‑1 (MQ=255)
                                                                                         cTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTcc                                                         <  2:2207628/150‑1 (MQ=255)
                                                                                          ttCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGGCCAGACAGCCATTTTCGATCTCCg                                                        >  1:1579560/1‑150 (MQ=255)
                                                                                            cGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCa                                                      <  1:1104563/150‑1 (MQ=255)
                                                                                              cTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGc                                                                                                                                              >  1:1689842/1‑60 (MQ=255)
                                                                                              cTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGc                                                                                                                                              <  2:1689842/60‑1 (MQ=255)
                                                                                               tACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCt                                                   <  1:260030/150‑1 (MQ=255)
                                                                                                     gggTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGa                                                                                                                                    <  1:2046327/63‑1 (MQ=38)
                                                                                                     gggTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGa                                                                                                                                    >  2:2046327/1‑63 (MQ=38)
                                                                                                           agagCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCNCGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAg                                       >  1:863997/1‑150 (MQ=255)
                                                                                                           agagCTCCGTGGCGTTCTCGGTGGTCATGGCGGCATCGGCGGTGTGCACGTTGGCTACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAg                                       <  2:863997/150‑1 (MQ=255)
                                                                                                            gagCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAgg                                      <  1:1434658/150‑1 (MQ=255)
                                                                                                             agCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggc                                     <  1:871495/150‑1 (MQ=255)
                                                                                                                tCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTgcgc                                                                                                    <  1:2243825/84‑1 (MQ=255)
                                                                                                                tCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTgcgc                                                                                                    >  2:2243825/1‑84 (MQ=255)
                                                                                                                tCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcggc                                  >  2:384410/1‑150 (MQ=255)
                                                                                                                tCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcggc                                  <  2:616194/150‑1 (MQ=255)
                                                                                                                  gcgGGCGTTCTCGATGGCCATGGTGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                >  2:1419054/1‑118 (MQ=255)
                                                                                                                  gcgGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                <  1:1419054/118‑1 (MQ=255)
                                                                                                                   cgGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGtt                               <  2:1780776/150‑1 (MQ=255)
                                                                                                                    gggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                          >  2:1483388/1‑58 (MQ=255)
                                                                                                                    gggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                          <  1:1483388/58‑1 (MQ=255)
                                                                                                                           ctcGATGGCCATGGCGGTATCGGCGGCGTGCACGCTGCTGACATGTTTGTGGTTAAGGCCGTGCTTGCTTCGCACGAGCACTTCTTAGTGGGTCCAGACGGACATTTTCGCTCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCt                       <  2:152605/150‑1 (MQ=255)
                                                                                                                            tcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcggc                                  >  1:1404663/1‑138 (MQ=255)
                                                                                                                            tcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcggc                                  <  2:1404663/138‑1 (MQ=255)
                                                                                                                            tcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTg                      >  1:934244/1‑150 (MQ=255)
                                                                                                                              gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGACCAAATCGTATAGGGTCCAGCAAAACACATTAGAACTCGTCTGCTATCGCCTCTCAGGCGCCGTTCTGGTTTTGtt                    >  1:540754/1‑150 (MQ=255)
                                                                                                                              gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                    >  1:338801/1‑86 (MQ=255)
                                                                                                                              gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGtt                    <  1:1818343/150‑1 (MQ=255)
                                                                                                                              gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGtt                    <  1:2124262/150‑1 (MQ=255)
                                                                                                                              gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCCACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                    <  2:338801/86‑1 (MQ=255)
                                                                                                                               aTGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCt                       >  2:2387953/1‑146 (MQ=255)
                                                                                                                               aTGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCt                       <  1:2387953/146‑1 (MQ=255)
                                                                                                                                tGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgc                  <  2:1181378/150‑1 (MQ=255)
                                                                                                                                    cATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgcttgc              >  2:94386/1‑150 (MQ=255)
                                                                                                                                     aTGGCGGCATCGGCGGCGTGCACGCTGCCGACAGGTTTGTGGTTAAGGCCGTGCTTGCTGCGCCCGAACACTTCGTAGAGGGTCCAGACAGAAATTCTCGATAACGTCATCTAGCGCCAGCCAGTCCCCTTCCTTGTTCGGCTGCTTgcg             >  2:1353550/1‑150 (MQ=255)
                                                                                                                                            cATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTTCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAg      <  1:534985/150‑1 (MQ=255)
                                                                                                                                             aTCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAgg     >  1:567938/1‑150 (MQ=255)
                                                                                                                                              tCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAgg     >  1:991888/1‑149 (MQ=255)
                                                                                                                                              tCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAgg     <  2:991888/149‑1 (MQ=255)
                                                                                                                                              tCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGc    <  2:934244/150‑1 (MQ=255)
                                                                                                                                              tCGGCAGCGTGCACGCCGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGc    >  2:1869759/1‑150 (MQ=255)
                                                                                                                                                ggcggcGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCa  >  1:1235120/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCA  >  NC_002947/3711796‑3712089

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