Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,892 G→A 5.3% 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,716,8920GA 5.3% 960.6 / 26.0 300A27A (GCC→GCTpaaBring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (148/135);  new base A (3/13);  total (151/149)
Fisher's exact test for biased strand distribution p-value = 9.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 316 total aligned reads displayed.

ACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACT  >  NC_002947/3716743‑3717040
                                                                                                                                                     |                                                                                                                                                    
 ccTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATcggcgg                                                                                                                                                     <  1:1360416/150‑1 (MQ=255)
  cTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCggcggc                                                                                                                                                    <  1:6187980/150‑1 (MQ=255)
    tGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGt                                                                                                                                                  <  1:6331570/150‑1 (MQ=255)
     gtcgtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTg                                                                                                                                                 >  1:5808826/1‑150 (MQ=255)
     gtcgtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTg                                                                                                                                                 >  1:531194/1‑150 (MQ=255)
      tcgtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGc                                                                                                                                                <  1:5721678/150‑1 (MQ=255)
        gtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCAc                                                                                                                                              >  1:3916110/1‑150 (MQ=255)
        gtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCAc                                                                                                                                              >  2:2359659/1‑150 (MQ=255)
        gtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCAc                                                                                                                                              >  1:1420268/1‑150 (MQ=255)
                gggCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGa                                                                                                                                      >  2:5311278/1‑150 (MQ=255)
                    gAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATg                                                                                                                                  >  1:254839/1‑150 (MQ=255)
                      aCAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGtt                                                                                                                                >  2:6562973/1‑150 (MQ=255)
                       cAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGttt                                                                                                                               <  2:6839173/150‑1 (MQ=255)
                        aGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtg                                                                                                                              >  2:1803107/1‑150 (MQ=255)
                         gCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgt                                                                                                                             >  1:3689531/1‑150 (MQ=255)
                         gCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgt                                                                                                                             >  1:2740400/1‑150 (MQ=255)
                         gCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgt                                                                                                                             <  1:5181356/150‑1 (MQ=255)
                           gggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCAGGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                           <  2:2068013/150‑1 (MQ=255)
                            ggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGt                                                                                                                          >  1:2954327/1‑150 (MQ=255)
                              tCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTa                                                                                                                        >  1:5137325/1‑150 (MQ=255)
                                 ttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAgg                                                                                                                     <  2:2993759/150‑1 (MQ=255)
                                  tctcGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGc                                                                                                                    >  2:6444824/1‑150 (MQ=255)
                                  tctcGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGc                                                                                                                    >  2:905380/1‑150 (MQ=255)
                                    tcgtcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCg                                                                                                                  >  1:5913836/1‑150 (MQ=255)
                                      gtcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGtg                                                                                                                >  1:908929/1‑150 (MQ=255)
                                       tcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGtgc                                                                                                               <  1:1974377/150‑1 (MQ=255)
                                            agagGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGtgcttgct                                                                                                          >  2:5160158/1‑150 (MQ=255)
                                            agagGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGGTTGTGGTTAAGGCCGtgcttgct                                                                                                          >  1:3794987/1‑150 (MQ=255)
                                                 aggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCa                                                                                                     <  1:6584389/150‑1 (MQ=255)
                                                  ggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCAc                                                                                                    >  1:6740544/1‑150 (MQ=255)
                                                   gCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACg                                                                                                   >  2:6300315/1‑150 (MQ=255)
                                                   gCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCGGCCGACCGGTTTGGGGGTCGGGCCGTGCTTGCTGCGCACg                                                                                                   >  2:857677/1‑150 (MQ=255)
                                                    cGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGa                                                                                                  <  2:5887156/150‑1 (MQ=255)
                                                          tCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGGCCGCTGCCGTGGGGGTGGTGGGGCCGGCCGTGCTTGGTGcgcgcggaggggg                                                                                            >  2:4438862/1‑143 (MQ=255)
                                                           cAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTc                                                                                           <  2:4202754/150‑1 (MQ=255)
                                                           cAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTc                                                                                           <  1:4146932/150‑1 (MQ=255)
                                                                ccGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGcc                                                                                                                                        >  2:1710747/1‑100 (MQ=255)
                                                                ccGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                      >  1:2374676/1‑150 (MQ=255)
                                                                ccGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                      >  1:1050630/1‑150 (MQ=255)
                                                                ccGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTaga                                                                                      >  1:1324387/1‑150 (MQ=255)
                                                                  gAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAgg                                                                                    <  2:4281465/150‑1 (MQ=255)
                                                                    aGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGt                                                                                  <  2:391741/150‑1 (MQ=255)
                                                                     ggCACCACCCACAGACTCACACCGTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTc                                                                                 <  1:6119440/150‑1 (MQ=255)
                                                                       caccacCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCca                                                                               <  2:1529149/150‑1 (MQ=255)
                                                                        accaccCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGtgcttg                                                                                                            <  2:3260780/120‑1 (MQ=255)
                                                                        accaccCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGGTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAGCACGTCGTAGAGGGTCcag                                                                              >  2:5096227/1‑150 (MQ=255)
                                                                         ccacccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcaga                                                                             <  2:4696018/150‑1 (MQ=255)
                                                                         ccacccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcaga                                                                             <  1:2867028/150‑1 (MQ=255)
                                                                          cacccacAGACTCACACCTTCGCTACGGCGGGTGTAGATCTCGCGGGCGTTCTCGATGGCCATGGCGGCTTCGGCGGCGTGCACGCTGCCGCCAGGGGTGGGGGTATGGCCGTGCGTGCTGCGCTCGAACACTTCGGGGAGGGTCcagac                                                                            >  2:1441570/1‑150 (MQ=255)
                                                                             ccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACCTGTTGGTGGTTCAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagacaga                                                                         >  2:5507428/1‑150 (MQ=255)
                                                                             ccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcagacaga                                                                         >  2:4782926/1‑150 (MQ=255)
                                                                              cacaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCagacagac                                                                        <  1:4921106/150‑1 (MQ=255)
                                                                              cacaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGGTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCagacagac                                                                        >  1:429997/1‑150 (MQ=255)
                                                                               acaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAgacagaca                                                                       >  1:5266853/1‑150 (MQ=255)
                                                                               acaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAgacagaca                                                                       >  1:6051812/1‑150 (MQ=255)
                                                                                  gACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACAttt                                                                    >  1:2413260/1‑150 (MQ=255)
                                                                                   aCTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACAtttt                                                                   <  1:3969069/150‑1 (MQ=255)
                                                                                     tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGTCATTTTCg                                                                 >  2:463387/1‑150 (MQ=255)
                                                                                     tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCg                                                                 >  1:1228321/1‑150 (MQ=255)
                                                                                     tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCAcgccgc                                                                                                                                         <  2:3054066/78‑1 (MQ=255)
                                                                                       acacCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACAGGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTGGAGGGTCCGGACAGACATTTTCGAt                                                               >  2:5002068/1‑150 (MQ=255)
                                                                                         acCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTACTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAtct                                                             >  1:211573/1‑150 (MQ=255)
                                                                                          ccTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAtctc                                                            >  1:5051016/1‑150 (MQ=255)
                                                                                          ccTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCCCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGCCATGTTGGTGGTTAAGGCCGTGCTTGCTGCGCACGAAGACTTCGTAGAGGGTCCAGACAGACAGTTTCGAtctc                                                            >  2:722385/1‑150 (MQ=255)
                                                                                           cTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTcc                                                           <  1:5523327/150‑1 (MQ=255)
                                                                                               gCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCAt                                                       >  1:3572908/1‑150 (MQ=255)
                                                                                                cTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGAGCGGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACGGCGTAGAGGGTCCAGACAGACAGTTTCGATCTCCGCGTc                                                      >  2:5997812/1‑150 (MQ=255)
                                                                                                      cggGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATcgc                                                >  1:2615587/1‑150 (MQ=255)
                                                                                                             agagCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAg                                         >  1:3092042/1‑150 (MQ=255)
                                                                                                              gagCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAgg                                        <  2:6070661/150‑1 (MQ=255)
                                                                                                              gagCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAgg                                        >  2:1212378/1‑150 (MQ=255)
                                                                                                                 cTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcgg                                     <  2:2635629/150‑1 (MQ=255)
                                                                                                                 cTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcgg                                     >  1:5946819/1‑150 (MQ=255)
                                                                                                                  tCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcggc                                    <  1:4393116/150‑1 (MQ=255)
                                                                                                                      gggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAggg                                                                                   <  2:4070188/99‑1 (MQ=255)
                                                                                                                      gggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTc                                <  2:1356786/150‑1 (MQ=255)
                                                                                                                      gggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTc                                <  2:2785479/150‑1 (MQ=255)
                                                                                                                       ggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCt                               <  1:3381769/150‑1 (MQ=255)
                                                                                                                        gCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCtt                              <  1:6359180/150‑1 (MQ=255)
                                                                                                                        gCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCtt                              <  1:5500524/150‑1 (MQ=255)
                                                                                                                         cGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTg                             <  2:5718780/150‑1 (MQ=255)
                                                                                                                           ttCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGtt                           <  1:3623140/150‑1 (MQ=255)
                                                                                                                            tctcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTc                          >  1:6276548/1‑150 (MQ=255)
                                                                                                                            tctcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTc                          <  2:3510938/150‑1 (MQ=255)
                                                                                                                            tctcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTAAGGCGGCGTTCTTGTTc                          <  1:2174657/150‑1 (MQ=255)
                                                                                                                            tctcGAAGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGCGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTc                          >  2:4950261/1‑150 (MQ=255)
                                                                                                                             ctcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCt                         >  1:719407/1‑150 (MQ=255)
                                                                                                                             ctcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCt                         <  2:3607087/150‑1 (MQ=255)
                                                                                                                                gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACCTGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGtt                      >  1:3020152/1‑150 (MQ=255)
                                                                                                                                   ggCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgct                   >  1:447395/1‑150 (MQ=255)
                                                                                                                                     ccATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgcttg                 >  2:4505758/1‑150 (MQ=255)
                                                                                                                                       aTGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTgcg               >  2:2889505/1‑150 (MQ=255)
                                                                                                                                        tGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTgcgc              <  2:3068928/150‑1 (MQ=255)
                                                                                                                                           cggcATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCa           <  1:1755422/150‑1 (MQ=255)
                                                                                                                                           cggcATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCa           <  2:3875899/150‑1 (MQ=255)
                                                                                                                                            ggcATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCAt          >  1:6777400/1‑150 (MQ=255)
                                                                                                                                               aTCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAgg       >  1:4617475/1‑150 (MQ=255)
                                                                                                                                                   gcggcgTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCAc   >  2:1400223/1‑150 (MQ=255)
                                                                                                                                                   gcggcgTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCAc   >  1:2285733/1‑150 (MQ=255)
                                                                                                                                                    cggcgTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACt  <  1:1593926/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACT  >  NC_002947/3716743‑3717040

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