Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,711,943 G→A 6.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,711,9430GA 6.3% 228.8 / 4.7 79A27A (GCC→GCTpaaBring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base A (3/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

ACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACT  >  NC_002947/3711794‑3712091
                                                                                                                                                     |                                                                                                                                                    
aCCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATcggcg                                                                                                                                                      >  1:272772/1‑150 (MQ=255)
 ccctGTCGTCCGAGGGGGCGAACAGCGCGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGAAGAGCTCGCTGGCGTTCTCGATGGCCATCGCAACATcagcga                                                                                                                                                     <  2:559342/147‑2 (MQ=255)
  cTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCggcggc                                                                                                                                                    >  2:146431/1‑150 (MQ=255)
  cTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCggcggc                                                                                                                                                    >  2:136787/1‑150 (MQ=255)
     gtcgtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATcggcg                                                                                                                                                      <  1:1257629/145‑1 (MQ=255)
     gtcgtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATcggcg                                                                                                                                                      >  2:1257629/1‑145 (MQ=255)
        gtcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCAc                                                                                                                                              >  1:1288288/1‑150 (MQ=255)
         tcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGGGGGCGTCCGCGACAGCCATGGCGGCATCGGCGGCGTGCACg                                                                                                                                             >  2:1331382/1‑150 (MQ=255)
         tcCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGAGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACg                                                                                                                                             >  1:1814773/1‑150 (MQ=255)
          ccGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGc                                                                                                                                            >  1:1166801/1‑150 (MQ=255)
            gAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCCCAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTg                                                                                                                                          >  1:991568/1‑150 (MQ=255)
                gggCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCCCAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGa                                                                                                                                      <  2:991568/150‑1 (MQ=255)
                       cAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGttt                                                                                                                               <  1:146431/150‑1 (MQ=255)
                       cAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGttt                                                                                                                               <  2:822979/150‑1 (MQ=255)
                       cAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGttt                                                                                                                               <  1:384678/150‑1 (MQ=255)
                         gCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgt                                                                                                                             >  1:500822/1‑150 (MQ=255)
                          cGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGCCTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                            <  1:1071565/150‑1 (MQ=255)
                          cGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                            >  1:568414/1‑150 (MQ=255)
                          cGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                            <  1:1657429/150‑1 (MQ=255)
                          cGGGTCTTTCTCGTCGGGAGAAGAGGAGGTGCACAGCGTCATATGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                            <  1:1066496/150‑1 (MQ=255)
                           gggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                           <  1:485918/150‑1 (MQ=255)
                           gggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                           >  1:1295451/1‑150 (MQ=255)
                           gggTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTgg                                                                                                                           <  2:1071849/150‑1 (MQ=255)
                             gTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGtt                                                                                                                         >  2:805854/1‑150 (MQ=255)
                               ttttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTaa                                                                                                                       <  2:1039970/149‑1 (MQ=255)
                               cTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTaa                                                                                                                       <  2:1295451/150‑1 (MQ=255)
                               cTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGAACGCTGCCGACATGTTTGTGGTTaa                                                                                                                       <  1:1698059/150‑1 (MQ=255)
                                tttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACgccgcc                                                                                                                                        >  1:1877766/1‑132 (MQ=255)
                                tttCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACgccgcc                                                                                                                                        <  2:1877766/132‑1 (MQ=255)
                                   cgcgtcgGGAGAGGGGGCGGTGATCAGCGCGGAAGGCACCACCCACAGTCTCACACCTTCGCTACGGCGGGTGAAGAGCTCGCGGGCGTTCTCGCTGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGcc                                                                                                                   <  2:1021488/148‑1 (MQ=255)
                                    tcgtcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCg                                                                                                                  <  2:163015/150‑1 (MQ=255)
                                                gaggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTgcgc                                                                                                      >  1:786290/1‑150 (MQ=255)
                                                gaggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGat                                                                                                                                     >  2:1853615/1‑118 (MQ=255)
                                                gaggCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCAGCGTGCACGCCGCCGat                                                                                                                                     <  1:1853615/119‑2 (MQ=255)
                                                   gCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACg                                                                                                   >  1:42004/1‑150 (MQ=255)
                                                    cGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGCTAAGGCCGTGCTTGCTGCGCACGa                                                                                                  <  2:1288288/150‑1 (MQ=255)
                                                     ggTGCTCAGCGCCGAGGGCACCCCCCAGAGTTTCACACCTTGGCGACGGCGGGTGTAGAGCTCGGGGGCGGTCTCGATGGCCATGGCGGCATCGGCGGCGGGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGaa                                                                                                 <  2:172453/150‑1 (MQ=255)
                                                      gTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAac                                                                                                >  2:1583012/1‑150 (MQ=255)
                                                          tCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACtt                                                                                            <  1:143952/150‑1 (MQ=255)
                                                           cAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTc                                                                                           <  2:992800/150‑1 (MQ=255)
                                                            aGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCg                                                                                          >  1:1742560/1‑150 (MQ=255)
                                                              cgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGt                                                                                                                                 <  1:1659884/109‑1 (MQ=255)
                                                              cgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGt                                                                                                                                 >  2:1659884/1‑109 (MQ=255)
                                                              cgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTa                                                                                        <  1:804631/150‑1 (MQ=255)
                                                              cgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTa                                                                                        <  2:1597997/150‑1 (MQ=255)
                                                              cgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTa                                                                                        <  1:805854/150‑1 (MQ=255)
                                                                  gAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAgg                                                                                    <  1:1488144/150‑1 (MQ=255)
                                                                   aaGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGATCTCGGTGGTCGCCGCATCATTAAAAACAACATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAggg                                                                                   >  2:1782/1‑150 (MQ=255)
                                                                      gCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGGTGGCCGCCTTATCATTAAAAAAAAAATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTcc                                                                                >  2:648531/1‑150 (MQ=255)
                                                                      gCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTcc                                                                                >  1:1642704/1‑150 (MQ=255)
                                                                      gCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGCTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGa                                                                                                                                      >  1:389541/1‑96 (MQ=255)
                                                                      gCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGCTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGa                                                                                                                                      <  2:389541/96‑1 (MQ=255)
                                                                       caccacCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCca                                                                               <  1:1583012/150‑1 (MQ=255)
                                                                       caccacCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCca                                                                               <  2:42004/150‑1 (MQ=255)
                                                                         ccacccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCcaga                                                                             >  1:315080/1‑150 (MQ=255)
                                                                          cacccacAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCccgac                                                                            <  2:1185310/150‑1 (MQ=255)
                                                                              cacaGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCagacagac                                                                        >  1:1783257/1‑150 (MQ=255)
                                                                                    cTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                  <  2:1330632/150‑1 (MQ=255)
                                                                                     tCACACCTTCGCTACGGCGGGTGTAGAGCTCGTGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                            >  1:1607867/1‑91 (MQ=255)
                                                                                     tCACACCTTCGCTACGGCGGGTGTAGAGCTCGTGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTtgtg                                                                                                                            <  2:1607867/91‑1 (MQ=255)
                                                                                     tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCg                                                                 >  1:1901256/1‑150 (MQ=255)
                                                                                       acacCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAt                                                               >  1:318441/1‑150 (MQ=255)
                                                                                          ccTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAtctc                                                            >  1:532885/1‑150 (MQ=255)
                                                                                           cctCGCTACGGCGGGTGTAGAGCTCGCGTGCGTTCTCGATGGCCCTGGCGGCATCGGCGGCGTGCACGCTGCCGACAGGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACCTTTTCGATCTcc                                                           <  2:170955/148‑1 (MQ=255)
                                                                                                           gtAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTc                                           >  2:426321/1‑150 (MQ=255)
                                                                                                                     cgGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGtt                                 <  1:940023/150‑1 (MQ=255)
                                                                                                                         cGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTg                             <  1:897915/150‑1 (MQ=255)
                                                                                                                           ttCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGtt                           <  2:1516328/150‑1 (MQ=255)
                                                                                                                              tcGATGGCCATGGCGGCATCGGCGGCGTGCGCGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTg                        >  1:1801246/1‑150 (MQ=255)
                                                                                                                               ggATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGt                       <  2:1864231/149‑1 (MQ=255)
                                                                                                                               cGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGt                       <  2:1742560/150‑1 (MQ=255)
                                                                                                                                gATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGATAGGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACCGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGtt                      >  1:346896/1‑150 (MQ=255)
                                                                                                                                  tGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACACGTTCGTAGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgc                    >  1:37585/1‑150 (MQ=255)
                                                                                                                                     ccATGGCGGCATCGGCGGCGTGCACGTTGCGGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgcttg                 >  1:50930/1‑150 (MQ=255)
                                                                                                                                     ccATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgcttg                 >  2:1548693/1‑150 (MQ=255)
                                                                                                                                     ccATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgcttg                 >  1:759947/1‑150 (MQ=255)
                                                                                                                                      cATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgcttgc                <  2:532885/150‑1 (MQ=255)
                                                                                                                                        tGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTgcgc              >  1:1169193/1‑150 (MQ=255)
                                                                                                                                              cATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAg        <  1:136787/150‑1 (MQ=255)
                                                                                                                                              cATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAg        <  2:1642704/150‑1 (MQ=255)
                                                                                                                                              cATCGGCGGCGTGCACGCGGCGTACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCATGTTCTGTTGCTTGCGCGCATAg        <  2:951947/150‑1 (MQ=255)
                                                                                                                                               aTCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAgg       >  2:285816/1‑150 (MQ=255)
                                                                                                                                                tCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGtt                                 >  2:578823/1‑123 (MQ=255)
                                                                                                                                                tCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGtt                                 <  1:578823/123‑1 (MQ=255)
                                                                                                                                                tCGGCAGCGTGCACGCCGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGc      >  2:1804145/1‑150 (MQ=255)
                                                                                                                                                   gcggcgTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCAc   >  1:1143760/1‑150 (MQ=255)
                                                                                                                                                    cggcgTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTATTGCTTGCGCGCATAGGCCACt  <  2:500822/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACT  >  NC_002947/3711794‑3712091

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